CN110107282A - A kind of paper shale grease horizontal well is with boring method for controlling trajectory - Google Patents

A kind of paper shale grease horizontal well is with boring method for controlling trajectory Download PDF

Info

Publication number
CN110107282A
CN110107282A CN201910461192.9A CN201910461192A CN110107282A CN 110107282 A CN110107282 A CN 110107282A CN 201910461192 A CN201910461192 A CN 201910461192A CN 110107282 A CN110107282 A CN 110107282A
Authority
CN
China
Prior art keywords
well
dessert
time
value
curve
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
Application number
CN201910461192.9A
Other languages
Chinese (zh)
Other versions
CN110107282B (en
Inventor
周立宏
肖敦清
韩国猛
马建英
杨飞
范军侠
石倩茹
陈双清
周可佳
王昌丽
付东立
梁晨
尹玲玲
孟立娜
周素彦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Research Institute Of Petroleum Exploration & Development Dagang Oil Field Of Cnpc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Research Institute Of Petroleum Exploration & Development Dagang Oil Field Of Cnpc filed Critical Research Institute Of Petroleum Exploration & Development Dagang Oil Field Of Cnpc
Priority to CN201910461192.9A priority Critical patent/CN110107282B/en
Publication of CN110107282A publication Critical patent/CN110107282A/en
Application granted granted Critical
Publication of CN110107282B publication Critical patent/CN110107282B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • E21B47/12Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Remote Sensing (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a kind of paper shale grease horizontal wells with brill method for controlling trajectory, comprising the following steps: 1) determines optimal artificial synthesized seismic channel, window when generating earthquake dessert according to horizontal well dessert section seismic velocity;2) dessert segment standard log response curve is established using the corresponding form of fixed well log, and prediction log is generated according to shale oil horizontal well design track;3) in drilling course, according to Real-time Logging curve, reference mark interval well logging absolute value and amplitude of variation are adjusted in time with brill track, beneficial effects of the present invention: 1) in human intervention, depth relationship when establishing more accurate;2) when the basis established of depth relationship on, earthquake profile position where calculating dessert top bottom determines dessert time window, and by live fluoroscopic drill bit position in drilling process, optimization is with boring track;3) the area stratum standard log curve is established using fixed well and stratigraphic dip, adjusts track in time, improve reservoir-encountered rate.

Description

A kind of paper shale grease horizontal well is with boring method for controlling trajectory
Technical field:
The present invention relates to oil exploration technology fields, more particularly to a kind of paper shale grease horizontal well with brill TRAJECTORY CONTROL The technical field of method.
Background technique:
The vertical thickness of marine facies shale oil dessert generally between 20~90m, stablize by plane distribution, and attitude of stratum is stablized, and is fitted Heshui horizontal well scale development.Based on terrestrial lake basin, terrestrial facies shale oil dessert spatial distribution stability is opposite for China oil basin Poor, for vertical thickness generally between 6~15m, plane distribution has apparent mid-depth, dessert thinner thickness and variation Fastly, plane distribution stability is poor, continuous in dessert section to how accurately to control well track during shale oil horizontal well construction Drilling brings problem.In view of the above problems, the present invention sums up a kind of paper shale grease horizontal well with boring method for controlling trajectory, land Practicable technical method is provided with brill TRAJECTORY CONTROL during phase lake basin paper shale oil gas horizontal well construction.
Summary of the invention:
The object of the invention is that providing a kind of paper shale grease horizontal well with boring track to solve the above-mentioned problems Control method solves the problems, such as existing at present.
To solve the above-mentioned problems, the present invention provides a kind of technical solutions:
A kind of paper shale grease horizontal well is with boring method for controlling trajectory, comprising the following steps:
1) optimal artificial synthesized seismic channel is determined, window when generating earthquake dessert according to horizontal well dessert section seismic velocity;
2) dessert segment standard log response curve is established using the corresponding form of fixed well log, and according to shale grease Horizontal well planned course generates prediction log;
3) in drilling course, according to Real-time Logging curve, reference mark interval well logging absolute value and amplitude of variation are adjusted in time It is whole with bore track.
Preferably, optimal artificial synthesized seismic channel is determined in the step 1), according to horizontal well dessert section seismic velocity Window determines 2 parts when window mainly includes holostrome Duan Jingzhen calibration and horizontal well earthquake dessert when generating earthquake dessert.
Preferably, holostrome Duan Jingzhen calibration the following steps are included:
A., optimal artificial synthesized seismic channel is determined;
Individual well geological layering by well head to shaft bottom number consecutively be G1, G2, Gn, utilize acoustic travel time logging bent Line (AC) artificial synthesized seismic channel correspond to time depth be Gt1, Gt2, Gtn, seismic interpretation interface corresponds to time depth From top to bottom number consecutively be T1, T2, Tn;
It defines well and shakes degree of fitting DF formula:
Wherein, Gt1, Gt2, Gtn be that artificial synthetic seismogram corresponds to time depth, T1, T2, Tn be Seismic interpretation interface corresponds to time depth;Well shake degree of fitting df value is smaller, and it is higher to represent well shake fitting degree.
Mean value is denoted as avgAC in interval transit time curve measuring point depth bounds;
WhenWhen, avgAC is bigger than normal or less than normal, causes artificial synthesized earthquake result and earthquake reflective data It is not inconsistent, needs artificially to adjust avgAC value at this time, artificial synthesized record is corrected, reduces DF value according to this, it is quasi- to improve well shake Conjunction degree is higher;
WhenWhen, avgAC is bigger than normal, and the corresponding df value of avgAC is denoted as DF0, and interval transit time value is integrally subtracted Small fixed value a (generally taking 5~10) re-starts artificial synthesized earthquake using the interval transit time curve newly obtained and notes down, right The df value answered is denoted as df1.The step is repeated, until Successively determine DF2, DF3, DFn, take DF The corresponding artificial synthesized seismic channel of minimum value is optimal synthetic seismogram;
WhenWhen, avgAC is less than normal, and the corresponding DF value of avgAC is denoted as DF0, and interval transit time value is integrally increased Add a fixed value a (generally taking 5~10), re-starts artificial synthesized earthquake using the interval transit time curve newly obtained and note down, it is right The df value answered is denoted as DF1.The step is repeated, until Successively determine DF2, DF3, DFn, take DF The corresponding artificial synthesized seismic channel of minimum value is optimal artificial synthesized seismic channel.
B, the artificial synthesized seismic channel of layer series optimization is divided to determine;
For lineups included in the interface seismic interpretation T1 and T2 from the bottom to top respectively number be t1, t2, Tm, take-off depth number consecutively corresponding to individual well be g1, g2, gm, it is deep that optimal artificial synthesized seismic channel correspond to the time Degree be gt1, gt2, gtm;
It defines well and shakes degree of fitting df formula:
Wherein, gt1, gt2, gtn be that artificial synthetic seismogram corresponds to time depth, t1, t2, tn be Seismic interpretation interface corresponds to time depth.Well shake degree of fitting df value is smaller, and it is higher to represent well shake fitting degree;
Mean value is denoted as avgac in interval transit time curve measuring point depth bounds.
WhenWhen, avgac is bigger than normal or less than normal, causes artificial synthesized earthquake result and earthquake reflective data It is not inconsistent, needs artificially to adjust avgac value at this time, artificial synthesized record is corrected, reduces df value according to this, it is quasi- to improve well shake Conjunction degree is higher.
WhenWhen, avgAC is bigger than normal, and the corresponding df value of avgAC is denoted as df0, and interval transit time value is integrally subtracted Small fixed value a (generally taking 3~6) re-starts artificial synthesized earthquake using the interval transit time curve newly obtained and notes down, corresponding Df value be denoted as df1.The step is repeated, until Successively determine df2, df3, dfn, take df most It is small to be worth corresponding artificial synthesized seismic channel optimal synthetic seismogram between the interface T1 and T2.
WhenWhen, avgAC is less than normal, and the corresponding df value of avgAC is denoted as df0, and interval transit time value is integrally increased Add a fixed value a (generally taking 3~6), re-starts artificial synthesized earthquake using the interval transit time curve newly obtained and note down, it is corresponding Df value be denoted as df1.The step is repeated, until Successively determine df2, df3, dfn, take df most It is small to be worth corresponding artificial synthesized seismic channel optimal artificial synthesized seismic channel between the interface T1 and T2.
Repeat the above steps, successively optimize T2~T3, T3~T4, the axis of reflection between Tn-1~Tn so that most Excellent artificial synthesized seismic channel is in optimum state, and establishes time and depth transfer relationship.
Preferably, when the earthquake dessert window it is determining the following steps are included:
A steady reflection lineups are found in the corresponding seismic reflector of dessert body, by Event tracking, Dessert steady reflection interface t, seismic velocity v in dessert section depth bounds are found in three-dimensional space.On steady reflection interface Dessert body is denoted as Su, and thickness is denoted as du, and dessert body is denoted as Sd under steady reflection interface, and thickness is denoted as dd, and dessert time window is (t-du/v, t+dd/v).If steady reflection lineups are in the top of dessert body, dessert time window is (t, t+dd/v). If steady reflection lineups are in the bottom of dessert body, dessert time window is (t-du/v, t).
Wherein, in the step 2), dessert segment standard log response is established using the corresponding form of fixed well log Curve, and prediction log is generated according to shale oil horizontal well design track.
1. standard log curve calculates;
The geologic body limited by G1 and G2, O point are the intersection points of well 1 Yu geological interface G1, and Q point is well 1 and geological interface G2 Intersection point, geologic dip α, P are any points in drilling well, which is β, and P ' is that P point (hangs down in the projection of vertical direction It is deep), P " is projection of the P point on ground layer vertical.OP is known length section, and OP is ground layer vertical projected length (OP "=OP Sin α/cos β), for any point P, unique subpoint P " can be found in OQ ".OQ " is the true thickness of geologic body, referred to as For standard paragraphs, assigns log value measured by all measuring point P to OP ", both obtained geologic body standard log curve.
The step is repeated, the standard log curve A of all geologic bodies is obtained.
2. generating horizontal wellbore logging curve according to planned course;
The intersection point A of planned well and geological interface G1, the intersection point with geological interface G2 are B, and P is any designed on well location Point, the point design hole angle are β, and P point point corresponding with standard log curve is P ", and geologic body inclination angle is that the α hole angle is β then has AP=OQ " cos β/sin α), for any point P, unique subpoint P " can be found therewith in standard paragraphs OQ " It is corresponding.By establishing the projection relation between OQ " and AB, the corresponding standard log curve values of P " point are assigned to P point Planned well is obtained in the intracorporal prediction log of the geology.
The step is repeated, prediction log is carried out to geologic body in all design well sections and is calculated, and then obtains all set It counts the full well section of well and predicts log B.
Preferably, in drilling course, the relationship between measured curve and prediction curve is compared in the step 3), and According to time and depth transfer relationship, optimization in time is as follows with Trace step is bored:
(1) well logging curve monitors in real time: real time contrast surveys log and predicts the relationship between log, According to the distance between measuring instrument and drill bit, drill bit position is grasped in time;
(2) earthquake trajectory real-time monitoring: utilizing time and depth transfer relationship, wellbore trace projected in seismic data body, and According to brill track position relative position in dessert body;
(3) optimize well track in time: according to measured curve position and earthquake trajectory position, being provided by reference lamina and continue to bore Into, increasing hole angle or drop angle instruction, until finishing drilling.
Beneficial effects of the present invention:
1) in human intervention, the relationship between geological layering and known seismic event is provided, it is quasi- by setting Conjunction condition, using the best fit position between Computer Automatic Search geological layering and seismic reflector, by seismic interpretation Control errors are within 2m, depth relationship when establishing more accurate.
2) when the basis established of depth relationship on, neighbouring seismic reflection lineups of positioning dessert body, and carry out same-phase Earthquake tracking, according to dessert thickness, earthquake profile position where calculating dessert top bottom determines dessert time window, in drilling well By live fluoroscopic drill bit position in journey, optimization is with brill track.
3) the area stratum standard log curve is established using fixed well and stratigraphic dip, by standard log curve projection to setting Count horizontal well path, generate prediction log, in Horizontal Well Drilling into the process, by comparison with bore curve and prediction curve it Between relationship, adjust track in time, improve reservoir-encountered rate.
Detailed description of the invention:
Detailed description will be given by the following detailed implementation and drawings by the present invention for ease of explanation,.
Fig. 1 is the broken line schematic diagram of seismic event of the present invention and take-off depth;
Window determines schematic diagram when Fig. 2 is earthquake dessert of the invention;
Fig. 3 is seismic interpretation schematic diagram of the invention;
Fig. 4 is calculating GR and resistivity schematic diagram of the invention;
Fig. 5 is earthquake of the invention and track schematic diagram;
Fig. 6 is the structural schematic diagram of log A of the invention;
Fig. 7 is the structural schematic diagram of log B of the invention.
Specific embodiment:
First embodiment: as shown in figs. 1-7, present embodiment uses following technical scheme: a kind of paper shale oil Horizontal well is with brill method for controlling trajectory, comprising the following steps:
1) optimal artificial synthesized seismic channel is determined, window when generating earthquake dessert according to horizontal well dessert section seismic velocity;
2) dessert segment standard log response curve is established using the corresponding form of fixed well log, and according to shale grease Horizontal well planned course generates prediction log;
3) in drilling course, according to Real-time Logging curve, reference mark interval well logging absolute value and amplitude of variation are adjusted in time It is whole with bore track.
Wherein, it in the step 1), determines optimal artificial synthesized seismic channel, is generated according to horizontal well dessert section seismic velocity Window determines 2 parts when window mainly includes holostrome Duan Jingzhen calibration and horizontal well earthquake dessert when earthquake dessert.
Wherein, holostrome Duan Jingzhen calibration the following steps are included:
A., optimal artificial synthesized seismic channel is determined;
Individual well geological layering by well head to shaft bottom number consecutively be G1, G2, Gn, utilize acoustic travel time logging bent Line (AC) artificial synthesized seismic channel correspond to time depth be Gt1, Gt2, Gtn, seismic interpretation interface corresponds to time depth From top to bottom number consecutively be T1, T2, Tn.
It defines well and shakes degree of fitting DF formula:
Wherein, Gt1, Gt2, Gtn be that artificial synthetic seismogram corresponds to time depth, T1, T2, Tn be Seismic interpretation interface corresponds to time depth.Well shake degree of fitting df value is smaller, and it is higher to represent well shake fitting degree.
Mean value is denoted as avgAC in interval transit time curve measuring point depth bounds.
WhenWhen, avgAC is bigger than normal or less than normal, causes artificial synthesized earthquake result and earthquake reflective data It is not inconsistent, needs artificially to adjust avgAC value at this time, artificial synthesized record is corrected, reduces DF value according to this, it is quasi- to improve well shake Conjunction degree is higher.
WhenWhen, avgAC is bigger than normal, and the corresponding df value of avgAC is denoted as DF0, and interval transit time value is integrally subtracted Small fixed value a (generally taking 5~10) re-starts artificial synthesized earthquake using the interval transit time curve newly obtained and notes down, right The df value answered is denoted as df1.The step is repeated, until Successively determine DF2, DF3, DFn, take DF The corresponding artificial synthesized seismic channel of minimum value is optimal synthetic seismogram.
WhenWhen, avgAC is less than normal, and the corresponding DF value of avgAC is denoted as DF0, and interval transit time value is integrally increased Add a fixed value a (generally taking 5~10), re-starts artificial synthesized earthquake using the interval transit time curve newly obtained and note down, it is right The df value answered is denoted as DF1.The step is repeated, until Successively determine DF2, DF3, DFn, take DF The corresponding artificial synthesized seismic channel of minimum value is optimal artificial synthesized seismic channel.
B, the artificial synthesized seismic channel of layer series optimization is divided to determine;
For lineups included in the interface seismic interpretation T1 and T2 from the bottom to top respectively number be t1, t2, Tm, take-off depth number consecutively corresponding to individual well be g1, g2, gm, it is deep that optimal artificial synthesized seismic channel correspond to the time Degree be gt1, gt2, gtm.
It defines well and shakes degree of fitting df formula:
Wherein, gt1, gt2, gtn be that artificial synthetic seismogram corresponds to time depth, t1, t2, tn be Seismic interpretation interface corresponds to time depth.Well shake degree of fitting df value is smaller, and it is higher to represent well shake fitting degree.
Mean value is denoted as avgac in interval transit time curve measuring point depth bounds.
WhenWhen, avgac is bigger than normal or less than normal, causes artificial synthesized earthquake result and earthquake reflective data not Symbol, needs artificially adjust avgac value at this time, is corrected to artificial synthesized record, according to this to reduce df value, improves well and shakes and is fitted Degree is higher.
WhenWhen, avgAC is bigger than normal, and the corresponding df value of avgAC is denoted as df0, and interval transit time value is integrally subtracted Small fixed value a (generally taking 3~6) re-starts artificial synthesized earthquake using the interval transit time curve newly obtained and notes down, corresponding Df value be denoted as df1.The step is repeated, until Successively determine df2, df3, dfn, take df most It is small to be worth corresponding artificial synthesized seismic channel optimal synthetic seismogram between the interface T1 and T2.
WhenWhen, avgAC is less than normal, and the corresponding df value of avgAC is denoted as df0, and interval transit time value is integrally increased Add a fixed value a (generally taking 3~6), re-starts artificial synthesized earthquake using the interval transit time curve newly obtained and note down, it is corresponding Df value be denoted as df1.The step is repeated, until Successively determine df2, df3, dfn, take df most It is small to be worth corresponding artificial synthesized seismic channel optimal artificial synthesized seismic channel between the interface T1 and T2.
Repeat the above steps, successively optimize T2~T3, T3~T4, the axis of reflection between Tn-1~Tn so that most Excellent artificial synthesized seismic channel is in optimum state, and establishes time and depth transfer relationship.
Wherein, when the earthquake dessert window it is determining the following steps are included:
A steady reflection lineups are found in the corresponding seismic reflector of dessert body, by Event tracking, Dessert steady reflection interface t, seismic velocity v in dessert section depth bounds are found in three-dimensional space.On steady reflection interface Dessert body is denoted as Su, and thickness is denoted as du, and dessert body is denoted as Sd under steady reflection interface, and thickness is denoted as dd, and dessert time window is (t-du/v, t+dd/v).If steady reflection lineups are in the top of dessert body, dessert time window is (t, t+dd/v). If steady reflection lineups are in the bottom of dessert body, dessert time window is (t-du/v, t).
Wherein, in the step 2), dessert segment standard log response is established using the corresponding form of fixed well log Curve, and prediction log is generated according to shale oil horizontal well design track.
1. standard log curve calculates;
The geologic body limited by G1 and G2, O point are the intersection points of well 1 Yu geological interface G1, and Q point is well 1 and geological interface G2 Intersection point, geologic dip α, P are any points in drilling well, which is β, and P ' is that P point (hangs down in the projection of vertical direction It is deep), P " is projection of the P point on ground layer vertical.OP is known length section, and OP is ground layer vertical projected length (OP "=OP Sin α/cos β), for any point P, unique subpoint P " can be found in OQ ".OQ " is the true thickness of geologic body, referred to as For standard paragraphs, assigns log value measured by all measuring point P to OP ", both obtained geologic body standard log curve.
The step is repeated, the standard log curve A of all geologic bodies is obtained.
2. generating horizontal wellbore logging curve according to planned course;
The intersection point A of planned well and geological interface G1, the intersection point with geological interface G2 are B, and P is any designed on well location Point, the point design hole angle are β, and P point point corresponding with standard log curve is P ", and geologic body inclination angle is that the α hole angle is β then has AP=OQ " cos β/sin α), for any point P, unique subpoint P " can be found therewith in standard paragraphs OQ " It is corresponding.By establishing the projection relation between OQ " and AB, the corresponding standard log curve values of P " point are assigned to P point Planned well is obtained in the intracorporal prediction log of the geology.
The step is repeated, prediction log is carried out to geologic body in all design well sections and is calculated, and then obtains all set It counts the full well section of well and predicts log B.
Wherein, in the step 3), in drilling course, the relationship between measured curve and prediction curve is compared, and according to Time and depth transfer relationship, in time optimization are as follows with Trace step is bored:
(1) well logging curve monitors in real time: real time contrast surveys log and predicts the relationship between log, According to the distance between measuring instrument and drill bit, drill bit position is grasped in time;
(2) earthquake trajectory real-time monitoring: utilizing time and depth transfer relationship, wellbore trace projected in seismic data body, and According to brill track position relative position in dessert body;
(3) optimize well track in time: according to measured curve position and earthquake trajectory position, being provided by reference lamina and continue to bore Into, increasing hole angle or drop angle instruction, until finishing drilling.
Second embodiment:
1) optimal artificial synthesized seismic channel is determined, window when generating earthquake dessert according to horizontal well dessert section seismic velocity.
1. holostrome Duan Jingzhen is demarcated.
A. optimal artificial synthesized seismic channel is determined.
The well section co-exists in G1、G2And G3Totally 3 geological layerings, artificial synthesizedly using acoustic travel time logging curve (AC) It is G that shake road, which corresponds to time depth,t1、Gt2And Gtn, seismic interpretation interface corresponds to time depth, and from top to bottom number consecutively is T1、T2With Tn, as shown in attached drawing 1-A.
It defines well and shakes degree of fitting DF formula:
Wherein, Gt1、Gt2、···、GtnTime depth, T are corresponded to for artificial synthetic seismogram1、T2、···、TnFor ground Shake explains that interface corresponds to time depth.
By computer search DF minimum value, optimal location corresponding to geological interface is determined, as shown in attached drawing 1-B.
B. optimize T1And T2Artificial synthesized seismic channel between reflecting interface.
Seismic interpretation T1And T2Number is t to included lineups respectively from the bottom to top in interface1、t2、···、t15, single Take-off depth number consecutively corresponding to well is g1、g2、···、g15, optimal artificial synthesized seismic channel corresponds to time depth and is gt1、gt2、···、gt15
It defines well and shakes degree of fitting df formula:
Wherein, gt1、gt2、···、gtnTime depth, t are corresponded to for artificial synthetic seismogram1、t2、···、tnFor ground Shake explains that interface corresponds to time depth.By computer search df minimum value, optimal location corresponding to each lineups is determined, such as Shown in attached drawing 1-C.
2. window determines when earthquake dessert.
(buried depth 3104.46m) and seismic reflection zero phase phase (time depth 25063.62ms) are right at the top of the well dessert It answers, dessert thickness d d is 6.1m, buried depth 3110.56m, time depth 2509.88ms, and dessert body seismic velocity is 1.87m/ms (double-pass reflection time).The dessert time, area's dessert plane distribution was stablized with a thickness of 3.26ms, window when earthquake dessert For (zero phase, zero phase+3.26ms), window is as shown in Fig. 2 when the earthquake dessert of generation.
2) dessert segment standard log response curve is established using the corresponding form of fixed well log, and according to shale grease Horizontal well planned course generates prediction log.
For C1 on C2 dessert section top, stratigraphic dip is about 10 ° of α, is converted by projection, C1 and C2 segment standard log As shown in Fig. 3.According to the track design level Jing Jing, C1 and C2 layers is projected on horizontal well, as a result as shown in Fig. 4.
Wherein, in step 3), in drilling course, compare the relationship between measured curve and prediction curve, and according to when Deep transformational relation, optimization in time is with brill track.
In drilling process, prediction curve and measured curve similarity demonstrate the reliable of prediction log up to 90% Property, and 80 meters of predictions in advance enter and leave window point position, realize that horizontal well smoothly enters window.
Stratigraphic dip change point is found in earthquake dessert cabinet as TRAJECTORY CONTROL OF HORIZONTAL WELL point.Incline on 1 stratum of control point Angle increases to 10.4 ° by 9.3 °, when rotary head is close to control point 1, issues drop angle instruction, hole angle is reduced to by 81.7 ° 89.6°.2 stratigraphic dip of control point increases to 8.2 ° by 10.4 °, when rotary head is close to control point 2, increasing hole angle instruction is issued, by well Oblique angle increases to 91.8 ° by 89.6 °.3 stratigraphic dip of control point increases to 9.5 ° by 8.2 °, when rotary head is close to control point 3, hair Drop angle instructs out, and hole angle is reduced to 90.5 ° by 91.8 °.
SL1 horizontal well is satisfactorily completed by the above method with track optimizing work is bored, horizontal segment length 1019m, oil reservoir, which bores, meets Rate 100%, completes scheduled target, which obtains high yield commercial oil by volume fracturing, obtains better effects.
The above shows and describes the basic principles and main features of the present invention and the advantages of the present invention, the technology of the industry Personnel are it should be appreciated that the present invention is not limited to the above embodiments, and the above embodiments and description only describe this The principle of invention, without departing from the spirit and scope of the present invention, various changes and improvements may be made to the invention, these changes Change and improvement all fall within the protetion scope of the claimed invention, the claimed scope of the invention by appended claims and its Equivalent thereof.

Claims (5)

1. a kind of paper shale grease horizontal well is with brill method for controlling trajectory, which comprises the following steps:
1) optimal artificial synthesized seismic channel, is determined, window when generating earthquake dessert according to horizontal well dessert section seismic velocity;
2) dessert segment standard log response curve, is established using the corresponding form of fixed well log, and according to shale oil level Well planned course generates prediction log;
3), in drilling course, according to Real-time Logging curve, reference mark interval well logging absolute value and amplitude of variation adjust in time with Bore track.
2. a kind of paper shale grease horizontal well according to claim 1 is with brill method for controlling trajectory, it is characterised in that: described In step 1), optimal artificial synthesized seismic channel is determined, window is main when generating earthquake dessert according to horizontal well dessert section seismic velocity 2 parts are determined with window when horizontal well earthquake dessert including holostrome Duan Jingzhen calibration.
3. a kind of paper shale grease horizontal well according to claim 2 is with brill method for controlling trajectory, it is characterised in that: described Holostrome Duan Jingzhen calibration the following steps are included:
A., optimal artificial synthesized seismic channel is determined;
Individual well geological layering by well head to shaft bottom number consecutively be G1, G2, Gn, utilize acoustic travel time logging curve (AC) artificial synthesized seismic channel correspond to time depth be Gt1, Gt2, Gtn, seismic interpretation interface correspond to time depth by Up to lower number consecutively be T1, T2, Tn;
It defines well and shakes degree of fitting DF formula:
Wherein, Gt1, Gt2, Gtn be that artificial synthetic seismogram correspond to time depth, T1, T2, Tn be earthquake Explain that interface corresponds to time depth.Well shake degree of fitting df value is smaller, and it is higher to represent well shake fitting degree;
Mean value is denoted as avgAC in interval transit time curve measuring point depth bounds;
WhenWhen, avgAC is bigger than normal or less than normal, and cause artificial synthesized earthquake result not to be inconsistent with earthquake reflective data, It needs artificially to adjust avgAC value at this time, artificial synthesized record is corrected, reduce DF value according to this, improve well shake fitting journey It spends higher;
WhenWhen, avgAC is bigger than normal, and the corresponding df value of avgAC is denoted as DF0, and interval transit time value is integrally reduced one Fixed value a (generally takes 5~10), re-starts artificial synthesized earthquake using the interval transit time curve newly obtained and notes down, corresponding Df value is denoted as df1.The step is repeated, until Successively determine DF2, DF3, DFn, take DF minimum Being worth corresponding artificial synthesized seismic channel is optimal synthetic seismogram;
WhenWhen, avgAC is less than normal, and the corresponding DF value of avgAC is denoted as DF0, and interval transit time value is integrally increased by one Fixed value a (generally takes 5~10), re-starts artificial synthesized earthquake using the interval transit time curve newly obtained and notes down, corresponding Df value is denoted as DF1;The step is repeated, until Successively determine DF2, DF3, DFn, take DF minimum Being worth corresponding artificial synthesized seismic channel is optimal artificial synthesized seismic channel;
B, the artificial synthesized seismic channel of layer series optimization is divided to determine;
For lineups included in the interface seismic interpretation T1 and T2 from the bottom to top respectively number be t1, t2, tm, Take-off depth number consecutively corresponding to individual well be g1, g2, gm, optimal artificial synthesized seismic channel corresponds to time depth For gt1, gt2, gtm;
It defines well and shakes degree of fitting df formula:
Wherein, gt1, gt2, gtn be that artificial synthetic seismogram correspond to time depth, t1, t2, tn be earthquake Explain that interface corresponds to time depth.Well shake degree of fitting df value is smaller, and it is higher to represent well shake fitting degree;
Mean value is denoted as avgac in interval transit time curve measuring point depth bounds;
WhenWhen, avgac is bigger than normal or less than normal, and cause artificial synthesized earthquake result not to be inconsistent with earthquake reflective data, It needs artificially to adjust avgac value at this time, artificial synthesized record is corrected, reduce df value according to this, improve well shake fitting journey It spends higher;
WhenWhen, avgAC is bigger than normal, and the corresponding df value of avgAC is denoted as df0, and interval transit time value is integrally reduced one Fixed value a (generally takes 3~6), re-starts artificial synthesized earthquake using the interval transit time curve newly obtained and notes down, corresponding df Value is denoted as df1;The step is repeated, until Successively determine df2, df3, dfn, take df minimum value Corresponding artificial synthesized seismic channel optimal synthetic seismogram between the interface T1 and T2;
WhenWhen, avgAC is less than normal, and the corresponding df value of avgAC is denoted as df0, and interval transit time value is integrally increased by one Fixed value a (generally takes 3~6), re-starts artificial synthesized earthquake using the interval transit time curve newly obtained and notes down, corresponding df Value is denoted as df1;The step is repeated, until Successively determine df2, df3, dfn, take df minimum value Corresponding artificial synthesized seismic channel optimal artificial synthesized seismic channel between the interface T1 and T2;
Repeat the above steps, successively optimize T2~T3, T3~T4, the axis of reflection between Tn-1~Tn so that optimal people Work synthetic seismogram is in optimum state, and establishes time and depth transfer relationship.
4. a kind of paper shale grease horizontal well according to claim 2 is with brill method for controlling trajectory, it is characterised in that: described When earthquake dessert window it is determining the following steps are included:
A steady reflection lineups are found in the corresponding seismic reflector of dessert body, by Event tracking, in three-dimensional Dessert steady reflection interface t, seismic velocity v in dessert section depth bounds are found in space.Dessert on steady reflection interface Body is denoted as Su, and thickness is denoted as du, and dessert body is denoted as Sd under steady reflection interface, and thickness is denoted as dd, and dessert time window is (t- Du/v, t+dd/v);If steady reflection lineups are in the top of dessert body, dessert time window is (t, t+dd/v);If steady Determine the bottom that reflection line-ups are in dessert body, then dessert time window is (t-du/v, t);
Wherein, in the step 2), dessert segment standard log response curve is established using the corresponding form of fixed well log, And prediction log is generated according to shale oil horizontal well design track;
1. standard log curve calculates;
The geologic body limited by G1 and G2, O point are the intersection points of well 1 Yu geological interface G1, and Q point is the friendship of well 1 Yu geological interface G2 Point, geologic dip α, P are any points in drilling well, which is β, and P ' is projection (vertical depth) of the P point in vertical direction, P " It is projection of the P point on ground layer vertical;OP is known length section, and OP is ground layer vertical projected length (OP "=OPsin α/cos β), for any point P, unique subpoint P " can be found in OQ ";OQ " is the true thickness of geologic body, referred to as standard Section assigns log value measured by all measuring point P to OP ", had both obtained geologic body standard log curve;
The step is repeated, the standard log curve graph A of all geologic bodies is obtained;
2. generating horizontal wellbore logging curve according to planned course;
The intersection point A of planned well and geological interface G1, the intersection point with geological interface G2 are B, and P is any point designed on well location, should Point design hole angle is β, and P point point corresponding with standard log curve is P ", and it is β that geologic body inclination angle, which is the α hole angle, then has AP=OQ " cos β/sin α), for any point P, unique subpoint P " can be found in contrast in standard paragraphs OQ " It answers;By establishing the projection relation between OQ " and AB, assigns the corresponding standard log curve values of P " point to P point, can be obtained Planned well is in the intracorporal prediction log of the geology;
The step is repeated, prediction log is carried out to geologic body in all design well sections and is calculated, and then obtains all planned wells Full well section predicts log B.
5. a kind of paper shale grease horizontal well according to claim 3 is with brill method for controlling trajectory, it is characterised in that: described In step 3), in drilling course, the relationship between measured curve and prediction curve is compared, and according to time and depth transfer relationship, in time Optimization is as follows with Trace step is bored:
(1) well logging curve monitors in real time: real time contrast surveys log and predicts the relationship between log, according to The distance between measuring instrument and drill bit grasp drill bit position in time;
(2) earthquake trajectory real-time monitoring: utilizing time and depth transfer relationship, wellbore trace projected in seismic data body, and according to With brill track position relative position in dessert body;
(3) in time optimization well track: according to measured curve position and earthquake trajectory position, provided by reference lamina continue creep into, Increasing hole angle or drop angle instruction, until finishing drilling.
CN201910461192.9A 2019-05-30 2019-05-30 Method for controlling while-drilling track of thin-layer shale oil horizontal well Active CN110107282B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910461192.9A CN110107282B (en) 2019-05-30 2019-05-30 Method for controlling while-drilling track of thin-layer shale oil horizontal well

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910461192.9A CN110107282B (en) 2019-05-30 2019-05-30 Method for controlling while-drilling track of thin-layer shale oil horizontal well

Publications (2)

Publication Number Publication Date
CN110107282A true CN110107282A (en) 2019-08-09
CN110107282B CN110107282B (en) 2022-05-20

Family

ID=67492922

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910461192.9A Active CN110107282B (en) 2019-05-30 2019-05-30 Method for controlling while-drilling track of thin-layer shale oil horizontal well

Country Status (1)

Country Link
CN (1) CN110107282B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111308558A (en) * 2020-04-08 2020-06-19 中国石油天然气集团有限公司 Shale gas horizontal well longitudinal wave time difference correction method
CN113126165A (en) * 2020-01-15 2021-07-16 中国石油天然气集团有限公司 Mosaic display method and device for two-dimensional inclined shaft synthetic seismic record
CN113550734A (en) * 2020-04-07 2021-10-26 中国石油天然气集团有限公司 Shale gas horizontal well earthquake tracking method and device while drilling
CN113743175A (en) * 2020-10-21 2021-12-03 中国石油天然气集团有限公司 Stratum inclination angle analysis method and device, computer equipment and storage medium
CN113738345A (en) * 2020-11-25 2021-12-03 中国石油天然气集团有限公司 Method for judging drilling condition based on conventional while-drilling tool
CN113738276A (en) * 2020-05-27 2021-12-03 中国石油化工股份有限公司 Control method and system for drilling horizontal well and drilling system for horizontal well
CN114112787A (en) * 2020-09-01 2022-03-01 中国石油化工股份有限公司 Method for identifying shale gas sweet spot sections of single well

Citations (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5642327A (en) * 1994-06-08 1997-06-24 Exxon Production Research Company Method for creating a gain function for seismic data and method for processing seismic data
CN1710446A (en) * 2005-06-21 2005-12-21 中国石油大学(北京) Method for inversion constituting virtual well data using before-folded seismic wave form
US20060090934A1 (en) * 2004-10-28 2006-05-04 Williams Danny T Formation dip geo-steering method
US20100312534A1 (en) * 2008-02-28 2010-12-09 Shiyu Xu Rock Physics Model For Simulating Seismic Response In Layered Fractured Rocks
CN102140889A (en) * 2011-03-09 2011-08-03 邱世军 Horizontal well logging while drilling geo-steering method
CN102943669A (en) * 2012-10-09 2013-02-27 中国石油化工股份有限公司 High-accuracy depth domain stochastic simulation inversion method suitable for designing track of horizontal well
CN103774989A (en) * 2013-12-10 2014-05-07 刘俊 Horizontal well geosteering analysis method for designing track of well to be drilled through dynamic modeling during drilling
US20140149042A1 (en) * 2012-11-23 2014-05-29 Fugro Geoconsulting, Inc. Method and System for identification of gas hydrates and free gas in geologic beds
CN104975808A (en) * 2015-01-15 2015-10-14 中石化西南石油工程有限公司地质录井分公司 Deep horizontal well drilling track regulating method
CN105093308A (en) * 2014-05-19 2015-11-25 中国石油化工股份有限公司 Well trajectory design method and system
CN105464592A (en) * 2015-03-31 2016-04-06 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 Shale gas horizontal well geosteering method
CN105899749A (en) * 2013-11-08 2016-08-24 普拉德研究及开发股份有限公司 Stratigraphic and structural interpretation of deviated and horizontal wellbores
CN106199727A (en) * 2016-08-30 2016-12-07 中国华能集团清洁能源技术研究院有限公司 A kind of method applying seismic interval velocity identification shale gas dessert
US20160377754A1 (en) * 2015-06-24 2016-12-29 Chevron U.S.A. Inc. Well Placement Using Closure Stress Based Landing Map
US20170167254A1 (en) * 2015-12-11 2017-06-15 Schlumberger Technology Corporation Sedimentary characterization from seismic data
CN106869790A (en) * 2017-02-24 2017-06-20 中石化重庆涪陵页岩气勘探开发有限公司 A kind of quick fine geology guidance method of shale gas horizontal well
CN106894761A (en) * 2017-01-13 2017-06-27 武汉时代地智科技股份有限公司 Using the geosteering method of the geosteering model of time-domain seismic volume
US20170248719A1 (en) * 2014-09-03 2017-08-31 The Board Of Regents For Oklahoma State University Methods of generation of fracture density maps from seismic data
CN108363100A (en) * 2018-01-16 2018-08-03 太原理工大学 Based on the coal bed gas dessert area Seismology and Geology recognition methods for sealing unit and rough set up for safekeeping
CN108442882A (en) * 2018-03-01 2018-08-24 中石化江汉石油工程有限公司 Shale gas extended horizontal well based on 2-d seismic data is with brill geosteering method
CN108625782A (en) * 2018-05-30 2018-10-09 徐芝香 A kind of horizontal well inclination section geosteering method
CN109138837A (en) * 2018-09-12 2019-01-04 中国石油集团西部钻探工程有限公司 Reservoir geology guidance method
CN109190318A (en) * 2018-10-31 2019-01-11 中国石油化工股份有限公司 Inclined shaft well track calculating method based on time measure

Patent Citations (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5642327A (en) * 1994-06-08 1997-06-24 Exxon Production Research Company Method for creating a gain function for seismic data and method for processing seismic data
US20060090934A1 (en) * 2004-10-28 2006-05-04 Williams Danny T Formation dip geo-steering method
CN1710446A (en) * 2005-06-21 2005-12-21 中国石油大学(北京) Method for inversion constituting virtual well data using before-folded seismic wave form
US20100312534A1 (en) * 2008-02-28 2010-12-09 Shiyu Xu Rock Physics Model For Simulating Seismic Response In Layered Fractured Rocks
CN102140889A (en) * 2011-03-09 2011-08-03 邱世军 Horizontal well logging while drilling geo-steering method
CN102943669A (en) * 2012-10-09 2013-02-27 中国石油化工股份有限公司 High-accuracy depth domain stochastic simulation inversion method suitable for designing track of horizontal well
US20140149042A1 (en) * 2012-11-23 2014-05-29 Fugro Geoconsulting, Inc. Method and System for identification of gas hydrates and free gas in geologic beds
CN105899749A (en) * 2013-11-08 2016-08-24 普拉德研究及开发股份有限公司 Stratigraphic and structural interpretation of deviated and horizontal wellbores
CN103774989A (en) * 2013-12-10 2014-05-07 刘俊 Horizontal well geosteering analysis method for designing track of well to be drilled through dynamic modeling during drilling
CN105093308A (en) * 2014-05-19 2015-11-25 中国石油化工股份有限公司 Well trajectory design method and system
US20170248719A1 (en) * 2014-09-03 2017-08-31 The Board Of Regents For Oklahoma State University Methods of generation of fracture density maps from seismic data
CN104975808A (en) * 2015-01-15 2015-10-14 中石化西南石油工程有限公司地质录井分公司 Deep horizontal well drilling track regulating method
CN105464592A (en) * 2015-03-31 2016-04-06 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 Shale gas horizontal well geosteering method
US20160377754A1 (en) * 2015-06-24 2016-12-29 Chevron U.S.A. Inc. Well Placement Using Closure Stress Based Landing Map
US20170167254A1 (en) * 2015-12-11 2017-06-15 Schlumberger Technology Corporation Sedimentary characterization from seismic data
CN106199727A (en) * 2016-08-30 2016-12-07 中国华能集团清洁能源技术研究院有限公司 A kind of method applying seismic interval velocity identification shale gas dessert
CN106894761A (en) * 2017-01-13 2017-06-27 武汉时代地智科技股份有限公司 Using the geosteering method of the geosteering model of time-domain seismic volume
CN106869790A (en) * 2017-02-24 2017-06-20 中石化重庆涪陵页岩气勘探开发有限公司 A kind of quick fine geology guidance method of shale gas horizontal well
CN108363100A (en) * 2018-01-16 2018-08-03 太原理工大学 Based on the coal bed gas dessert area Seismology and Geology recognition methods for sealing unit and rough set up for safekeeping
CN108442882A (en) * 2018-03-01 2018-08-24 中石化江汉石油工程有限公司 Shale gas extended horizontal well based on 2-d seismic data is with brill geosteering method
CN108625782A (en) * 2018-05-30 2018-10-09 徐芝香 A kind of horizontal well inclination section geosteering method
CN109138837A (en) * 2018-09-12 2019-01-04 中国石油集团西部钻探工程有限公司 Reservoir geology guidance method
CN109190318A (en) * 2018-10-31 2019-01-11 中国石油化工股份有限公司 Inclined shaft well track calculating method based on time measure

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
吴俊军等: "吉木萨尔凹陷二叠系芦草沟组上甜点地震预测方法", 《地质评论》 *
吴宗国等: "三维地质导向在地质工程一体化实践中的应用", 《中国石油勘探》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113126165A (en) * 2020-01-15 2021-07-16 中国石油天然气集团有限公司 Mosaic display method and device for two-dimensional inclined shaft synthetic seismic record
CN113550734A (en) * 2020-04-07 2021-10-26 中国石油天然气集团有限公司 Shale gas horizontal well earthquake tracking method and device while drilling
CN111308558A (en) * 2020-04-08 2020-06-19 中国石油天然气集团有限公司 Shale gas horizontal well longitudinal wave time difference correction method
CN113738276A (en) * 2020-05-27 2021-12-03 中国石油化工股份有限公司 Control method and system for drilling horizontal well and drilling system for horizontal well
CN113738276B (en) * 2020-05-27 2024-06-07 中国石油化工股份有限公司 Control method and system for drilling horizontal well and drilling system for horizontal well
CN114112787A (en) * 2020-09-01 2022-03-01 中国石油化工股份有限公司 Method for identifying shale gas sweet spot sections of single well
CN113743175A (en) * 2020-10-21 2021-12-03 中国石油天然气集团有限公司 Stratum inclination angle analysis method and device, computer equipment and storage medium
CN113743175B (en) * 2020-10-21 2024-05-31 中国石油天然气集团有限公司 Stratum inclination analysis method and device, computer equipment and storage medium
CN113738345A (en) * 2020-11-25 2021-12-03 中国石油天然气集团有限公司 Method for judging drilling condition based on conventional while-drilling tool
CN113738345B (en) * 2020-11-25 2024-05-28 中国石油天然气集团有限公司 Method for judging drilling situation based on conventional while-drilling tool

Also Published As

Publication number Publication date
CN110107282B (en) 2022-05-20

Similar Documents

Publication Publication Date Title
CN110107282A (en) A kind of paper shale grease horizontal well is with boring method for controlling trajectory
CN103046868B (en) Integrated geosteering method for horizontal well
CN105464592B (en) A kind of shale gas horizontal well geosteering method
US8875806B2 (en) Formation dip geo-steering method
CN100487488C (en) Zero hypocentral distance vertical seismic section compressional-shear wave data depth field corridor stacked section processing method
CN106894761B (en) Utilize the geosteering method of the geosteering model of time-domain seismic volume
CN106154322B (en) Log curve correction method and apparatus
CN108442882B (en) Shale gas large-displacement horizontal well while-drilling geosteering method based on two-dimensional seismic data
US9182510B2 (en) Methods and systems of incorporating pseudo-surface pick locations in seismic velocity models
CN108316859B (en) Shale gas horizontal well enters method for controlling trajectory before target with geosteering is bored
CN101236257A (en) Oil well location determination technical method
CN109138837A (en) Reservoir geology guidance method
CN106940450B (en) Geosteering method for establishing model based on time-domain seismic volume
CN108625782A (en) A kind of horizontal well inclination section geosteering method
CN106526671B (en) Carbonate reservoir seismic wave response characteristic analysis method
CN105911592B (en) A kind of real three dimensional seismic data long wavelength's static correcting method for boring constraint
US6073079A (en) Method of maintaining a borehole within a multidimensional target zone during drilling
CN111257946B (en) Geophysical drilling guiding method and method for updating stratum seismic velocity
CN107942384A (en) Prediction method for realizing intelligent guiding of directional well by using three-dimensional seismic data
CN109991663A (en) Work area seismic velocity body correction method and device
CN107765309B (en) Depth prediction approach before a kind of brill based on vertical seismic data
CN108732621A (en) Depth prediction technique when a kind of fine with brill based on FFC- resistivity
CN113514876B (en) Design well depth calculation method based on stratum velocity correction
CN112901070B (en) Method for tracking horizontal well while drilling and control method for controlling drilling direction of drill bit
Kolbikov et al. Integrated approach for the horizontal well design and geosteering in thin Jurassic sands

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right

Effective date of registration: 20201202

Address after: 100000 Beijing city Dongcheng District No. 16 Andrew

Applicant after: PetroChina Company Limited

Address before: 300280 Dagang Oilfield Exploration and Development Research Institute, No. 1278 Happy Road, Binhai New Area, Tianjin

Applicant before: DAGANG OIL FIELD OF CNPC

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant