CN105604541B - 一种生产测井多臂井径斜井校正处理的方法 - Google Patents
一种生产测井多臂井径斜井校正处理的方法 Download PDFInfo
- Publication number
- CN105604541B CN105604541B CN201511001098.3A CN201511001098A CN105604541B CN 105604541 B CN105604541 B CN 105604541B CN 201511001098 A CN201511001098 A CN 201511001098A CN 105604541 B CN105604541 B CN 105604541B
- Authority
- CN
- China
- Prior art keywords
- depth
- data
- relative variation
- arm caliper
- inclined shaft
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000012937 correction Methods 0.000 title claims abstract description 25
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 20
- 238000001914 filtration Methods 0.000 claims abstract description 28
- 238000005259 measurement Methods 0.000 claims description 10
- 101100382321 Caenorhabditis elegans cal-1 gene Proteins 0.000 claims description 4
- 101100096319 Drosophila melanogaster Spc25 gene Proteins 0.000 claims description 4
- 238000012545 processing Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 1
Classifications
-
- 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/08—Measuring diameters or related dimensions at the borehole
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)
- Control Of Metal Rolling (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明公开了一种生产测井多臂井径斜井校正处理的方法,包括以下步骤:1)获取滤波数据B40FilterDepth;2)将原始数据B40Depth与步骤1)得到的滤波数据B40FilterDepth相减,得相对变化量B40ResDepth;3)对相对变化量B40ResDepth的L个井径测量值进行5点中值滤波,然后再进行椭圆拟合,得拟合参数,其中,所述拟合参数包括椭圆的长、短轴及偏心距;4)根据步骤3)得到的拟合参数对相对变化量B40FilterDepth进行椭圆校正,得校正后的相对变化量B40DEVCorrDepth;5)将校正后的相对变化量B40DEVCorrDepth与相对变量B40ResDepth相加,再将相加得到的结果作为生产测井多臂井径斜井的结果。本发明能够实现对生产测井多臂井径斜井进行准确校正。
Description
技术领域
本发明属于石油开发和工程领域,涉及一种生产测井多臂井径斜井校正处理的方法。
背景技术
随着生产测井技术的发展,多臂井径仪器在套管损伤、扭曲变形检测方面有着广泛的应用。每个深度点沿井周测量多个井径值,利用测量的井径值,可以非常直观的反应套管扭曲变形情况。测量时,一般要加扶正器,保证仪器居中,这样测量的数据就不需要做校正处理。
但是在井况较差情况下,仪器测量时无法加上扶正器,测量时就存在偏心情况。另外,斜井和大斜度井越来越多,即使加上扶正器,由于仪器的重量影响,也很难保证仪器居中,在这两种情况需要进行校正处理,以反映真实的套管变形情况。
但是套管没有变形、而是由于仪器偏心测量导致的“套管变形”和套管扭曲变形、仪器居中测量得到的真实的“套管变形”是很难区分的,如果不加区别的进行椭圆拟合校正处理,会将套管扭曲变形导致的偏心进行误校正,不利于其后的套损解释。
发明内容
本发明的目的在于克服上述现有技术的缺点,提供了一种生产测井多臂井径斜井校正处理的方法,该方法能够实现对生产测井多臂井径斜井进行准确校正。
为达到上述目的,本发明所述的生产测井多臂井径斜井校正处理的方法包括以下步骤:
1)当输入标示FLAG=1时,则对多臂井径B40当前处理深度段StartDepth:EndDepth之间的每个深度点数据利用上下N个深度点数据进行平滑滤波,得到滤波数据B40FilterDepth;
当输入标示FLAG=2时,对多臂井径B40当前处理深度段StartDepth:EndDepth之间的深度点数据利用多臂井径B40StartDepth和多臂井径B40EndDepth进行线性插值,然后对线性插值的结果利用上下N个深度点数据进行平滑滤波,得到滤波数据B40FilterDepth;
2)将原始数据B40Depth与步骤1)得到的滤波数据B40FilterDepth相减,得相对变化量B40ResDepth;
3)对相对变化量B40ResDepth的L个井径测量值进行5点中值滤波,然后再进行椭圆拟合,得拟合参数,其中,所述拟合参数包括椭圆的长、短轴及偏心距;
4)根据步骤3)得到的拟合参数对相对变化量B40FilterDepth进行椭圆校正,得校正后的相对变化量B40DEVCorrDepth;
5)将校正后的相对变化量B40DEVCorrDepth与相对变量B40ResDepth相加,再将相加得到的结果作为生产测井多臂井径斜井的结果。
步骤1中,当处理井段套管扭曲变形小于预设值,且仅存在偏心,则输入标示FLAG=1。
用上下N个深度点数据对当前深度点数据进行N点平滑滤波的具体操作为:每个深度点数据都有共Cal0,Cal1,Cal2.....CalL-1L个测量井径值,对N个深度点数据的N个Calj进行平滑滤波,其中,0≤j≤L-1,
其中,Calji为第i个深度点数值的Calj。
当套管扭曲变形的幅度大于等于预设值时,则输入标示FLAG=2。
设多臂井径仪器的井径数据为L点的井径数据(Ri,θi),i=0,2,....L-1,则将所述L点的井径数据(Ri,θi)利用下式(3)转为直角坐标系:
xi=Ri*Cosθi,yi=Ri*Sinθi (3)
设椭圆方程
其中,A、B、C、D及E为椭圆的方程系数;
根据最小二乘原理得方程组:
求解式(5),得到椭圆的方程系数(A,B,C,D,E),然后根据式(6)、式(7)、式(8)得椭圆的中心(x0,y0)、长轴a、短轴b及偏心距,其中,
则校正后的相对变化量B40DEVCorrDepth为:
本发明具有以下有益效果:
本发明所述的生产测井多臂井径斜井校正处理的方法在具体操作时,通过插值、滤波后的滤波数据B40FilterDepth与原始数据B40Depth进行相减,得相对变化量B40ResDepth,通过所述相对变化量B40ResDepth反映真实的套管值变化情况,再通过椭圆拟合校正对相对变化量B40ResDepth进行校正,然后再根据校正后的相对变化量B40DEVCorrDepth与相对变量B40ResDepth进行相加,相加的结果作为生产测井多臂井径斜井的结果,其中,通过椭圆拟合校正的方法对多臂井径仪器偏心测量的数据进行校正,解决了多臂井径仪器偏心测量数据无法校正的问题,扩大了多臂井径仪器的应用范围。
附图说明
图1为本发明的流程图。
具体实施方式
下面结合附图对本发明做进一步详细描述:
参考图1,本发明所述的生产测井多臂井径斜井校正处理的方法包括以下步骤:
1)当输入标示FLAG=1时,则对多臂井径B40当前处理深度段StartDepth:EndDepth之间的每个深度点数据利用上下N个深度点数据进行平滑滤波,得到滤波数据B40FilterDepth;
当输入标示FLAG=2时,对多臂井径B40当前处理深度段StartDepth:EndDepth之间的深度点数据利用多臂井径B40StartDepth和多臂井径B40EndDepth进行线性插值,然后对线性插值的结果利用上下N个深度点数据进行平滑滤波,得到滤波数据B40FilterDepth;
2)将原始数据B40Depth与步骤1)得到的滤波数据B40FilterDepth相减,得相对变化量B40ResDepth;
3)对相对变化量B40ResDepth的L个井径测量值进行5点中值滤波,然后再进行椭圆拟合,得拟合参数,其中,所述拟合参数包括椭圆的长、短轴及偏心距;
4)根据步骤3)得到的拟合参数对相对变化量B40FilterDepth进行椭圆校正,得校正后的相对变化量B40DEVCorrDepth;
5)将校正后的相对变化量B40DEVCorrDepth与相对变量B40ResDepth相加,再将相加得到的结果作为生产测井多臂井径斜井的结果。
当处理井段套管扭曲变形小于预设值,且仅存在偏心,则输入标示FLAG=1。
用上下N个深度点数据对当前深度点数据进行N点平滑滤波的具体操作为:每个深度点数据都有共Cal0,Cal1,Cal2.....CalL-1L个测量井径值,对N个深度点数据的N个Calj进行平滑滤波,其中,0≤j≤L-1,
其中,Calji为第i个深度点数值的Calj。
当套管扭曲变形的幅度大于等于预设值时,则输入标示FLAG=2,对多臂井径B40当前处理深度段StartDepth:EndDepth之间的深度点数据,利用B40StartDepth和B40EndDepth进行线性插值的具体操作为:对Cal0,Cal1,Cal2.....CalL-1分别进行线性插值;
设多臂井径仪器的井径数据为L点的井径数据(Ri,θi),i=0,2,....L-1,则将所述L点的井径数据(Ri,θi)利用下式(3)转为直角坐标系:
xi=Ri*Cosθi,yi=Ri*Sinθi (3)
设椭圆方程
其中,A、B、C、D及E为椭圆的方程系数;
根据最小二乘原理得方程组:
求解式(5),得到椭圆的方程系数(A,B,C,D,E),然后根据式(6)、式(7)、式(8)得椭圆的中心(x0,y0)、长轴a、短轴b及偏心距,其中,
则校正后的相对变化量B40DEVCorrDepth为:
以上对本发明实施例所提供的实施方案进行了详细介绍,本文中应用了具体个例对本发明实施例的原理以及实施方式进行了阐述,以上实施例的说明只适用于帮助理解本发明实施例的原理;同时,对于本领域的一般技术人员,依据本发明实施例,在具体实施方式以及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (4)
1.一种生产测井多臂井径斜井校正处理的方法,其特征在于,包括以下步骤:
1)当输入标示FLAG=1时,则对多臂井径B40当前处理深度段StartDepth:EndDepth之间的每个深度点数据利用上下N个深度点数据进行平滑滤波,得到滤波数据B40FilterDepth;
当输入标示FLAG=2时,对多臂井径B40当前处理深度段StartDepth:EndDepth之间的深度点数据利用多臂井径B40StartDepth和多臂井径B40EndDepth进行线性插值,然后对线性插值的结果利用上下N个深度点数据进行平滑滤波,得到滤波数据B40FilterDepth;
2)将原始数据B40Depth与步骤1)得到的滤波数据B40FilterDepth相减,得相对变化量B40ResDepth;
3)对相对变化量B40ResDepth的L个井径测量值进行5点中值滤波,然后再进行椭圆拟合,得拟合参数,其中,所述拟合参数包括椭圆的长、短轴及偏心距;
4)根据步骤3)得到的拟合参数对相对变化量B40FilterDepth进行椭圆校正,得校正后的相对变化量B40DEVCorrDepth;
5)将校正后的相对变化量B40DEVCorrDepth与相对变化量B40ResDepth相加,再将相加得到的结果作为生产测井多臂井径斜井的结果;
步骤1中,当处理井段套管扭曲变形小于预设值,且仅存在偏心,则输入标示FLAG=1;
当套管扭曲变形的幅度大于等于预设值时,则输入标示FLAG=2。
2.根据权利要求1所述的生产测井多臂井径斜井校正处理的方法,其特征在于,用上下N个深度点数据对当前深度点数据进行N点平滑滤波的具体操作为:每个深度点数据都有共Cal0,Cal1,Cal2.....CalL-1L个测量井径值,对N个深度点数据的N个Calj进行平滑滤波,其中,0≤j≤L-1,
其中,Calji为第i个深度点数值的Calj。
3.根据权利要求1所述的生产测井多臂井径斜井校正处理的方法,其特征在于,
设多臂井径仪器的井径数据为L点的井径数据(Ri,θi),i=0,2,....L-1,则将所述L点的井径数据(Ri,θi)利用下式(3)转为直角坐标系:
xi=Ri*Cosθi,yi=Ri*Sinθi (3)
设椭圆方程
其中,A、B、C、D及E为椭圆的方程系数;
根据最小二乘原理得方程组:
求解式(5),得到椭圆的方程系数(A,B,C,D,E),然后根据式(6)、式(7)、式(8)得椭圆的中心(x0,y0)、长轴a、短轴b及偏心距,其中,
4.根据权利要求3所述的生产测井多臂井径斜井校正处理的方法,其特征在于,校正后的相对变化量B40DEVCorrDepth为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511001098.3A CN105604541B (zh) | 2015-12-28 | 2015-12-28 | 一种生产测井多臂井径斜井校正处理的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511001098.3A CN105604541B (zh) | 2015-12-28 | 2015-12-28 | 一种生产测井多臂井径斜井校正处理的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105604541A CN105604541A (zh) | 2016-05-25 |
CN105604541B true CN105604541B (zh) | 2018-11-16 |
Family
ID=55984826
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511001098.3A Active CN105604541B (zh) | 2015-12-28 | 2015-12-28 | 一种生产测井多臂井径斜井校正处理的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105604541B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107299826B (zh) * | 2016-04-14 | 2019-02-15 | 中国石油天然气股份有限公司 | 表层套管固井居中方法 |
CN107780926B (zh) * | 2017-09-20 | 2020-05-22 | 南方科技大学 | 井眼截面形状测量方法、装置和终端 |
CN109372497B (zh) * | 2018-08-20 | 2022-03-29 | 中国石油天然气集团有限公司 | 一种超声成像动态均衡处理的方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4803667A (en) * | 1986-11-13 | 1989-02-07 | Atlantic Richfield Company | Televiewer processing system |
US5737277A (en) * | 1996-08-01 | 1998-04-07 | Western Atlas International, Inc. | Method for computing borehole geometry from ultrasonic pulse echo data |
CA2336039C (en) * | 1998-06-26 | 2004-09-14 | Dresser Industries, Inc. | Method and apparatus for quick determination of the ellipticity of an earth borehole |
GB2415257B (en) * | 2004-06-18 | 2008-10-08 | Pathfinder Energy Services Inc | Estimation of borehole geometry parameters and lateral tool displacements |
GB2473561A (en) * | 2008-06-11 | 2011-03-16 | Baker Hughes Inc | Multi-resolution borehole profiling |
CN103225502A (zh) * | 2013-04-16 | 2013-07-31 | 中国石油天然气集团公司 | 一种测量井径的系统及其方法 |
CA2947596A1 (en) * | 2014-05-09 | 2015-11-12 | Probe Holdings, Inc. | Caliper tool with in-situ temperature compensation |
-
2015
- 2015-12-28 CN CN201511001098.3A patent/CN105604541B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4803667A (en) * | 1986-11-13 | 1989-02-07 | Atlantic Richfield Company | Televiewer processing system |
US5737277A (en) * | 1996-08-01 | 1998-04-07 | Western Atlas International, Inc. | Method for computing borehole geometry from ultrasonic pulse echo data |
CA2336039C (en) * | 1998-06-26 | 2004-09-14 | Dresser Industries, Inc. | Method and apparatus for quick determination of the ellipticity of an earth borehole |
GB2415257B (en) * | 2004-06-18 | 2008-10-08 | Pathfinder Energy Services Inc | Estimation of borehole geometry parameters and lateral tool displacements |
GB2473561A (en) * | 2008-06-11 | 2011-03-16 | Baker Hughes Inc | Multi-resolution borehole profiling |
CN103225502A (zh) * | 2013-04-16 | 2013-07-31 | 中国石油天然气集团公司 | 一种测量井径的系统及其方法 |
CA2947596A1 (en) * | 2014-05-09 | 2015-11-12 | Probe Holdings, Inc. | Caliper tool with in-situ temperature compensation |
Non-Patent Citations (3)
Title |
---|
基于最小二乘椭圆拟合的井径测量面偏角修正方法;刘春光等;《机械制造》;20150531;第53卷(第609期);第76-79页 * |
基于最小二乘法的成像测井中井径偏心图像校正方法研究;张健;《科学技术与工程》;20140228;第14卷(第6期);第77-81页 * |
超声成像测井井径的椭圆拟合及偏心时间图像的校正;余厚全等;《测井技术》;20111031;第35卷(第5期);第418-421页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105604541A (zh) | 2016-05-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105604541B (zh) | 一种生产测井多臂井径斜井校正处理的方法 | |
SA516380312B1 (ar) | معايرة ذاتية لمستشعر توجيه مقيد في الوقت الفعلي | |
KR101783050B1 (ko) | 3축 자계 센서의 교정 방법 및 장치 | |
CN105823504B (zh) | 一种编码器的越零点处理方法 | |
US20110010098A1 (en) | Method of pre-stack two-dimension-like transformation of three-dimensional seismic record | |
CN109142792B (zh) | 一种石英挠性加速度计温度误差标定补偿方法 | |
CN110544255B (zh) | 面向3d打印的三角网格孔洞的分割方法及系统 | |
Choudhury et al. | Tripled coincidence point theorems for nonlinear contractions in partially ordered metric spaces | |
CN107270943B (zh) | 一种相关两维的惯性导航落点精度评估方法 | |
CN104122440B (zh) | 晶振频率漂移感知方法 | |
CN104903750A (zh) | 一种电磁数据的保幅静校正方法 | |
CN105528735B (zh) | 基于量测风速与空间相关性的异常数据点的校正方法 | |
EP3287907A3 (en) | Mobile machine-generated data correction | |
CN105300280A (zh) | 连接器尺寸视觉测量方法 | |
CN105336637A (zh) | 量测晶圆变形的方法 | |
CN205433693U (zh) | 髋臼外杯的检测卡尺结构 | |
CN103604404A (zh) | 基于数值积分的加速度信号测取位移方法 | |
CN104933205A (zh) | 土地利用空间数据处理中基于地理本体的属性匹配方法 | |
CN110954041B (zh) | 三维螺纹测量仪的校准方法 | |
CN106209092A (zh) | 一种基于gps秒脉冲信号的精确数字分频方法及装置 | |
CN106556807A (zh) | 基于内部偏置场和传感方向的磁阻静态特性优化方法 | |
Gupta et al. | A new type of tripled fixed point theorem in partially ordered complete metric space | |
Zhang | A note on ternary cyclotomic polynomials | |
CN103839851A (zh) | 终点判断方法 | |
CN110778306A (zh) | 一种气井产能方程异常的修正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |