CN108509675A - It is a kind of with boring density log to the computational methods of hole diameter and gap - Google Patents

It is a kind of with boring density log to the computational methods of hole diameter and gap Download PDF

Info

Publication number
CN108509675A
CN108509675A CN201810125774.5A CN201810125774A CN108509675A CN 108509675 A CN108509675 A CN 108509675A CN 201810125774 A CN201810125774 A CN 201810125774A CN 108509675 A CN108509675 A CN 108509675A
Authority
CN
China
Prior art keywords
density
gap
hole diameter
following
mud
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
CN201810125774.5A
Other languages
Chinese (zh)
Other versions
CN108509675B (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.)
China University of Petroleum East China
Original Assignee
China University of Petroleum East China
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 China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN201810125774.5A priority Critical patent/CN108509675B/en
Publication of CN108509675A publication Critical patent/CN108509675A/en
Application granted granted Critical
Publication of CN108509675B publication Critical patent/CN108509675B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00—Survey of boreholes or wells
    • E21B47/12—Means 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
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • E21B49/005—Testing the nature of borehole walls or the formation by using drilling mud or cutting data
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling

Landscapes

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

Abstract

The invention discloses a kind of with density log is bored to the computational methods of hole diameter and gap, belongs to oil reservoir logging technical field, specifically comprises the following steps：X is sought firstLAnd XLFunctional relation；Next seeks gap size tm；Finally seek wellbore radius D.The present invention can cover the measurement blind area using ultrasonic measurement hole diameter in the range of calculated gap, and also wider than ultrasonic measurement gap and hole diameter for the scope of application of mud；In addition, from the point of view of engineering, the present invention solves the problems such as inconvenient maintenance of ultrasonic instrument, maintenance cost are high and instrument expense is very expensive.

Description

It is a kind of with boring density log to the computational methods of hole diameter and gap
Technical field
The invention belongs to oil reservoir logging technical fields, and in particular to a kind of calculating with brill density log to hole diameter and gap Method.
Background technology
In well logging, the variation of hole diameter (borehole size) affects Instrument measuring precision, is not only to eliminate porosity It logs well, the important parameter of litho-density log wellbore effect, the cement amount needed for the well cementation of hole diameter data-evaluation can also be utilized, it is fixed The position of slit band, reflects to visual pattern the contracting occurred in borehole wall stability, such as drilling process during the oil exploitation of position Diameter, landslide bit freezing, wellbore expansion etc..Traditional cable well logging is that caliper logging is carried out after drilling well is complete, in the process that instrument rises Borehole size is measured by the backup arm contact borehole wall.And well logging is that underground letter is measured in the rotation progradation of drilling well Breath, this has prodigious difference with the apparatus structure and form of construction work of wireline logging and in measuring environment, therefore with brill It surveys borehole measurement borehole size and needs contactless measurement method.
Currently, usually using ultrasound wave calliper logging in well logging, using the ultrasonic echo time difference determine instrument with Borehole wall distance simultaneously calculates hole diameter size, but the increase of mud balance can make ultrasonic signal decaying serious, and the reliability of measurement is just It can decline, therefore conventional measurement technology is only applicable in mud balance and is less than 1.3g/cm3Condition, and mud balance be more than 1.3g/ cm3Under conditions of the reliability that measures it is very poor.In addition, ultrasonic wave caliper is popped one's head at work and Instrument shell will produce aftershock And the electric current moment of circuit can make by force apparatus measures generate blind area, the measurement of present ultrasonic wave caliper very much when ultrasonic reflections Blind area is substantially within 3cm.This reduces the precision of the parameter measurements such as density, neutron porosity.Therefore ultrasound is utilized There are significant limitations for wave measurement hole diameter.
Invention content
For the above-mentioned technical problems in the prior art, the present invention propose it is a kind of with bore density log to hole diameter and The computational methods in gap, reasonable design overcome the deficiencies in the prior art, have good effect.
To achieve the goals above, the present invention adopts the following technical scheme that：
It is a kind of with bore density log to the computational methods of hole diameter and gap, specifically comprise the following steps：
Step 1：Seek mud flexible strategy XLAnd XLFunctional relation specifically comprises the following steps：
Step 1.1：Using Monte Carlo simulation to being simulated with brill density log, remote density p is obtainedfar, density of earth formations ρb, mud density ρmAnd gap size tmValue；
Step 1.2：Utilize formula XL=(ρb-ρfar)/(ρb-ρm), calculate XLValue；
Step 1.3：Change gap size tm, density of earth formations ρbWith mud density ρm, many experiments are carried out, scatterplot group is obtained Close (XL1, (ρb-ρm)2, tm1)、(XL2, (ρb-ρm)2, tm2)……(XLi, (ρb-ρm)i, tmi), wherein i is pendulous frequency；
Step 1.4：Obtain XL=f (ρb-ρm, tm) functional relation；
Step 2：Seek gap size tm, specifically comprise the following steps：
Step 2.1：By XL=f (ρb-ρm, tm) functional relation deformed to obtain tm=g (XL, ρb-ρm)；
Step 2.2：With multiple linear regression analysis method, calculated gap size tm, specific formula is as follows：
tm=aXL+bXL 2+c(ρb-ρm)+d(ρb-ρm)2+eXL(ρb-ρm)
Wherein, a, b, c, d, e are coefficient respectively；
Step 3：Wellbore radius D is sought, is specifically comprised the following steps：
Step 3.1：In the size for arbitrarily taking three gaps every 120 ° of position；
Step 3.2：Using following equation find out instrument axle center to the borehole wall distance S1、S2、S3；
S1=r+tm1；
S2=r+tm2；
S3=r+tm3；
Wherein, tm1、tm2、tm3It is gap size, r is the radius of instrument, S1、S2、S3Angle between any two is 120 °,
To angular using d as 3 three of common edge, following 3 trigonometric function equations can be listed：
Wherein, d is the distance in wellbore axle center and instrument axle center, α S1Angle between d；
To equation group find out the radius D of wellbore, it is as a result as follows：
Wherein, T is temporary variable.
General principles are described below：
The basic principle of the present invention is based on a kind of computational methods for boring density log progress hole diameter and gap size.Mud Starch flexible strategy XLIt is by density of earth formations ρbWith the mud density ρ of gap fillingmDifference and gap thickness tmIt is determined, mud flexible strategy XL With the difference (ρ of density of earth formations and mud densityb-ρm) and gap thickness tmBetween functional relation：
XL∝(ρb-ρm)tm
Wherein, ρbIt is density of earth formations value, ρmIt is mud density value, tmIt is the size in gap.
Wherein, to mud flexible strategy XLBe explained as follows：
In with density log is bored, measured using double detector, the density away from the remote detector measurement in source is referred to as remote Density, the density away from the close detector measurement in source are known as nearly density.Due to by gap and filling mud influenced, calculating it is remote Density and nearly density are all not equal to stratum real density, substantially corresponding at this time remote density and nearly density be all density of earth formations and The weighted value of mud density.
ρfar=XLρm+(1-XL)ρb
By deforming mud flexible strategy XLFor：
XL=(ρb-ρfar)/(ρb-ρm)
Wherein, ρfarIt is the density value of far detector measurement stratum.
Advantageous effects caused by the present invention：
Calculating of present invention on the basis of with density log is bored to gap and hole diameter, compared to utilization ultrasonic ranging Method can cover the blind area using ultrasonic measurement in the size range of calculated gap, and in the mud scope of application Also it improves to some extent；In addition, from the point of view of engineering, the present invention solve the inconvenient maintenance of ultrasonic instrument, maintenance cost it is high and The problems such as instrument expense is very expensive.
Description of the drawings
Fig. 1 is gap and wellbore, the geometrical relationship schematic diagram of instrument.
Wherein, the 1- boreholes wall；2- apparatus walls；3- wellbores axle center；4- instruments axle center；5- detectors；6- detects position of stratum；7- The angle of instrument rotation；8- instruments axle center is to borehole wall distance S1With wellbore axle center to the angle of instrument axle center distance；9- instrument shafts Distance S of the heart to the borehole wall1；Distance S of the 10- instruments axle center to the borehole wall2；Distance S of the 11- instruments axle center to the borehole wall3。
Specific implementation mode
Below in conjunction with the accompanying drawings and specific implementation mode invention is further described in detail：
Gap and wellbore, the geometrical relationship of instrument are as shown in Figure 1.
It is a kind of with bore density log to the computational methods of hole diameter and gap, specifically comprise the following steps：
Step 1：Seek mud flexible strategy XLAnd XLFunctional relation specifically comprises the following steps：
Step 1.1：Using Monte Carlo simulation to being simulated with brill density log, remote density p is obtainedfar, density of earth formations ρb, mud density ρmAnd gap size tmValue；
Step 1.2：Utilize formula XL=(ρb-ρfar)/(ρb-ρm), calculate XLValue；
Step 1.3：Change gap size tm, density of earth formations ρbWith mud density ρm, many experiments are carried out, scatterplot group is obtained Close (XL1, (ρb-ρm)2, tm1)、(XL2, (ρb-ρm)2, tm2)……(XLi, (ρb-ρm)i, tmi), wherein i is pendulous frequency；
Step 1.4：Obtain XL=f (ρb-ρm, tm) functional relation；
Step 2：Seek gap size tm, specifically comprise the following steps：
Step 2.1：By XL=f (ρb-ρm, tm) functional relation deformed to obtain tm=g (XL, ρb-ρm)；
Step 2.2：With multiple linear regression analysis method, calculated gap size tm, specific formula is as follows：
tm=aXL+bXL 2+c(ρb-ρm)+d(ρb-ρm)2+eXL(ρb-ρm)
Wherein, a, b, c, d, e are coefficient respectively；
Step 3：Wellbore radius D is sought, is specifically comprised the following steps：
Step 3.1：In the size for arbitrarily taking three gaps every 120 ° of position；
Step 3.2：Using following equation find out instrument axle center to the borehole wall distance S1、S2、S3；
S1=r+tm1；
S2=r+tm2；
S3=r+tm3；
Wherein, tm1、tm2、tm3It is gap size, r is the radius of instrument, S1、S2、S3Angle between any two is 120 °,
To angular using d as 3 three of common edge, following 3 trigonometric function equations can be listed：
Wherein, d is the distance in wellbore axle center and instrument axle center, α S1Angle between d；
To equation group find out the radius D of wellbore, it is as a result as follows：
Wherein, T is temporary variable.
The platform of this experiment is Monte Carlo simulation software, simulates data, is obtained with the computational methods processing of the present invention Gap and hole diameter (such as table 1, table 2) is compared with the gap of design and hole diameter, by analytical error come determine the present invention Feasibility, table 1 are the comparisons of design gaps and calculated gap, and table 2 is the wellbore radius of design and calculates the comparison of radius.
The comparison of table 1 design gaps and calculated gap
2 design radial of table and the comparison for calculating radius
Serial number Design radial cm Rotation angle/degree Gap size/cm Calculate radius/cm Error
1 12.9 0 0.093188 12.99319 0.09
2 12.9 22.5 0.310813 13.05561 0.16
3 12.9 45 0.542498 12.83261 0.07
4 12.9 67.5 1.358426 12.92321 0.02
5 12.9 90 1.9991 12.64862 0.25
6 12.9 112.5 3.032269 12.68757 0.21
7 12.9 135 3.976884 12.73243 0.17
8 12.9 157.5 4.691901 12.81796 0.08
9 12.9 180 5.011733 12.91173 0.01
From table 1 above it can be seen that the gap width for the gap width and design being calculated using the present invention is made comparisons, the two Difference is up to 0.323794 (i.e. error), mean error 0.1641107cm, between being determined with the utilization ultrasonic echo time difference It is very nearly the same both in precision for gap, but the range of calculated gap of the present invention covers the measurement blind area of ultrasonic caliper, And for the limitation also opposite reduction that mud property requires；The error of the hole diameter calculated as can be seen from Table 2 by the present invention Very little, substantially close to the borehole size of design.By Tables 1 and 2 it is found that the method for the present invention can effectively calculate gap thickness With hole diameter size, it can be seen that, method provided by the invention disclosure satisfy that the requirement that result of calculation is obtained in engineering.
Certainly, above description is not limitation of the present invention, and the present invention is also not limited to the example above, this technology neck The variations, modifications, additions or substitutions that the technical staff in domain is made in the essential scope of the present invention should also belong to the present invention's Protection domain.

Claims (1)

1. a kind of with boring density log to the computational methods of hole diameter and gap, it is characterised in that：Specifically comprise the following steps：
Step 1：Seek mud flexible strategy XLAnd XLFunctional relation specifically comprises the following steps：
Step 1.1：Using Monte Carlo simulation to being simulated with brill density log, remote density p is obtainedfar, density of earth formations ρb, mud Pulp density ρmAnd gap size tmValue；
Step 1.2：Utilize formula XL=(ρb-ρfar)/(ρb-ρm), calculate XLValue；
Step 1.3：Change gap size tm, density of earth formations ρbWith mud density ρm, many experiments are carried out, scatterplot combination is obtained (XL1, (ρb-ρm)2, tm1)、(XL2, (ρb-ρm)2, tm2)……(XLi, (ρb-ρm)i, tmi), wherein i is pendulous frequency；
Step 1.4：Obtain XL=f (ρb-ρm, tm) functional relation；
Step 2：Seek gap size tm, specifically comprise the following steps：
Step 2.1：By XL=f (ρb-ρm, tm) functional relation deformed to obtain tm=g (XL, ρb-ρm)；
Step 2.2：With multiple linear regression analysis method, calculated gap size tm, specific formula is as follows：
tm=aXL+bXL 2+c(ρb-ρm)+d(ρb-ρm)2+eXL(ρb-ρm)
Wherein, a, b, c, d, e are coefficient respectively；
Step 3：Wellbore radius D is sought, is specifically comprised the following steps：
Step 3.1：In the size for arbitrarily taking three gaps every 120 ° of position；
Step 3.2：Using following equation find out instrument axle center to the borehole wall distance S1、S2、S3；
S1=r+tm1；
S2=r+tm2；
S3=r+tm3；
Wherein, tm1、tm2、tm3It is gap size, r is the radius of instrument, S1、S2、S3Angle between any two is 120 °, to d 3 three for common edge are angular, can list following 3 trigonometric function equations：
Wherein, d is the distance in wellbore axle center and instrument axle center, α S1Angle between d；
To equation group find out the radius D of wellbore, it is as a result as follows：
Wherein, T is temporary variable.
CN201810125774.5A 2018-02-08 2018-02-08 Method for calculating hole diameter and clearance by density logging while drilling Active CN108509675B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810125774.5A CN108509675B (en) 2018-02-08 2018-02-08 Method for calculating hole diameter and clearance by density logging while drilling

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810125774.5A CN108509675B (en) 2018-02-08 2018-02-08 Method for calculating hole diameter and clearance by density logging while drilling

Publications (2)

Publication Number Publication Date
CN108509675A true CN108509675A (en) 2018-09-07
CN108509675B CN108509675B (en) 2021-10-26

Family

ID=63375594

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810125774.5A Active CN108509675B (en) 2018-02-08 2018-02-08 Method for calculating hole diameter and clearance by density logging while drilling

Country Status (1)

Country Link
CN (1) CN108509675B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117148331A (en) * 2023-10-31 2023-12-01 中海油田服务股份有限公司 Test system and test method of while-drilling ultrasonic imaging instrument

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5513528A (en) * 1994-01-14 1996-05-07 Schlumberger Technology Corporation Logging while drilling method and apparatus for measuring standoff as a function of angular position within a borehole
US6044326A (en) * 1999-01-15 2000-03-28 Dresser Industries, Inc. Measuring borehole size
US20070040110A1 (en) * 2005-08-22 2007-02-22 Schlumberger Technology Corporation Measuring wellbore diameter with an LWD instrument using compton and photoelectric effects
US20120166087A1 (en) * 2009-04-17 2012-06-28 Michael Evans Method of Determining Density of Underground Formations Using Neutron-Gamma Ray Measurements
CN103233727A (en) * 2013-05-13 2013-08-07 中国石油大学(华东) Inversion method of stratum shear wave velocity radial sections
CN106646642A (en) * 2016-12-29 2017-05-10 中国石油大学(华东) Scanning type radioactive hole diameter measurement device and method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5513528A (en) * 1994-01-14 1996-05-07 Schlumberger Technology Corporation Logging while drilling method and apparatus for measuring standoff as a function of angular position within a borehole
US6044326A (en) * 1999-01-15 2000-03-28 Dresser Industries, Inc. Measuring borehole size
US20070040110A1 (en) * 2005-08-22 2007-02-22 Schlumberger Technology Corporation Measuring wellbore diameter with an LWD instrument using compton and photoelectric effects
US20120166087A1 (en) * 2009-04-17 2012-06-28 Michael Evans Method of Determining Density of Underground Formations Using Neutron-Gamma Ray Measurements
CN103233727A (en) * 2013-05-13 2013-08-07 中国石油大学(华东) Inversion method of stratum shear wave velocity radial sections
CN106646642A (en) * 2016-12-29 2017-05-10 中国石油大学(华东) Scanning type radioactive hole diameter measurement device and method

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
M MICKAEL 等: "Standoff compensation and hole size correction of a new LWD density/neutron logging system", 《SPE TECHNICAL CONFERENCE AND EXHIBITION》 *
P.COOPER 等: "Borehole and Invasion Effects of Formate-Based Mud Systems on LWD Density,Neutron,PE,and Gamma-rays Logs", 《SPE ANNUAL TECHNICAL CONFERENCE AND EXHIBITION》 *
于华伟 等: "蒙特卡罗模拟研究水平井密度测井中岩屑层的影响", 《核技术》 *
倪卫宁 等: "基于超声测距的随钻井径测量系统设计", 《科技导报》 *
朱頔 等: "随钻密度测井间隙对探测器晶体影响的蒙特卡罗模拟", 《测井技术》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117148331A (en) * 2023-10-31 2023-12-01 中海油田服务股份有限公司 Test system and test method of while-drilling ultrasonic imaging instrument
CN117148331B (en) * 2023-10-31 2023-12-26 中海油田服务股份有限公司 Test system and test method of while-drilling ultrasonic imaging instrument

Also Published As

Publication number Publication date
CN108509675B (en) 2021-10-26

Similar Documents

Publication Publication Date Title
CN103233727B (en) Inversion method of stratum shear wave velocity radial sections
US9354050B2 (en) Borehole characterization
CN105550448B (en) Pre- drilling well three-dimensional well eye modeling method based on wellbore trace design parameter and device
CN103775074B (en) Method for correcting expanding influence of coalbed methane reservoir density logging
RU2599648C1 (en) Detection of location of boundaries of formation on basis of measurements at several depths of tool in well bore
RU2663653C1 (en) Improved estimation of well bore logging based on results of measurements of tool bending moment
US10954780B2 (en) Eccentricity correction algorithm for borehole shape and tool location computations from caliper data
CA2964228C (en) Methods and systems for modeling an advanced 3-dimensional bottomhole assembly
RU2528279C1 (en) Acoustic waveform stacking by grouping based on azimuth angles and/or deviations of logging probe
MX2012004797A (en) System and method for determining stretch or compression of a drill string.
WO2012166111A1 (en) Azimuthal brittleness logging systems and methods
CN105114063A (en) Inclined-well dual laterolog curve correction method
WO2017003517A1 (en) Correcting the effects of deviation and dispersion on sonic log measurements of deviated wells in laminated formations
CN108509675A (en) It is a kind of with boring density log to the computational methods of hole diameter and gap
RU2628660C2 (en) Range measurement using modulated signals
CN207296995U (en) A kind of high accuracy Integral wireless measurement-while-drilling system
Heisig et al. Continuous borehole curvature estimates while drilling based on downhole bending moment measurements
Kang et al. Review on the influence of complex stratum on the drilling trajectory of the drilling robot
NO20171979A1 (en) Passive ranging using acoustic energy originating from a target wellbore
CN106761666A (en) A kind of four probes scattering gamma logging and the method and device of nonlinear data inverting
CN107103165B (en) Method for calculating density of steep dip shear crack line of oil field reservoir
Øy et al. 4D Radar Imaging of Wellbore Geometry While Drilling
CN112906465A (en) Coal measure stratum acoustic curve reconstruction method and system based on stratum factors
Rassadkin et al. Precise Interlateral Spacing for Optimal Stimulation and Enhanced Production in North American Shale
Elahifar et al. Determining the Accuracy of Ultrasonic Caliper Measurements for Real Time Drilling Application

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant