CN109209346A - A method of reducing irregular oscillation influences continuous gyroscopic survey instrument precision - Google Patents

A method of reducing irregular oscillation influences continuous gyroscopic survey instrument precision Download PDF

Info

Publication number
CN109209346A
CN109209346A CN201811326118.8A CN201811326118A CN109209346A CN 109209346 A CN109209346 A CN 109209346A CN 201811326118 A CN201811326118 A CN 201811326118A CN 109209346 A CN109209346 A CN 109209346A
Authority
CN
China
Prior art keywords
survey instrument
gyroscopic survey
gyro
continuous gyroscopic
calibration
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.)
Pending
Application number
CN201811326118.8A
Other languages
Chinese (zh)
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.)
Guizhou Aerospace Kaishan Petroleum Instrument Co Ltd
Original Assignee
Guizhou Aerospace Kaishan Petroleum Instrument Co Ltd
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 Guizhou Aerospace Kaishan Petroleum Instrument Co Ltd filed Critical Guizhou Aerospace Kaishan Petroleum Instrument Co Ltd
Priority to CN201811326118.8A priority Critical patent/CN109209346A/en
Publication of CN109209346A publication Critical patent/CN109209346A/en
Pending legal-status Critical Current

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/02Determining slope or direction

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Geophysics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Gyroscopes (AREA)

Abstract

The invention discloses a kind of methods that reduction irregular oscillation influences continuous gyroscopic survey instrument precision, this method utilizes 3-axis acceleration information and the angular velocity information of combination optical fibre gyro Z axis, the interference outside instrument attitude measurement is rejected, to improve the measurement accuracy of continuous gyroscopic survey instrument comprehensively.This method is able to solve continuous gyroscopic survey instrument and the problem of irregular collision causes accuracy of instrument to decline occurs with tubing and casing during upper go into the well.

Description

A method of reducing irregular oscillation influences continuous gyroscopic survey instrument precision
Technical field
The present invention relates to the applications for applying gyrolevel to carry out continuous dynamic test trails in a kind of oil well logging industry.
Background technique
Caused due to tube or cased well breakage, deformation, wax deposition etc. gyrolevel in upper motion process of going into the well with pipe Irregular collision occurs for wall, transports in the angular movement of inertial space and linearly to the sensitive instrumentation of optical fibre gyro and accelerometer composition Dynamic to interfere, the well track for clearing of navigating in continuous motion process can introduce error, so as to cause accuracy of instrument decline.Often The gyrolevel seen only simply carries out low-pass filtering to 3-axis acceleration information and smoothly obtains new acceleration information, filters The band logical frequency of wave device is fixed, attenuation coefficient and filter order fixation can not change according to hole condition, and method is simple and crude.It is practical The speed gone into the well in the instrument course of work is impossible really at the uniform velocity constant, while random variation can also occur with oil jacket tube wall Impact force, the interference changed at random will necessarily be brought, to acceleration information so as to cause the azimuth and tool-face cleared out Angle error is very big, and precision is also more far short of what is expected than fixed point gyrolevel.
Summary of the invention
The object of the invention is that solve, instrument cannot achieve uniform motion in oil well and tube wall incrustation deformation is irregular Caused by random posture vibratory impulse influence, the side that a kind of reduction irregular oscillation influences continuous gyroscopic survey instrument precision is provided Method.
Method of the invention is that this method is using 3-axis acceleration information and the angular speed of optical fibre gyro Z axis to be combined to believe Breath rejects the interference outside instrument attitude measurement.
Moreover, this method is to constitute inclinometer using double-shaft optical fiber gyroscope and three axis quartz accelerometers, it is used using strapdown Guiding method realizes the real-time measurement of instrument posture and speed, to obtain hole angle, azimuth and tool face azimuth;And utilize oil Well track can not mutate within certain distance, establish data model and reject random disturbances, to effectively improve instrument Measuring accuracy.
Wherein, to accelerometer signal and gyro signal not instead of not simply low-pass filtering treatment, according to external environment Sef-adapting filter is established in variation.
The present invention specifically includes calibration modeling and modification method;The calibration modeling includes location position, angular speed calibration It is demarcated with random vibration.
Wherein, location position above-mentioned is to determine that the calibration factor of table, installation is added to miss using the method for static 8 location positions Difference and gyro and the installation error for adding table, complete the calibration to whole system.
Angular speed calibration above-mentioned is according to the variation range of continuous gyroscopic survey instrument angular speed during upper go into the well point Not She Ding high-precision electric turntable speed, then it is two axis of gyro are horizontal fixed respectively, record the output of two axis gyros; Error coefficient is found out according to gyro error output equation, to show that the constant value drift of gyro, scale factor and gyro installation miss Difference.
Random vibration calibration above-mentioned is to be fixed on continuous gyroscopic survey instrument on three axis electric vibration tables with fixture, simulation Instrument is arranged to random vibration mode in the situation of movement of tube or cased well, according to minimal sampling time interval continuous sampling data; The vibration of electric vibration table interval is spaced static alternate cycles;Gyroscope and accelerometer are analyzed in difference by the data of acquisition Input noise characteristic under vibration condition, model determines the band connection frequency coefficient range of filter according to the collected data With the scale factor of vibration frequency.
Modification method above-mentioned is using denoising is forced in the selection of filtering algorithm, and continuous gyroscopic survey instrument is selecting By time domain evaluation error in judgement characteristic in rate of going into the well on fixed variation range;Continuous gyroscopic survey instrument occurs in tube or cased well Attitude range is about scheduled on signal specific amplitude range within the unit time, is adjusted and is carried out in conjunction with random vibration peg model Fuzzy comparability identification, to have the function that the extraneous vibration interference for rejecting accelerometer and gyroscope.
The innovation of the invention consists in that:
(1), to accelerometer signal and gyro signal not instead of not simply low-pass filtering treatment, become according to external environment Sef-adapting filter is established in change.
(2), the selection of angular speed calibration turntable rate and the foundation of gyroscope error model.
(3), the establishment of random vibration scaling method and the foundation of error model.
(4), modification method.
Compared with prior art, method of the invention using 3-axis acceleration information and combines the angle of optical fibre gyro Z axis fast Information is spent, the interference outside instrument attitude measurement is rejected, to improve the measurement accuracy of continuous gyroscopic survey instrument comprehensively.This method energy It enough solves continuous gyroscopic survey instrument and causes asking for accuracy of instrument decline with the irregular collision of tubing and casing generation during upper go into the well Topic.
The present invention is the mark fixed sum data processing method summed up for many years to the practical experience of oil field necessary instrument.Calibration ginseng Several selections is also one of its core, while being designed using adaptive Real Time Digital Filter, bandpass center frequency f0, excessive band Wide coefficient ct, band logical side frequency fp=(1+n) * f0, band logical cutoff frequency wp=fp* 2 and stopband side frequency fs=fp*(1+ct) etc. bases The variation of external condition adjusts, and by carrying out analysis de-noising to signal, minimizes the distortion factor of signal.
Specific embodiment
Below with reference to embodiment, the invention will be further described, but not as any limitation of the invention.
Embodiment 1:
1. calibration modeling
1) location position
The installation for adding calibration factor, installation error and the gyro of table and adding table is determined using the method for static 8 location positions Error completes the calibration to whole system.The specific layout in 8 positions are as follows:
X-Y-Z Terrestrial gravitation
Northeast day 0 0 -g
Southwestern day 0 0 -g
East northeast 0 0 g
Southern west ground 0 0 g
Its southwest -g 0 0
Its northeast -g 0 0
Ground Nan Xi g 0 0
Ground east northeast g 0 0
When location position, in each position, static testing acceleration meter exports average value fi, obtain the recurrence meter of accelerometer Calculate model:
By the output model of accelerometer:
F=a0+Kaab
It is denoted as:
F=MA
In formula,
F=[f1 f2 … f8]TFor the output of a certain axis accelerometer,For terrestrial gravitation plus Speed input.Each parameter computation model of accelerometer then can be obtained according to Least square-fit:
M=FAT·(AAT)-1
2) angular speed is demarcated
According to continuous gyroscopic survey instrument, high precision electro turn is set separately in the variation range of angular speed during upper go into the well The speed of platform be ± 0.05 °/s, ± 0.1 °/s, ± 0.2 °/s, ± 0.5 °/s, ± 1 °/s, ± 2 °/s, ± 5 °/s, ± 10 °/s, ± 20 °/s, ± 30 °/s, then two axis of gyro are horizontal fixed respectively, the output of two axis gyros is recorded, the sampling time presses It is determined according to 360 ° divided by turntable rate, chronomere s, the sampling time corresponds to range (7200~12) s.It is defeated according to gyro error Equation finds out error coefficient out, to obtain the constant value drift of gyro, scale factor and gyro misalignment.Gyroscope exports mould Type are as follows:
In formula:
Dx、DyAnd DzIt is constant value drift, N for gyrox、NyAnd NzFor gyroscope output, Kx、KyAnd KzFor gyroscope scale Factor, kijFor the installation error of gyro i axis and j axis, ωiFor the input angular velocity of gyro i axis, in the present invention, due to inclinometer Only X and Z axis gyro, so the amount in relation to Y-axis gyro is set as 0.Finally, through-rate test data, which is fitted, to be obtained Obtain each coefficient entry in relation to gyro.
3) random vibration is demarcated
Continuous gyroscopic survey instrument is fixed on three axis electric vibration tables with fixture, movement of the analog meter in tube or cased well Situation, the swing displacement along instrument X, Y-axis is 0~(casing diameter-instrument diameters)/2mm, longitudinal to consider trembling for logging cable Swinging displacement is 0~20mm, and acceleration is 0~3g, and frequency is that 5Hz~100Hz is arranged to random vibration mode, is adopted according to minimum Sample time 0.5ms interval continuous sampling data.Electric vibration table operating mode is set as the vibration of 5s multiple interval, 5s multiple interval Static alternate cycles, stream time 30min.Gyroscope and accelerometer are analyzed in different vibrator bars by the data of acquisition Input noise characteristic under part, model determines band connection frequency coefficient range and the vibration of filter according to the collected data The scale factor of frequency.
2. modification method
Using forcing denoising in the selection of filtering algorithm, continuous gyroscopic survey instrument rate of going into the well on selected becomes Change and presses time domain evaluation error in judgement characteristic in range.The attitude range that continuous gyroscopic survey instrument occurs in tube or cased well is in list It is about scheduled on signal specific amplitude range in the time of position, is adjusted in conjunction with random vibration peg model and carries out fuzzy comparability identification, To have the function that reject the extraneous vibration interference of accelerometer and gyroscope.
Signal measurement model are as follows:
C is the measurement factor of signal, and v (k) is the external interference noise on signal superposition, and test discovery, inclinometer is even Under continuous working condition, the mean value of external interference noise is zero, variance σ2, then recursive estimation filter can be designed are as follows:
Formula first item is the estimated data that sensor weights under current pose, and Section 2 is that sensor is currently adopted in real time Sample data, the method for solving of a (k) and b (k) in formula are as follows:
A (k)=a [1-cb (k)]
Wherein:
According to the attitude data of a sampling instant on instrument, the posture of subsequent time can be estimated, passes through the posture estimated It can be concluded that gyroscope and accelerometer estimate output data in the case where this estimates posture, estimated data and test data carried out Recursive filtering obtains the credible measured value of sensor, and bringing the confidence values into attitude algorithm algorithm can be obtained instrument posture.
It is above specific application example of the invention, there are other embodiments of the invention, all to use equivalent replacement Or the technical solution that equivalent transformation is formed, it all falls within protection scope of the presently claimed invention.

Claims (8)

1. a kind of method for reducing irregular oscillation influence continuous gyroscopic survey instrument precision, it is characterised in that: this method is to utilize 3-axis acceleration information and the angular velocity information for combining optical fibre gyro Z axis reject the interference outside instrument attitude measurement.
2. the method according to claim 1 for reducing irregular oscillation influence continuous gyroscopic survey instrument precision, feature exist In: this method is to constitute inclinometer using double-shaft optical fiber gyroscope and three axis quartz accelerometers, is realized using inertial navigation method The real-time measurement of instrument posture and speed, to obtain hole angle, azimuth and tool face azimuth;And using oil well track one It can not mutate within set a distance, establish data model and reject random disturbances, to effectively improve instrument test precision.
3. the method according to claim 1 for reducing irregular oscillation influence continuous gyroscopic survey instrument precision, feature exist In: to accelerometer signal and gyro signal not instead of not simply low-pass filtering treatment, according to external environment variation establish from Adaptive filter.
4. the method according to claim 1 for reducing irregular oscillation influence continuous gyroscopic survey instrument precision, feature exist In: including calibration modeling and modification method;The calibration modeling includes location position, angular speed calibration and random vibration calibration.
5. the method according to claim 4 for reducing irregular oscillation influence continuous gyroscopic survey instrument precision, feature exist In: the location position is to utilize determining calibration factor, installation error and the gyro for adding table of method of static 8 location positions and add The installation error of table completes the calibration to whole system.
6. the method according to claim 4 for reducing irregular oscillation influence continuous gyroscopic survey instrument precision, feature exist In: the angular speed calibration is that height is set separately in the variation range of angular speed during upper go into the well according to continuous gyroscopic survey instrument The speed of precision electric turntable, then two axis of gyro are horizontal fixed respectively, record the output of two axis gyros;According to gyro Error output equation finds out error coefficient, to obtain the constant value drift of gyro, scale factor and gyro misalignment.
7. the method according to claim 4 for reducing irregular oscillation influence continuous gyroscopic survey instrument precision, feature exist In: the random vibration calibration is to be fixed on continuous gyroscopic survey instrument on three axis electric vibration tables with fixture, and analog meter exists The situation of movement of tube or cased well is arranged to random vibration mode, according to minimal sampling time interval continuous sampling data;Electronic vibration Dynamic platform interval vibration is spaced static alternate cycles;Gyroscope and accelerometer are analyzed in different vibrator bars by the data of acquisition Input noise characteristic under part, model determines band connection frequency coefficient range and the vibration of filter according to the collected data The scale factor of frequency.
8. the method according to claim 4 for reducing irregular oscillation influence continuous gyroscopic survey instrument precision, feature exist In: the modification method is using denoising is forced in the selection of filtering algorithm, and continuous gyroscopic survey instrument is upper what is selected By time domain evaluation error in judgement characteristic in rate of going into the well variation range;Continuous gyroscopic survey instrument becomes in the posture that tube or cased well occurs Change range and be about scheduled on signal specific amplitude range within the unit time, carries out fuzzy phase in conjunction with the adjustment of random vibration peg model It is identified like property, to have the function that the extraneous vibration interference for rejecting accelerometer and gyroscope.
CN201811326118.8A 2018-11-08 2018-11-08 A method of reducing irregular oscillation influences continuous gyroscopic survey instrument precision Pending CN109209346A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811326118.8A CN109209346A (en) 2018-11-08 2018-11-08 A method of reducing irregular oscillation influences continuous gyroscopic survey instrument precision

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811326118.8A CN109209346A (en) 2018-11-08 2018-11-08 A method of reducing irregular oscillation influences continuous gyroscopic survey instrument precision

Publications (1)

Publication Number Publication Date
CN109209346A true CN109209346A (en) 2019-01-15

Family

ID=64995171

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811326118.8A Pending CN109209346A (en) 2018-11-08 2018-11-08 A method of reducing irregular oscillation influences continuous gyroscopic survey instrument precision

Country Status (1)

Country Link
CN (1) CN109209346A (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU94045552A (en) * 1994-12-29 1996-10-10 Центральный научно-исследовательский институт "Дельфин" Gimbal-free gyroscopic inclinometer and method of determination of inclination
JP2004347587A (en) * 2003-04-28 2004-12-09 National Institute Of Advanced Industrial & Technology Dynamic matrix sensitivity measurement device for inertia sensor, and measurement method therefor
CN101566483A (en) * 2009-05-22 2009-10-28 哈尔滨工程大学 Vibration error compensation method of fiber optic gyro strap-down inertia measurement system
CN202255374U (en) * 2011-09-29 2012-05-30 西安格威石油仪器有限公司 System for measuring inclined angle and relative azimuth of oil logging instrument
CN102692239A (en) * 2012-06-14 2012-09-26 辽宁工程技术大学 Fiber optic gyroscope eight-position calibration method based on rotating mechanism
US20120259579A1 (en) * 2011-04-07 2012-10-11 Icefield Tools Corporation Method and apparatus for determining orientation using a plurality of angular rate sensors and accelerometers
CN104406608A (en) * 2014-11-28 2015-03-11 陕西宝成航空仪表有限责任公司 Vibration table for static balance of gyroscope
CN105134171A (en) * 2015-07-29 2015-12-09 北京航空航天大学 Realization method for dual-spindle optic-fiber gyro continuous inclination measuring system
CN105157723A (en) * 2015-07-14 2015-12-16 北京理工大学 Fiber gyro-based strapdown inertial navigation system calibration method

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU94045552A (en) * 1994-12-29 1996-10-10 Центральный научно-исследовательский институт "Дельфин" Gimbal-free gyroscopic inclinometer and method of determination of inclination
JP2004347587A (en) * 2003-04-28 2004-12-09 National Institute Of Advanced Industrial & Technology Dynamic matrix sensitivity measurement device for inertia sensor, and measurement method therefor
US20070073502A1 (en) * 2003-04-28 2007-03-29 National Inst. Of Adv. Industrial Science & Tech. Dynamic matrix sensitivity measuring instrument for inertial sensors, and measuring method therefor
CN101566483A (en) * 2009-05-22 2009-10-28 哈尔滨工程大学 Vibration error compensation method of fiber optic gyro strap-down inertia measurement system
US20120259579A1 (en) * 2011-04-07 2012-10-11 Icefield Tools Corporation Method and apparatus for determining orientation using a plurality of angular rate sensors and accelerometers
CN202255374U (en) * 2011-09-29 2012-05-30 西安格威石油仪器有限公司 System for measuring inclined angle and relative azimuth of oil logging instrument
CN102692239A (en) * 2012-06-14 2012-09-26 辽宁工程技术大学 Fiber optic gyroscope eight-position calibration method based on rotating mechanism
CN104406608A (en) * 2014-11-28 2015-03-11 陕西宝成航空仪表有限责任公司 Vibration table for static balance of gyroscope
CN105157723A (en) * 2015-07-14 2015-12-16 北京理工大学 Fiber gyro-based strapdown inertial navigation system calibration method
CN105134171A (en) * 2015-07-29 2015-12-09 北京航空航天大学 Realization method for dual-spindle optic-fiber gyro continuous inclination measuring system

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
刘宇等: "《固态振动陀螺与导航技术》", 30 September 2010 *
曹玉坤等: "模糊识别和随机森林算法在柴油机振动信号状态识别中的应用", 《装甲兵工程学院学报》 *
李家垒等: "光纤陀螺消噪方法研究", 《电光与控制》 *
李晨曦等: "两轴光纤陀螺连续测斜仪姿态解算方法研究", 《电子测量技术》 *
柴嘉薪等: "一种光纤陀螺随机振动误差高精度建模方法", 《航空兵器》 *
王保华: "《生物医学电子学》", 31 October 1988 *
苏游: "基于FOG捷联系统的误差分析与补偿技术研究", 《CNKI硕士学位论文数据库信息科技》 *

Similar Documents

Publication Publication Date Title
KR101988786B1 (en) Initial alignment of inertial navigation devices
CN109188026B (en) Automatic calibration deep learning method suitable for MEMS accelerometer
JP5068531B2 (en) Method and system for improving the accuracy of inertial navigation measurements using measured and stored gravity gradients
CN105698792A (en) Dynamic MEMS (micro-electromechanical systems) inertial attitude measuring system based on self-adaptive robust integration algorithm
Arias-Lara et al. Assessment of methodologies to estimate displacements from measured acceleration records
Roberts et al. High frequency deflection monitoring of bridges by GPS
CN104453857B (en) A kind of small hole deviation go into the well tiltedly and tool face azimuth dynamic measurement method and device
CN108197339B (en) Dynamic displacement monitoring method of real-time acceleration integral algorithm based on recursive least square method
CN109766812B (en) Post-event compensation method for motion error of gravity gradiometer of rotating accelerometer
CN103104251A (en) Method used for improving accuracy of azimuthal angle and tool face angle of small well inclined downward fiber-optic gyroscope inclinometer
CN110567492A (en) Low-cost MEMS inertial sensor system-level calibration method
CN106767686A (en) Structure sediment monitoring method and system based on inclination angle Yu structural joint deformation measurement
Leppäkoski et al. Error analysis of step length estimation in pedestrian dead reckoning
CN113153270A (en) Measurement-while-drilling method for near-bit dynamic well inclination angle and tool face angle
CN106197376B (en) Car body obliqueness measurement method based on uniaxial MEMS inertial sensor
CN106705995A (en) Calibration method of MEMS gyroscope g value sensitive coefficient
Yang et al. Research on identification and suppression of vibration error for MEMS inertial sensor in near-bit inclinometer
CN104500036B (en) A kind of method and device for being applied to dynamic measurement in hole deviation tool face azimuth under The Small Well slanted bar part
CN204788411U (en) Wall climbing robot navigational positioning system based on sound localization , INS
CN107063597A (en) Bridges Detection based on POS system
CN109212620A (en) Moving base rotating accelerometer gravity gradiometer error-compensating apparatus and method
Zhang et al. High precision tri-axial quartz flexible accelerometers resolution measurement method based on tri-axial turntable
CN109209346A (en) A method of reducing irregular oscillation influences continuous gyroscopic survey instrument precision
Escalona Vertical track geometry monitoring using inertial sensors and complementary filters
CN105444777B (en) A kind of fibre optic gyroscope error testing method under swinging condition

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20190115