CN112228046A - 一种随钻超声井径数据校正方法 - Google Patents
一种随钻超声井径数据校正方法 Download PDFInfo
- Publication number
- CN112228046A CN112228046A CN202011298787.6A CN202011298787A CN112228046A CN 112228046 A CN112228046 A CN 112228046A CN 202011298787 A CN202011298787 A CN 202011298787A CN 112228046 A CN112228046 A CN 112228046A
- Authority
- CN
- China
- Prior art keywords
- borehole diameter
- ultrasonic
- drilling
- borehole
- epsilon
- 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
Links
- 238000005553 drilling Methods 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000012937 correction Methods 0.000 title claims abstract description 9
- 238000005259 measurement Methods 0.000 claims abstract description 14
- 238000013461 design Methods 0.000 claims abstract description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 230000005284 excitation Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 239000012530 fluid Substances 0.000 claims description 4
- 239000011435 rock Substances 0.000 claims description 4
- 238000004088 simulation Methods 0.000 description 6
- 238000001914 filtration Methods 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 3
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP 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
- E21B47/14—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 using acoustic waves
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Theoretical Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics (AREA)
- Remote Sensing (AREA)
- Geochemistry & Mineralogy (AREA)
- Computer Hardware Design (AREA)
- Fluid Mechanics (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开了一种随钻超声井径数据校正方法，属于石油天然气钻探技术领域，其特征在于，包括以下步骤：a、采集超声井径波形；b、计算得到井径r，c、进行归一化，得到实际衰减系数k；d、建立物理模型，计算理论衰减系数k1；e、根据随钻超声井径仪器的设计误差ε，得到实际测量的置信区间[k‑ε，k+ε]；f、将理论衰减系数k1与置信区间进行对比；对深度点井径数据进行校正；g、遍历整个深度区间。本发明不仅可以评价数据的可靠性，也可以针对可靠性低的数据进行校正，能够极大提高超声井径数据的可靠性，为井孔稳定性评价提供保障。
Description
技术领域
本发明涉及到石油天然气钻探技术领域，尤其涉及一种随钻超声井径数据校正方法。
背景技术
井径仪的种类很多，测量原理不尽相同。主要有电阻式、感应式及声波井径仪。井径仪的结构主要有两种：一种是进行单独井径测量的张臂式井径仪；另一种就是利用某些测井仪器的推靠臂，如密度仪、井壁中子测井仪和微侧向仪，在这些仪器测井的同时测量。随着钻井技术的发展，在随钻过程中实时获取井径参数对于保障钻井安全至关重要。
目前井径的计算方法主要有包络线法和多次回波计算法。这两种方法均受岩屑散射干扰较大，在有岩屑信号干扰时不能准确获取井孔的真实井径。此外，测井系统发射的超声波信号经过井液到达井壁的过程中不可避免的受泥浆气泡以及仪器偏心的影响，这种影响将会使井径测量结果精度大大降低，从而使得测得的实时钻井井径值的可靠性降低。
公开号为CN 105604541A，公开日为2016年05月25日的中国专利文献公开了一种生产测井多臂井径斜井校正处理的方法，其特征在于，包括以下步骤：
1)当输入标示FLAG＝1时，则对多臂井径B40当前处理深度段 StartDepth:EndDepth之间的每个深度点数据利用上下N个深度点数据进行平滑滤波，得到滤波数据B40Filter；
当输入标示FLAG＝2时，对多臂井径B40当前处理深度段 StartDepth:EndDepth之间的深度点数据利用多臂井径B40和多臂井径B40 进行线性插值，然后对线性插值的结果利用上下N个深度点数据进行平滑滤波，得到滤波数据B40Filter；
2)将原始数据B40与步骤1)得到的滤波数据B40Filter相减，得相对变化量B40Res；
3)对相对变化量B40Res的L个井径测量值进行5点中值滤波，然后再进行椭圆拟合，得拟合参数，其中，所述拟合参数包括椭圆的长、短轴及偏心距；
4)根据步骤3)得到的拟合参数对相对变化量B40Filter进行椭圆校正，得校正后的相对变化量B40DEVCorr；
5)将校正后的相对变化量B40DEVCorr与相对变化量B40Res相加，再将相加得到的结果作为生产测井多臂井径斜井的结果；
步骤1中，当处理井段套管扭曲变形小于预设值，且仅存在偏心，则输入标示FLAG＝1；
当套管扭曲变形的幅度大于等于预设值时，则输入标示FLAG＝2。
该专利文献公开的井径斜井校正处理的方法，虽然能够实现对生产测井多臂井径斜井进行校正。但是，仍然存在校正精度低，评价数据可靠性差的问题。
发明内容
本发明为了克服上述现有技术的缺陷，提供一种随钻超声井径数据校正方法，本发明基于实测的数据建立物理模型，利用有限差分进行数值模拟，将实测衰减系数与数值模拟结果进行对比，不仅可以评价数据的可靠性，也可以针对可靠性低的数据进行校正，能够极大提高超声井径数据的可靠性，为井孔稳定性评价提供保障。
本发明通过下述技术方案实现：
一种随钻超声井径数据校正方法，其特征在于，包括以下步骤：
a、在一深度区间利用随钻超声井径仪器进行井径测量，采集超声井径波形；
b、在一深度点通过能量包络法提取井壁反射波到时t及井壁反射波能量E1，根据井壁反射波到时t和超声波在井液中传播波速v计算得到井径 r，
c、利用发射波能量E2对井壁反射波能量E1进行归一化，得到实际衰减系数k，
k＝E1/E2 式2
d、根据实际随钻井径测井时的激发频率、泥浆参数、地层参数、仪器尺寸和井径r建立物理模型，利用有限差分数值模拟理论条件下随钻超声井径波形，并根据步骤c中的方法计算理论衰减系数k1；
e、根据随钻超声井径仪器的设计误差ε，得到实际测量的置信区间 [k-ε,k+ε]；
f、将步骤d中计算得到的理论衰减系数k1与步骤e中的置信区间进行对比；若k1在[k-ε,k+ε]范围内，则随钻超声井径仪器测得的井径数据可靠；若k1不在[k-ε,k+ε]范围内，则随钻超声井径仪器测得的井径数据不可靠，考虑为泥浆岩屑反射信号，则在步骤b包络中重新选定反射信号，重复步骤c-步骤f，对深度点井径数据进行校正；
g、重复步骤b-步骤f，直至遍历整个深度区间。
所述步骤d中，泥浆参数包括泥浆速度和泥浆密度，地层参数包括纵波速度、横波速度和密度。
本发明的有益效果是：
基于实测的数据建立物理模型，利用有限差分进行数值模拟，将实测衰减系数与数值模拟结果进行对比，不仅可以评价数据的可靠性，也可以针对可靠性低的数据进行校正，能够极大提高超声井径数据的可靠性，为井孔稳定性评价提供保障。
附图说明
下面将结合说明书附图和具体实施方式对本发明作进一步的具体说明：
图1为本发明正演模拟得到的不同激发频率下衰减系数随井径变化的示意图。
具体实施方式
实施例1
一种随钻超声井径数据校正方法，包括以下步骤：
a、在一深度区间利用随钻超声井径仪器进行井径测量，采集超声井径波形；
b、在一深度点通过能量包络法提取井壁反射波到时t及井壁反射波能量E1，根据井壁反射波到时t和超声波在井液中传播波速v计算得到井径 r，
c、利用发射波能量E2对井壁反射波能量E1进行归一化，得到实际衰减系数k；
k＝E1/E2 式2
d、根据实际随钻井径测井时的激发频率、泥浆参数、地层参数、仪器尺寸和井径r建立物理模型，利用有限差分数值模拟理论条件下随钻超声井径波形，并根据步骤c中的方法计算理论衰减系数k1；
e、根据随钻超声井径仪器的设计误差ε，得到实际测量的置信区间 [k-ε,k+ε]；
f、将步骤d中计算得到的理论衰减系数k1与步骤e中的置信区间进行对比；若k1在[k-ε,k+ε]范围内，则随钻超声井径仪器测得的井径数据可靠；若k1不在[k-ε,k+ε]范围内，则随钻超声井径仪器测得的井径数据不可靠，考虑为泥浆岩屑反射信号，则在步骤b包络中重新选定反射信号，重复步骤c-步骤f，对深度点井径数据进行校正；
g、重复步骤b-步骤f，直至遍历整个深度区间。
实施例2
一种随钻超声井径数据校正方法，包括以下步骤：
a、在一深度区间利用随钻超声井径仪器进行井径测量，采集超声井径波形；
b、在一深度点通过能量包络法提取井壁反射波到时t及井壁反射波能量E1，根据井壁反射波到时t和超声波在井液中传播波速v计算得到井径 r，
c、利用发射波能量E2对井壁反射波能量E1进行归一化，得到实际衰减系数k，
k＝E1/E2 式2
d、根据实际随钻井径测井时的激发频率、泥浆参数、地层参数、仪器尺寸和井径r建立物理模型，利用有限差分数值模拟理论条件下随钻超声井径波形，并根据步骤c中的方法计算理论衰减系数k1；
e、根据随钻超声井径仪器的设计误差ε，得到实际测量的置信区间 [k-ε,k+ε]；
f、将步骤d中计算得到的理论衰减系数k1与步骤e中的置信区间进行对比；若k1在[k-ε,k+ε]范围内，则随钻超声井径仪器测得的井径数据可靠；若k1不在[k-ε,k+ε]范围内，则随钻超声井径仪器测得的井径数据不可靠，考虑为泥浆岩屑反射信号，则在步骤b包络中重新选定反射信号，重复步骤c-步骤f，对深度点井径数据进行校正；
g、重复步骤b-步骤f，直至遍历整个深度区间。
所述步骤d中，泥浆参数包括泥浆速度和泥浆密度，地层参数包括纵波速度、横波速度和密度。
基于实测的数据建立物理模型，利用有限差分进行数值模拟，将实测衰减系数与数值模拟结果进行对比，不仅可以评价数据的可靠性，也可以针对可靠性低的数据进行校正，能够极大提高超声井径数据的可靠性，为井孔稳定性评价提供保障。
下面结合图1来说明本发明的实施案例。针对某一型号仪器，泥浆参数和地层参数如表1所示，设定不同井径、不同频率进行正演模拟，所得的衰减系数如图1所示，从图1中可以看出在井径一定时，随着频率的升高衰减系数逐渐降低；频率一定时，随着井径的增大，衰减系数逐渐降低。在实际运用中，仪器的激发频率为200kHz，测量的井径为0.14m，仪器的测量误差为±5％，实际衰减系数k为0.08，理论衰减系数k1为0.095，此时根据步骤f可判断此时的井径数据不可靠，需要重新选定回波的起跳点。
| 模型参数 | 纵波速度(m/s) | 横波速度(m/s) | 密度(kg/m<sup>3</sup>) |
| 地层参数 | 4000 | 2300 | 2500 |
| 泥浆参数 | 1480 | 0 | 1100 |
表1 。
Claims (2)
1.一种随钻超声井径数据校正方法，其特征在于，包括以下步骤：
a、在一深度区间利用随钻超声井径仪器进行井径测量，采集超声井径波形；
b、在一深度点通过能量包络法提取井壁反射波到时t及井壁反射波能量E1，根据井壁反射波到时t和超声波在井液中传播波速v计算得到井径r，
c、利用发射波能量E2对井壁反射波能量E1进行归一化，得到实际衰减系数k，
k＝E1/E2 式2
d、根据实际随钻井径测井时的激发频率、泥浆参数、地层参数、仪器尺寸和井径r建立物理模型，利用有限差分数值模拟理论条件下随钻超声井径波形，并根据步骤c中的方法计算理论衰减系数k1；
e、根据随钻超声井径仪器的设计误差ε，得到实际测量的置信区间[k-ε，k+ε]；
f、将步骤d中计算得到的理论衰减系数k1与步骤e中的置信区间进行对比；若k1在[k-ε，k+ε]范围内，则随钻超声井径仪器测得的井径数据可靠；若k1不在[k-ε，k+ε]范围内，则随钻超声井径仪器测得的井径数据不可靠，考虑为泥浆岩屑反射信号，则在步骤b包络中重新选定反射信号，重复步骤c-步骤f，对深度点井径数据进行校正；
g、重复步骤b-步骤f，直至遍历整个深度区间。
2.根据权利要求1所述的一种随钻超声井径数据校正方法，其特征在于：所述步骤d中，泥浆参数包括泥浆速度和泥浆密度，地层参数包括纵波速度、横波速度和密度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011298787.6A CN112228046B (zh) | 2020-11-18 | 2020-11-18 | 一种随钻超声井径数据校正方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011298787.6A CN112228046B (zh) | 2020-11-18 | 2020-11-18 | 一种随钻超声井径数据校正方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112228046A true CN112228046A (zh) | 2021-01-15 |
| CN112228046B CN112228046B (zh) | 2024-04-26 |
Family
ID=74123654
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011298787.6A Active CN112228046B (zh) | 2020-11-18 | 2020-11-18 | 一种随钻超声井径数据校正方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112228046B (zh) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004059126A1 (en) * | 2002-12-31 | 2004-07-15 | Services Petroliers Schlumberger | Methods and apparatus for ultrasound velocity measurements in drilling fluids |
| CN102454399A (zh) * | 2010-10-26 | 2012-05-16 | 中国石油化工股份有限公司 | 测井声波时差信号校正方法 |
| CN206000534U (zh) * | 2016-09-12 | 2017-03-08 | 中国海洋石油总公司 | 一种随钻超声井径测量仪器 |
| CN107092031A (zh) * | 2017-06-01 | 2017-08-25 | 中国石油天然气集团公司 | 一种声波测井数据校正方法及装置 |
| CN108571317A (zh) * | 2017-03-09 | 2018-09-25 | 中国石油化工股份有限公司 | 一种随钻测井井径测量系统 |
| CN109765622A (zh) * | 2019-03-06 | 2019-05-17 | 北京华瑞九州能源科技有限公司 | 一种基于井径测量仪的扇区测量数据校正方法 |
| CN110173259A (zh) * | 2019-05-14 | 2019-08-27 | 中石化石油工程技术服务有限公司 | 一种斯通利波能量补偿校正的方法 |
| WO2020036596A1 (en) * | 2018-08-14 | 2020-02-20 | Halliburton Energy Services, Inc. | Eccentricity correction algorithm for borehole shape and tool location computations from caliper data |
| CN111119869A (zh) * | 2018-11-01 | 2020-05-08 | 中国石油化工股份有限公司 | 超声换能器检测装置及随钻超声换能器检测方法 |
-
2020
- 2020-11-18 CN CN202011298787.6A patent/CN112228046B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004059126A1 (en) * | 2002-12-31 | 2004-07-15 | Services Petroliers Schlumberger | Methods and apparatus for ultrasound velocity measurements in drilling fluids |
| CN102454399A (zh) * | 2010-10-26 | 2012-05-16 | 中国石油化工股份有限公司 | 测井声波时差信号校正方法 |
| CN206000534U (zh) * | 2016-09-12 | 2017-03-08 | 中国海洋石油总公司 | 一种随钻超声井径测量仪器 |
| CN108571317A (zh) * | 2017-03-09 | 2018-09-25 | 中国石油化工股份有限公司 | 一种随钻测井井径测量系统 |
| CN107092031A (zh) * | 2017-06-01 | 2017-08-25 | 中国石油天然气集团公司 | 一种声波测井数据校正方法及装置 |
| WO2020036596A1 (en) * | 2018-08-14 | 2020-02-20 | Halliburton Energy Services, Inc. | Eccentricity correction algorithm for borehole shape and tool location computations from caliper data |
| CN111119869A (zh) * | 2018-11-01 | 2020-05-08 | 中国石油化工股份有限公司 | 超声换能器检测装置及随钻超声换能器检测方法 |
| CN109765622A (zh) * | 2019-03-06 | 2019-05-17 | 北京华瑞九州能源科技有限公司 | 一种基于井径测量仪的扇区测量数据校正方法 |
| CN110173259A (zh) * | 2019-05-14 | 2019-08-27 | 中石化石油工程技术服务有限公司 | 一种斯通利波能量补偿校正的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112228046B (zh) | 2024-04-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104345346B (zh) | 一种获取裂缝宽度的方法 | |
| CN110348135B (zh) | 一种随钻声波测井评价地层渗透率的方法 | |
| CN111980676B (zh) | 一种阵列声波测井评价固井质量的方法及处理装置 | |
| EP2887101A1 (en) | Systems and methods for removing coherent noise in log data | |
| CN112593922B (zh) | 一种阵列声波测井评价固井二界面胶结质量的方法及装置 | |
| CN109958432B (zh) | 利用超声回波测井评价固井ⅱ界面胶结质量方法和装置 | |
| CN106548456B (zh) | 校正超声图像轮廓中的偏心的方法 | |
| CN107506556B (zh) | 一种确定新鲜完整岩块声波纵波速度值的方法 | |
| CN114233277B (zh) | 基于套管井模式波衰减系数的固井质量评价方法 | |
| CN106770668B (zh) | 一种用于单孔的基桩质量声波透射法检测方法 | |
| CN112034257A (zh) | 一种井下电阻率的计算方法 | |
| CN106767583B (zh) | 用于基桩检测声波透射法的纵向剖面等效桩径计算方法 | |
| CN111520126A (zh) | 密度测井扩径校正方法及装置 | |
| CN106383173A (zh) | 一种水泥声阻抗计算方法和装置 | |
| CN112228046A (zh) | 一种随钻超声井径数据校正方法 | |
| CN110273675B (zh) | 一种瞬变电磁差分测井方法及其系统 | |
| CN114740532B (zh) | 一种水泥声阻抗、泥浆声阻抗的计算方法 | |
| CN113982562B (zh) | 一种基于阵列声波测井的固井质量评价方法 | |
| CN106761666A (zh) | 一种四探头散射伽马测井及非线性数据反演的方法及装置 | |
| CN111089904B (zh) | 一种考虑地层特性的室内测量波速围压、频散校正方法 | |
| CN113189672B (zh) | 一种基于多属性反演的隧道超前地质预报方法 | |
| CN112835124B (zh) | 基于成像测井和阵列声波测井资料的裂缝有效性评价方法 | |
| CN111456726B (zh) | 声波时差校正方法及装置 | |
| CN109061737B (zh) | 一种基于合成地震记录的储层预测方法及装置 | |
| CN110219642B (zh) | 基于声波传播路径的声波时差校正方法 |
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 |