CN106150484B - A kind of drilling well leakage prediction technique based on digital rock mass - Google Patents
A kind of drilling well leakage prediction technique based on digital rock mass Download PDFInfo
- Publication number
- CN106150484B CN106150484B CN201510163437.1A CN201510163437A CN106150484B CN 106150484 B CN106150484 B CN 106150484B CN 201510163437 A CN201510163437 A CN 201510163437A CN 106150484 B CN106150484 B CN 106150484B
- Authority
- CN
- China
- Prior art keywords
- density
- parameter
- confidence level
- drilling
- pressure
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The drilling well leakage prediction technique based on digital rock mass that the present invention provides a kind of, this method comprises: the borehole track data of current well are acquired, to obtain the coordinate at borehole track each position；According to the coordinate in borehole track at each position, from the geologic parameter information obtained in the digital rock mass of current well from each position；According to the coordinate in borehole track at each position, the engineering parameter information at each position is obtained；According to the geologic parameter information at each position and the drilling well leakage value-at-risk at the engineering parameter acquisition of information position.The method of the present invention can be based on related geologic parameter and drilling engineering parameter, leakage risk profile is carried out to the full well section of drilling well, not only it can be used to carry out the evaluation and optimization of Drilling Design scheme, but also can identify in time the leakage symptom of a trend during wellbore construction, to avoid the generation of leakage to greatest extent.
Description
Technical field
The invention belongs to oil gas well drilling risk profile technical field more particularly to a kind of drilling well wells based on digital rock mass
Leak prediction technique.
Background technique
Leakage is one of most common complexity event in deep & ultra-deep well drilling engineering.In order to avoid leakage and guarantor occurs
Protect oil-gas Layer, it is necessary to a possibility that leakage risk, occurs for each interval of the full well of look-ahead, carry out in advance drilling plan optimization or
It takes preventive measures, pays close attention to high risk well section in the construction process, or identify imminent leakage risk in time, have
Put arrow, it is ensured that drilling safety.
To solve the above problems, domestic and international experts and scholars have had carried out some research, the Forecasting recognition method of existing leakage,
Or live phenomenon is observed by sensor, or determine according to the various characteristics of leakage after leakage generation or right
The geologic feature for causing leakage to occur carries out fuzzy evaluation macroscopic view and obtains the potential leakage risk in somewhere, or utilizes to be predicted
The partial parameters of well predict leakage with the method that some " characteristic values " of offset well leakage match.
However, most prior art is the method judged afterwards, although can control risk, partial loss is retrieved,
It is to inevitably lead to drilling fluid leakage and hydrocarbon contamination in this way.And the method predicted in advance does not account for specific drilling well
Interaction relationship between engineering design information and geological condition, thus can not be taken for specific well specific aim measure come
The prevention optimized.
Summary of the invention
The first technical problem to be solved by the present invention is to need to provide a kind of drilling well leakage prediction based on digital rock mass
Method, this method can carry out the potential leakage risk of the full well section of drilling well in conjunction with geology parameter information and engineering parameter information pre-
It surveys.
In order to solve the above-mentioned technical problem, it is pre- to provide a kind of drilling well leakage based on digital rock mass for embodiments herein
Survey method, this method comprises: the borehole track data of current well are acquired, to obtain the coordinate at borehole track each position；Root
According to the coordinate at each position in borehole track, from the geologic parameter letter obtained in the digital rock mass of current well from each position
Breath；According to the coordinate in borehole track at each position, the engineering parameter information at each position is obtained；At each position
Geologic parameter information and the engineering parameter acquisition of information position at drilling well leakage value-at-risk.
Preferably, in the step of geologic parameter information at each position is obtained in the digital rock mass from current well, packet
It includes: from the numerical value for extracting the geologic parameter from each position in the digital rock mass of current well；According to the data source of digital rock mass
And building process, obtain the confidence value of the geologic parameter at each position, wherein the geologic parameter includes the following: solution cavity
Information, fault information, crack information, rock type, permeability, formation pore pressure equal yield density and formation fracture pressure equivalent
Density.
Preferably, in the step of obtaining the engineering parameter information at each position, comprising: from drilling engineering design scheme
Or with the numerical value for extracting the engineering parameter at each position in the data of drilling collection；Every kind of engineering parameter is determined according to data source
Confidence level, wherein the engineering parameter includes the following: drilling fluid density, drilling fluid plastic viscosity, discharge capacity, drilling speed, landwaste ruler
Very little, landwaste, borehole diameter, drilling tool internal diameter, drilling tool outer diameter, hydrophthalmia diameter and surge pressure coefficient or surge pressure equivalent are close
Degree；The equivalent circulating density information at each position is calculated in engineering parameter information at each position based on extraction.
Preferably, in the pressure equivalent circulating density information being calculated at each position the step of, comprising: according to every
Engineering parameter information at one position selects respective algorithms to calculate the pressure equivalent circulating density at each position；According to meter
The parameter calculated algorithm used by pressure equivalent circulating density and participate in calculating obtains the confidence level of pressure equivalent circulating density.
Preferably, pressure is obtained in the algorithm according to used by calculating pressure equivalent circulating density and the parameter for participating in calculating
In the step of confidence level of equivalent circulating density ECD, further comprise: assuming that participating in the n parameter E1 ... En that have calculated, calculating
The independent confidence level CF (Ei) that each parameter calculates ECD=CF (E) [i] * CF (ECD algorithm), CF (ECD algorithm) indicate meter
Calculate the confidence level of the algorithm of ECD, the confidence level of CF (E) [i] expression parameter Ei；Calculate the synthetic reliability CF (E1E2) of E1 and E2
CF (E1E2) is then used as independent confidence level, it is counted with CF (E3) by=CF (E1)+CF (E2)-CF (E1) * CF (E2)
Calculation obtains CF (E1E2E3), and successively iterative calculation is up to the last one parameter En, and then obtains final synthetic reliability CF
(ECD), as the confidence level of pressure equivalent circulating density ECD.
Preferably, according to the geologic parameter information at each position and the drilling well at the engineering parameter acquisition of information position
In the step of leakage value-at-risk, comprising: according at each position geologic parameter information and engineering parameter information judge whether it is full
Sufficient preset condition: if satisfied, then determining initial risk values according to preset condition, and the confidence level of the initial risk values, base are calculated
Final risk value is obtained in initial risk values and its confidence level；If being unsatisfactory for preset condition, it is determined that final risk value 0.
Preferably, the preset condition includes at least one of: there are solution cavity or there are tomographies；Pressure equivalent circulation is close
Degree is greater than formation fracture pressure equal yield density with the sum of surge pressure coefficient or surge pressure equal yield density；Crack is to send out very much
It educates, and the sum of pressure equivalent circulating density and surge pressure coefficient or surge pressure equal yield density are worked as greater than formation pore pressure
Metric density；Crack is general development, and pressure equivalent circulating density and surge pressure coefficient or surge pressure equal yield density it
Be greater than formation pore pressure equal yield density；Rock type is sandstone, conglomerate or glutenite, and permeability is greater than certain value, pressure
Power equivalent circulating density and the sum of surge pressure coefficient or surge pressure equal yield density are greater than formation pore pressure equal yield density；
Rock type is sandstone, conglomerate or glutenite, and pressure equivalent circulating density and surge pressure coefficient or surge pressure equivalent
The sum of density is greater than formation pore pressure equal yield density；Rock type is sandstone, conglomerate or glutenite, and pressure equivalent circulation is close
Degree is more than or equal to formation pore pressure equal yield density；Pressure equivalent circulating density and surge pressure coefficient or surge pressure equivalent
The sum of density is greater than formation pore pressure equal yield density；Pressure equivalent circulating density is greater than formation pore pressure equal yield density.
Preferably, the confidence level of initial risk values is calculated according to the confidence level of parameter each in the preset condition met.
Preferably, further includes: geologic parameter information, engineering parameter information and drilling well leakage value-at-risk are based on, with same
Leakage risk profile section is shown with curve form on the basis of well depth.
Preferably, further includes: when the drilling well leakage value-at-risk at a position is more than or equal to setting thresholding, at the position
Engineering parameter information be adjusted so that the drilling well leakage value-at-risk at the position be less than setting thresholding.
Compared with prior art, one or more embodiments in above scheme can have following advantage or beneficial to effect
Fruit.
Drilling well leakage prediction technique based on digital rock mass of the invention, can be based on related geologic parameter and drilling engineering
Parameter carries out leakage risk profile to the full well section of drilling well, had not only been embodied leakage severity but also had embodied the quantization wind of probability of happening
Dangerous rank can not only be used to carry out the evaluation and optimization of Drilling Design scheme, but also can identify in time well during wellbore construction
The symptom of a trend is leaked, to avoid the generation of leakage to greatest extent.
Other features and advantages of the present invention will be illustrated in the following description, also, partly becomes from specification
It obtains it is clear that being understood by implementing technical solution of the present invention.The objectives and other advantages of the invention can by
Specifically noted structure and/or process are achieved and obtained in specification, claims and attached drawing.
Detailed description of the invention
Attached drawing is used to provide to the technical solution of the application or further understanding for the prior art, and constitutes specification
A part.Wherein, the attached drawing for expressing the embodiment of the present application is used to explain the technical side of the application together with embodiments herein
Case, but do not constitute the limitation to technical scheme.
Fig. 1 is the flow diagram of the drilling well leakage prediction technique based on digital rock mass of the embodiment of the present invention one.
Fig. 2 is the flow diagram of the drilling well leakage prediction technique based on digital rock mass of an example of the invention.
Specific embodiment
Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings and examples, how to apply to the present invention whereby
Technological means solves technical problem, and the realization process for reaching relevant art effect can fully understand and implement.This Shen
Please each feature in embodiment and embodiment, can be combined with each other under the premise of not colliding, be formed by technical solution
It is within the scope of the present invention.
In addition, the process of attached drawing can be in the computer system of such as a group of computer-executable instructions the step of illustrating
Middle execution.Also, although logical order is shown in flow charts, and it in some cases, can be to be different from herein
Sequence executes shown or described step.
(embodiment one)
Fig. 1 is the flow diagram of the drilling well leakage prediction technique based on digital rock mass of the embodiment of the present invention, is joined below
Examine each step that Fig. 1 carrys out process in detail.
In step s 110, the borehole track data of current well are acquired, to obtain the coordinate at borehole track each position.
Borehole track data include the information such as well depth, hole angle, azimuth.Generally, it is preferably taken a bit every relatively closely spaced,
1 point is taken every 1 meter in this example, successively acquires the data of each point, is then converted to the data of each point by calculating
Cartesian coordinate value X, Y, Z, finally obtain the coordinate sequence of a borehole track in relation to current well.
In the step s 120, it according to the coordinate in borehole track at each position, is obtained from the digital rock mass of current well
Geologic parameter information (including geologic parameter value and its confidence value) at each position.
It should be noted that digital rock mass refers to serial three-dimensional mesh data body, each grid node have coordinate value and
Attribute value (such as solution cavity, crack, porosity, pore pressure) with rock mass.Generally, it is based on existing geological structure mould
Type, it will be able to construct the digital rock mass for drilling well.Moreover, obtaining the various engineerings along well track from digital rock mass
Matter parameter, can fine optimization wellbore construction scheme, improve prediction drilling well leakage precision.
Specifically, in the step of geologic parameter information at each position is obtained in the digital rock mass from current well, packet
It includes: from the numerical value for extracting the geologic parameter from each position in the digital rock mass of current well, according to the data source of digital rock mass
And building process, obtain the confidence value of the geologic parameter at each position.Wherein, 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
Deng.
In step s 130, according to the coordinate in borehole track at each position, the engineering parameter at each position is obtained
Information (including engineering parameter values and its confidence value).
In particular it is required that executing following steps: being extracted in the data from drilling engineering design scheme or with drilling collection each
Engineering parameter at position determines the confidence level of every kind of engineering parameter according to data source.Wherein, 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 outer diameter,
The surge pressure equal yield density etc. that the surge pressure coefficient or calculating of hydrophthalmia diameter and setting obtain.Then, according to extraction
Engineering parameter information at each position is calculated the pressure equivalent circulating density at each position and (also joins for an engineering
Number) information.
In the pressure equivalent circulating density information being calculated at each position the step of, specifically, according to each
The engineering parameter information at place is set, selects respective algorithms to calculate the pressure equivalent circulating density at each position, then according to meter
The parameter calculated algorithm used by pressure equivalent circulating density and participate in calculating obtains the confidence level of pressure equivalent circulating density.
It is close that pressure equivalent circulation is obtained in the algorithm according to used by calculating pressure equivalent circulating density and the parameter used
In the step of spending the confidence level of ECD, further comprise: assuming that participate in calculating has n parameter E1 ... En, then first calculating each ginseng
Independent confidence level CF (Ei)=CF (E) [i] * CF (ECD algorithm) that number calculates ECD, CF (ECD algorithm) indicate to calculate ECD
Algorithm confidence level, the confidence level of CF (E) [i] expression parameter Ei.Then the synthetic reliability CF (E1E2) of E1 and E2 is calculated
=CF (E1)+CF (E2)-CF (E1) * CF (E2), then by CF (E1E2) as independent confidence level with similar method and CF
(E3) it is calculated CF (E1E2E3), i.e. CF (E1E2E3)=CF (E1E2)+CF (E3)-CF (E1E2) * CF (E3) ... ... is successively
Iterative calculation obtains final synthetic reliability CF (ECD), as pressure equivalent up to the position of the last one parameter En
The confidence level of circulating density ECD.
It should be noted that existing method assumes that the confidence level of used data is 100% accurate, it is used
It calculates and judgment method is also 100% reliable, however these data and method non-ideal in actual operation, compare existing skill
Art, the present invention consider the reliability of data and calculation method, further improve the accuracy of prediction drilling well leakage.
In step S140, according to the geologic parameter information and engineering parameter acquisition of information at each position of borehole track
Drilling well leakage value-at-risk at the position.
Specifically, according at each position geologic parameter information and engineering parameter information judge whether to meet default item
Part calculates the confidence level of the initial risk values if satisfied, then determining initial risk values according to preset condition, is based on initial risks
Value and its confidence level obtain final risk value.If being unsatisfactory for preset condition, it is determined that final risk value 0.In this example, according to
The confidence level of each parameter calculates the confidence levels of initial risk values in the preset condition met.
Further, above-mentioned preset condition may include at least one of: (1) there are solution cavity or there are tomographies；(2) pressure
Equivalent circulating density and the sum of surge pressure coefficient or surge pressure equal yield density are greater than formation fracture pressure equal yield density；
(3) crack is to develop very much, and pressure equivalent circulating density and the sum of surge pressure coefficient or surge pressure equal yield density are big
In formation pore pressure equal yield density；(4) crack is general development, and pressure equivalent circulating density and surge pressure coefficient or
The sum of surge pressure equal yield density is greater than formation pore pressure equal yield density；(5) rock type be sandstone, conglomerate or glutenite,
And permeability is greater than certain value, pressure equivalent circulating density is greater than with the sum of surge pressure coefficient or surge pressure equal yield density
Formation pore pressure equal yield density；(6) rock type be sandstone, conglomerate or glutenite, and pressure equivalent circulating density and excitement
The sum of pressure coefficient or surge pressure equal yield density are greater than formation pore pressure equal yield density；(7) rock type is sandstone, gravel
Rock or glutenite, and pressure equivalent circulating density is more than or equal to formation pore pressure equal yield density；(8) pressure equivalent circulating density
It is greater than formation pore pressure equal yield density with the sum of surge pressure coefficient or surge pressure equal yield density；(9) pressure equivalent follows
Ring density is greater than formation pore pressure equal yield density.
Since the present invention not only allows for initial risk values when determining value-at-risk, it is also contemplated that initial risk values it is credible
Degree, therefore the accuracy of obtained drilling well leakage value-at-risk is higher.Although the method judged in the prior art can be use up afterwards
Amount is controlled risk, and retrieves partial loss, but inevitably lead to drilling fluid leakage and hydrocarbon contamination in this way, and if adopt
It takes measure not in time, is easy to deteriorate the state of affairs, lead to a disaster.And the present invention can accurately predict leakage at design initial stage
The risk of generation can early avoid the generation of risk.
In addition, as shown in Figure 1, this method can also include step S150.
In step S150, it is based on geologic parameter information, engineering parameter information and drilling well leakage value-at-risk, with the same well
Leakage risk profile section is shown with curve form on the basis of deep axis.It can be shown in the leakage risk profile sectional view
Geologic parameter curve, engineering parameter curve and the risk class curve of the borehole track of the well.If being provided with alarm threshold, that
The leakage risk profile section can also show risk alarm prompt information.
When drilling well leakage value-at-risk at a position is more than or equal to setting thresholding, to the engineering parameter information at the position
It is adjusted, so that the drilling well leakage value-at-risk at the position is less than setting thresholding.
By the step can support technician check the geology, project situation and corresponding risk of full well section comprehensively
The high well section of rank, especially risk class formulates improvement or counter-measure, can control wind by further analyzing reason
The generation of danger.
Fig. 2 is the flow diagram of the drilling well leakage prediction technique based on digital rock mass of an example of the invention.The example
It is that leakage risk profile is carried out to Drilling Design scheme, is described in detail below with reference to Fig. 2 and obtains drilling well leakage value-at-risk
Method.
In step 1, since 0m, using well depth 1m as interval, the seat of each point is extracted from wellbore designed path data
Scale value (well depth, northern coordinate, eastern coordinate), obtains a coordinate array T [N] [3] (assuming that one shares N number of coordinate points).
In step 2, the coordinate array obtained according to step 1 extracts the geology ginseng of each coordinate points from digital rock mass
Number data, form geologic parameter array G [N] [I] (I=7), wherein i respectively represents solution cavity (with or without) from 0 to 6, tomography (has
Or nothing), crack (very development, general development, nothing), rock type, permeability K, formation pore pressure equal yield density rp, stratum
Fracture pressure equal yield density rf.Meanwhile according to the primary data source of digital rock mass and building process, it is each to give each depth
The confidence level CF (G) [n] [i] of attribute, obtains confidence level array CF (G) [N] [I] (I=7).
It should be noted that digital rock mass is mainly based upon the historical data of offset well to construct, therefore offset well quantity, neighbour
The average distance of well and this well, offset well the data precision all will affect the confidence level of digital rock mass.
Generally, the independent confidence level of three can be set according to expertise and certain rule: 1. constructing digital rock mass
The confidence level of offset well the data precision (comprehensively considers the original accuracy of these offset well all datas, calculation method by expert
Maturity and give), for example the confidence level of block permeability K is defined as CF (K0).2. offset well quantity confidence level CF (W) can root
Degree of membership delimited according to experience, set if W < 3 CF (W)=0.4,3≤W < 10 then sets CF (W)=0.5, and 10≤W < 20 then sets CF
(W)=0.6,20≤W < 50 is then set CF (W)=0.7, and 50≤W < 200 then sets CF (W)=0.8, and 200≤W < 800 then sets CF
(W)=0.9, W > 800 item set CF (W)=1.3. average distance confidence level CF (D) can rule of thumb delimit degree of membership, such as D < 100m
It then sets CF (D)=1,100≤D < 2000 is then set CF (D)=0.9, and 2000≤D < 5000 then sets CF (D)=0.8,5000≤D <
10000 set CF (D)=0.7, and 10000≤D < 20000 is then set CF (D)=0.6,20000≤D < 50000 then set CF (D)=
0.5,50000≤D < 100000 then sets CF (D)=0.4 ..., finally calculates geometrical mean (or the weighted geometric mean of three
Value) obtain the confidence level CF (K) of the point.The different interval divisions that are subordinate to will affect the size of CF (K) value, but not influence its phase
To changing rule, therefore it will not influence its engineering application value.
In step 3, the coordinate array obtained according to step 1, from drilling engineering design scheme or with drilling collection
It extracts the engineering parameter of each coordinate points in data, is formed engineering parameter array E [N] [I] (I=11), wherein i is from 0 to 10 point
Drilling fluid density, drilling fluid plastic viscosity, design discharge capacity, design drilling speed, landwaste size, landwaste, borehole diameter, drilling tool are not represented
Internal diameter, drilling tool outer diameter, hydrophthalmia diameter, the surge pressure coefficient S of settingg.Every kind of engineering parameter is given according to data source situation
Confidence value CF (E) [i].
If these data are from engineering design, therefore, it is considered that all parameters are completely credible, set each parameter can
Certainty value CF (E) [i]=1；If from brill certain parameters can be reset according to the reliability of collection in worksite data
The value of CF (E) [i].
It should be noted that data source mode is identical, therefore for same for the design of a well or real brill
A parameter does not have difference in the confidence level of different depth, therefore ignores [N] in array here, i.e., CF is all used in any depth point
(E0) [I] this group of confidence level.
In step 4, it based on the data in step 3, (can also be used using the algorithm that SY/T 6613-2005 standard provides
Other similar algorithm), the pressure equivalent circulating density (i.e. ECD) of each coordinate points is calculated, r is usedcIt indicates, forms ECD array rc
[N]。
If using professional standard algorithm, then it is assumed that the algorithm is very credible, and the confidence level that can set the algorithm (is advised
Then intensity) CF (ECD algorithm)=0.98；If using other algorithms corresponding CF can be given according to the degree of recognition to the algorithm
Value, CF ∈ [0,1].
For each engineering parameter that this algorithm is used, the confidence level of final ECD is calculated as follows:
Assuming that participate in calculating has n parameter E1 ... En, then the independent confidence level that each parameter calculates ECD is first calculated
CF (Ei)=CF (E) [i] * CF (ECD algorithm), CF (ECD algorithm) indicate to calculate the confidence level of the algorithm of ECD, CF (E) [i] table
Show the confidence level of parameter Ei.Then synthetic reliability CF (E1E2)=CF (E1)+CF (E2)-CF (E1) * CF of E1 and E2 is calculated
(E2), CF (E1E2E3) then is calculated with similar method and CF (E3) as independent confidence level in CF (E1E2), successively
Iterative calculation obtains final summation confidence level CF (ECD) until the last one parameter position.Compared to other methods, in utilization
The confidence level precision that the method for stating calculates is higher.
In steps of 5, the prediction for starting leakage risk, enables n=0.
In step 6, according to n value, data G [n] is extracted from geologic parameter array G [N] [I] respectively, from ECD array rc
R is extracted in [N]c[n] extracts array CF (G) [n] from geologic parameter confidence level array CF (G) [N] [I], from engineering parameter number
Array E [n] is extracted in group E [N] [I].
In step 7, if (solution cavity=have) or (tomography=have), initial risk values R0=5, it is then based on solution cavity and disconnected
The confidence level (i.e. CF (solution cavity) and CF (tomography)) of layer calculates the confidence level CF (R0), final risk value R=R0*CF of R0
(R0), step 17 is then executed.Wherein, CF (R0)=max (CF (solution cavity), CF (tomography)).
In step 8, if rc+Sg>rf, then R0=5, is then based on rc、Sg、rfConfidence level calculate the confidence level CF of R0
(R0), final risk value R=R0*CF (R0) then executes step 17.Due to S hereingFor design value, therefore it is not related to it
The problem of confidence level, then CF (R0)=CF (rc)+CF(rf)-CF(rc)*CF(rf).In other examples, if the lower brill bored in fact
Process, then SgIt, at this time need to be according to the iterative calculation r of step 4 for the surge pressure equal yield density actually calculatedc、Sg、rfThis three
The synthetic reliability of a factor.
In step 9, if (crack=develop very much) and (rc+Sg>rp), then R0=5, is then based on crack, rc、Sg、rf
Confidence level calculate confidence level CF (R0), the final risk value R=R0*CF (R0) of R0, then execute step 17.Here R0
Confidence level can iterate to calculate crack, r according to the method for step 4c、rfThe synthetic reliability of these three factors.In real drill-through journey
Under lower brill state, S is substituted with surge pressure equal yield densityg, its confidence level is considered at this time, and the confidence level of R0 can be according to step 4
Method iterates to calculate crack, rc、Sg、rfThe synthetic reliability of this four factors.
In step 10, if (crack=general development) and (rc+Sg>rp), then R0=4, is then based on crack, rc、Sg、rf
Confidence level calculate confidence level CF (R0), the final risk value R=R0*CF (R0) of R0, then execute step 17.R0's is credible
The circular for spending CF (R0) can be identical as 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 is then based on rock type, K, rc、Sg、rfConfidence level calculate the confidence level CF (R0), final risk value R=of R0
R0*CF (R0) then executes step 17.The circular of the confidence level CF (R0) of R0 can be similar with step 10.
By the step can be seen that this example not only allow for fracture hole leakage, but also can to permeability loss into
Row identification.
In step 12, if (rock type=sandstone or conglomerate or glutenite) and (rc+Sg>rp), then R0=3, then base
In rock type, rc、Sg、rfConfidence level calculate confidence level CF (R0), the final risk value R=R0*CF (R0) of R0, then
Execute step 17.The circular of the confidence level CF (R0) of R0 can be similar with step 10.
In step 13, if (rock type=sandstone or conglomerate or glutenite) and rc≥rp, then R0=2, is then based on rock
Stone type, rc、rfConfidence level calculate the confidence level CF (R0) of R0, final risk value R=R0*CF (R0) then executes step
17.Here the confidence level of R0 can iterate to calculate rock type, r according to the method for step 4c、rfThe synthesis of these three factors is credible
Degree.
At step 14, if rc+Sg>rp, then R0=2, is then based on rc、Sg、rfConfidence level calculate the confidence level of R0
CF (R0), final risk value R=R0*CF (R0)；Execute step 17.Due to S hereingFor design value, therefore not being related to it can
The problem of reliability, then CF (R0)=CF (rc)+CF(rf)-CF(rc)*CF(rf).If the going down process bored in fact, then SgFor reality
The surge pressure equal yield density of calculating, at this time need to be according to the iterative calculation r of step 4c、Sg、rfThe synthesis of these three factors is credible
Degree.
In step 15, if rc>rp, then R0=1, is then based on rc、rfConfidence level calculate the confidence level CF of R0
(R0), final risk value R=R0*CF (R0) then executes step 17.Wherein, CF (R0)=max (CF (rc),CF(rf))。
In step 16, if above-mentioned condition is all unsatisfactory for, R=0.
In step 17, R [n]=R is enabled.
In step 18, if n < N, n=n+1,6~17 is returned to step, is otherwise terminated.
Confidence level CF (R0) in addition to calculating R0 using above-mentioned alternative manner can also utilize the confidence level of each factor
Geometrical mean obtains the confidence level CF (R0) of R0, is not limited thereto.
It should be noted that if being for brill, engineering parameter to be real brill data rather than design data, part geologic parameter
Information also can be by measuring and calculating is updated.During boring in fact, when in " lower to bore " state, above-mentioned SgUnified reality
The surge pressure equal yield density of calculating substitutes.
The present invention is suitable for drilling engineering design personnel, field technician, rear decision expert and production management personnel
The potential leakage risk profile of full well section, and real-time monitoring risk symptom of a trend during boring in fact are carried out, to reduce construction risk.
Example 1
In the YB*** well Drilling Design stage, according to conventional design cycle, by borehole track, casing programme, drill set
Conjunction, drill bit, Hydraulics design complete (well depth 0m~6756m).The engineering needed for extracting this method in this design scheme
Parameter set, while from the digital rock mass of YB block where the well, the rock mass attribute value that the borehole track is passed through is extracted, i.e.,
Matter parameter set.Above-mentioned engineering parameter collection and geologic parameter collection are respectively with the extraction of identical depth interval, such as 1m.Then this is utilized
The invention method differentiates to obtain the risk class of each point on borehole track, finds in 4801m to 5112m well section leakage wind
Dangerous rank R > 4, designer think that this potential risk is too high, then pass through and carefully study, and repeatedly adjust relevant design parameter,
Risk profile is re-started every time, until value-at-risk R < 3, reach expected design, it is just final to finalize a text.
This example considers the interaction relationship between specific drilling engineering design information and geological condition, therefore can be with
Specific aim measure is taken to optimize and prevent for specific well.
Example 2
During SB*** well wellbore construction, using the well track information and mud logging technique of real-time measurement, obtain in real time
Engineering parameter collection, and ECD and surge pressure are calculated in real time；Where from the well in the digital rock mass of block, the reality drilling well is extracted
The rock mass attribute value that eye track passes through or passing through, while formation pore pressure is calculated using the engineering parameter of real-time measurement
To correct the value formation pressure that former digital rock mass extracts.Then before differentiating to obtain well track most using method of the present invention
The leakage risk class of end or drill bit position or entire well section.One day, in the going down process that makes up a joint, Field Force's setting
Good running speed, starts lower brill, is judged by this method, shows the leakage value-at-risk R=of current bottom hole location (5263.2m)
4.6, indicate risk it is very high, technical staff think may because under it is drill-through cause pressure oscillation fastly caused by, slow down running speed immediately,
More new data differentiates again, and risk class is down to 3 hereinafter, effectively preventing malignant event.
The present invention devises a kind of drilling well leakage prediction technique based on digital rock mass, the geology extracted from digital rock mass
Parameter, and drilling engineering parameter coupling analysis, complete the potential leakage risk profile of well section or local well section is carried out before brill, with brill
Identification early warning in real time is carried out to the symptom of a trend of leakage risk in the process, energy aided drilling designer and construction technical staff control are bored
Well risk.
The present invention can be used for the Drilling Design stage and carry out leakage prediction, and then optimization design to full well, and propose reply wind
The prediction scheme of danger；In the construction process, automatic amendment and early warning in real time are carried out to the risk of prediction before drilling, can utmostly controls well
The generation of leakage finds risk in first time and takes the measure of being effectively treated, and loss is minimized.
Although disclosed herein embodiment it is as above, the content only technical solution of the present invention for ease of understanding
And the embodiment used, it is not intended to limit the invention.Technical staff in any fields of the present invention is not departing from this
Under the premise of the disclosed spirit and scope of invention, any modification and variation can be carried out in the form and details of implementation,
But scope of patent protection of the invention, still should be subject to the scope of the claims as defined in the appended claims.
Those of ordinary skill in the art will appreciate that implement the method for the above embodiments be can be with
Relevant hardware is instructed to complete by program, the program can be stored in a computer readable storage medium,
The program when being executed, includes the steps that method in the embodiment, the storage medium, such as: ROM/RAM, magnetic disk, CD
Deng.
Claims (7)
1. a kind of drilling well leakage prediction technique based on digital rock mass, this method comprises:
The borehole track data of current well are acquired, to obtain the coordinate at borehole track each position；
According to the coordinate in borehole track at each position, from the geology ginseng obtained in the digital rock mass of current well from each position
Number information；
According to the coordinate in borehole track at each position, the engineering parameter information at each position is obtained；
According to the geologic parameter information at each position and the drilling well leakage value-at-risk at the engineering parameter acquisition of information position；
Wherein, in the step of obtaining the engineering parameter information at each position, comprising:
The numerical value of the engineering parameter from each position is extracted in data from drilling engineering design scheme or with drilling collection；
The confidence level of every kind of engineering parameter is determined according to data source, wherein the engineering parameter includes the following: that drilling fluid is close
Degree, drilling fluid plastic viscosity, discharge capacity, drilling speed, landwaste size, landwaste, borehole diameter, drilling tool internal diameter, drilling tool outer diameter, hydrophthalmia diameter
With surge pressure coefficient or surge pressure equal yield density；
The equivalent circulating density information at each position is calculated in engineering parameter information at each position based on extraction；
Wherein, in the pressure equivalent circulating density information being calculated at each position the step of, comprising:
According to the engineering parameter information at each position, it is close to calculate pressure equivalent circulation at each position to select respective algorithms
Degree；
The algorithm according to used by calculating pressure equivalent circulating density and the parameter for participating in calculating obtain pressure equivalent circulating density
Confidence level；
Wherein, it obtains pressure equivalent in the algorithm according to used by calculating pressure equivalent circulating density and the parameter for participating in calculating and follows
In the step of confidence level of ring density ECD, further comprise:
Assuming that participate in calculating has n parameter E1 ... En, calculate independent confidence level CF (Ei) that each parameter calculates ECD=
CF (E) [i] * CF (ECD algorithm), CF (ECD algorithm) indicate to calculate the confidence level of the algorithm of ECD, CF (E) [i] expression parameter Ei
Confidence level；
Synthetic reliability CF (E1E2)=CF (E1)+CF (E2)-CF (E1) * CF (E2) for calculating E1 and E2, then by CF
(E1E2) it is used as independent confidence level, it is carried out CF (E1E2 E3) is calculated with CF (E3), successively iterative calculation is until last
One parameter En, and then final synthetic reliability CF (ECD) is obtained, as the credible of pressure equivalent circulating density ECD
Degree.
2. the method according to claim 1, wherein being obtained at each position in the digital rock mass from current well
Geologic parameter information the step of in, comprising:
From the numerical value for extracting the geologic parameter from each position in the digital rock mass of current well；
According to the data source and building process of digital rock mass, the confidence value of the geologic parameter at each position is obtained,
Wherein, the geologic parameter includes the following: solution cavity information, fault information, crack information, rock type, permeability, stratum
Pore pressure equal yield density and formation fracture pressure equal yield density.
3. method described in any one of -2 according to claim 1, which is characterized in that according to the geologic parameter at each position
In the step of drilling well leakage value-at-risk at information and the engineering parameter acquisition of information position, comprising:
According at each position geologic parameter information and engineering parameter information judge whether to meet preset condition:
If satisfied, then determining initial risk values according to preset condition, and the confidence level of the initial risk values is calculated, is based on initial wind
Danger value and its confidence level obtain final risk value；If being unsatisfactory for preset condition, it is determined that final risk value 0.
4. according to the method described in claim 3, it is characterized in that, the preset condition includes at least one of:
There are solution cavity or there are tomographies；
The sum of pressure equivalent circulating density and surge pressure coefficient or surge pressure equal yield density are worked as greater than formation fracture pressure
Metric density；
Crack is to develop very much, and pressure equivalent circulating density and the sum of surge pressure coefficient or surge pressure equal yield density are big
In formation pore pressure equal yield density；
Crack is general development, and pressure equivalent circulating density and the sum of surge pressure coefficient or surge pressure equal yield density are big
In formation pore pressure equal yield density；
Rock type is sandstone, conglomerate or glutenite, and permeability is greater than certain value, pressure equivalent circulating density and surge pressure
The sum of coefficient or surge pressure equal yield density are greater than formation pore pressure equal yield density；
Rock type is sandstone, conglomerate or glutenite, and pressure equivalent circulating density and surge pressure coefficient or surge pressure
The sum of equal yield density is greater than formation pore pressure equal yield density；
Rock type is sandstone, conglomerate or glutenite, and pressure equivalent circulating density is close more than or equal to formation pore pressure equivalent
Degree；
The sum of pressure equivalent circulating density and surge pressure coefficient or surge pressure equal yield density are worked as greater than formation pore pressure
Metric density；
Pressure equivalent circulating density is greater than formation pore pressure equal yield density.
5. according to the method described in claim 3, it is characterized in that,
The confidence level of initial risk values is calculated according to the confidence level of parameter each in the preset condition met.
6. the method according to claim 1, wherein further include:
Based on geologic parameter information, engineering parameter information and drilling well leakage value-at-risk, with curved shape on the basis of the same well depth
Formula shows leakage risk profile section.
7. method according to claim 1 or 6, which is characterized in that further include:
When drilling well leakage value-at-risk at a position is more than or equal to setting thresholding,
Engineering parameter information at the position is adjusted, so that the drilling well leakage value-at-risk at the position is less than setting door
Limit.
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 CN106150484A (en) | 2016-11-23 |
| CN106150484B true 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) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108694459B (en) * | 2017-04-10 | 2021-11-02 | 中国石油化工股份有限公司 | Drilling scheme optimization method based on virtual simulation |
| CN108804722B (en) * | 2017-04-26 | 2021-07-27 | 中国石油化工股份有限公司 | Parameter calculation method and device for drilling simulation |
| CN110656932B (en) * | 2019-10-24 | 2020-04-24 | 西南石油大学 | Drilling fluid loss position estimation method in drilling process |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101139925B (en) * | 2006-09-08 | 2012-07-04 | 西南石油大学 | Method for while-drilling testing reservoir parameter property and adjusting well drilling action in real time |
| CN101726761B (en) * | 2008-10-15 | 2012-03-07 | 中国石油天然气股份有限公司 | Risk-constrained oil and gas resource spatial distribution forecasting method |
| CN102128022B (en) * | 2010-12-30 | 2013-06-12 | 中国电子科技集团公司第二十二研究所 | Drilling engineering early warning method and system thereof |
| GB201204815D0 (en) * | 2012-03-19 | 2012-05-02 | Halliburton Energy Serv Inc | 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 |
-
2015
- 2015-04-08 CN CN201510163437.1A patent/CN106150484B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106150484A (en) | 2016-11-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106150485B (en) | A kind of drilling well leakage forecasting system based on digital rock mass | |
| CN108388921B (en) | Overflow loss real-time identification method based on random forest | |
| CN104790926B (en) | A kind of fracture-pore reservoir waterflooding development effect evaluation method | |
| AU2013274606B2 (en) | Methods and related systems of building models and predicting operational outcomes of a drilling operation | |
| KR101904278B1 (en) | Method for decline curve analysis according to cumulative production incline rate in unconventional gas field | |
| BR112015009197B1 (en) | METHOD AND SYSTEM FOR PERFORMING A DRILLING OPERATION | |
| US10168447B2 (en) | Automatic geosteering and evolutionary algorithm for use with same | |
| US20150276978A1 (en) | System, Method and Computer Program Product For Determining Placement of Perforation Intervals Using Facies, Fluid Boundaries, Geobodies and Dynamic Fluid Properties | |
| CN107895092B (en) | Inter-well communication quantitative evaluation method based on complex nonlinear injection-production modeling | |
| CN106150484B (en) | A kind of drilling well leakage prediction technique based on digital rock mass | |
| Aljubran et al. | Deep learning and time-series analysis for the early detection of lost circulation incidents during drilling operations | |
| CN111612899B (en) | Carbonate fracture-cavity type oil reservoir geological modeling method | |
| CN105484724A (en) | Drilling downhole anomaly monitoring method | |
| CN108952657A (en) | A kind of horizontal well platform fracturing fracture length determining method and device | |
| CN107153881B (en) | Drilling risk prediction method | |
| Dekkers et al. | Resource assessment: estimating the potential of a geothermal reservoir | |
| Bahari et al. | Intelligent drilling rate predictor | |
| CN112836442A (en) | Method for determining hydraulic fracture injection amount of old well of shale oil horizontal well pattern | |
| Liang et al. | RESEARCH ON DRILLING KICK AND LOSS MONITORING METHOD BASED ON BAYESIAN CLASSIFICATION. | |
| Rafiei | Improved oil production and waterflood performance by water allocation management | |
| CN117150875A (en) | Pre-drilling logging curve prediction method based on deep learning | |
| Mouedden et al. | A screening methodology using fuzzy logic to improve the well stimulation candidate selection | |
| CN112241835A (en) | Deep shaft project water inrush disaster multi-source information evaluation method | |
| Artun et al. | A pattern-based approach to waterflood performance prediction using knowledge management tools and classical reservoir engineering forecasting methods | |
| Lamberghini et al. | Performance Based Reservoir Characterization in a Tight Gas Reservoir System-Case Study from Lajas and Punta Rosada Formations in the Neuquen Basin, Argentina |
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 |