CN116816338A - 基于三分量电子罗盘振动监测的随钻岩层识别方法 - Google Patents
基于三分量电子罗盘振动监测的随钻岩层识别方法 Download PDFInfo
- Publication number
- CN116816338A CN116816338A CN202310533501.5A CN202310533501A CN116816338A CN 116816338 A CN116816338 A CN 116816338A CN 202310533501 A CN202310533501 A CN 202310533501A CN 116816338 A CN116816338 A CN 116816338A
- Authority
- CN
- China
- Prior art keywords
- vibration
- rock
- drilling
- signals
- wavelet packet
- 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
Links
- 239000011435 rock Substances 0.000 title claims abstract description 72
- 238000005553 drilling Methods 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000012544 monitoring process Methods 0.000 title claims abstract description 19
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 26
- 238000005070 sampling Methods 0.000 claims abstract description 22
- 230000008569 process Effects 0.000 claims abstract description 12
- 230000001133 acceleration Effects 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 9
- 239000012530 fluid Substances 0.000 claims description 3
- 230000006870 function Effects 0.000 claims description 3
- 230000009131 signaling function Effects 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 description 7
- 238000005755 formation reaction Methods 0.000 description 7
- 230000009471 action Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000010720 hydraulic oil Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 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
- 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/003—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 by analysing drilling variables or conditions
-
- 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
- E21B47/13—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 by electromagnetic energy, e.g. radio frequency
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D21/00—Measuring or testing not otherwise provided for
- G01D21/02—Measuring two or more variables by means not covered by a single other subclass
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- Geochemistry & Mineralogy (AREA)
- Remote Sensing (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Electromagnetism (AREA)
- Geophysics (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Earth Drilling (AREA)
Abstract
本发明公开了一种基于三分量电子罗盘振动监测的随钻岩层识别方法,包括:通过三分量电子罗盘采集钻杆在Y、X和Z方向的振动信号,并将振动信号通过WIFI模块发送给监测主机处理,包括对监测到的振动信号采用短时傅里叶变换进行时频分析;对每个采样时间段tn内的振动信号数据组进行方差、均方根、波峰因素和峭度指数分析;对每个采样时间段tn内的振动信号数据组进行波峰因素和峭度指数分析;对振动信号分别采用小波包分解,计算小波包分解后重构信号的能量。本发明通过对振动信号分析得到的方差、均方根值、波峰因素和峭度等信号特征以及小波包分解后重构信号的能量实现了在钻探过程中对岩层界面与岩体结构特征的识别与区分。
Description
技术领域
本发明涉及地质勘探技术领域,特别涉及一种随钻震源装置。
背景技术
钻孔取芯是地质勘探中常用的勘测手段,取芯完成后,技术人员通过人工编录的方法获取岩体结构特征,最终形成钻孔柱状图和地质勘察报告,这种方法由于岩芯数量多,耗费量时间久,人工编录的工作量相当庞大,需要的技术人员众多,并且在编录过程中,由于不同技术人员的经验不同,常常导致汇编出的结果出现偏差。若能从地质钻探数量庞大的钻孔施工过程中实时快速获得岩体结构特征等信息,无疑对地质勘探能力和水平的提升具有重要的意义。
研究发现,钻机钻头在钻进破岩过程中,钻杆会受到给进压力和钻机扭矩施加的作用力,同时还受到来自于岩层抵抗切削破碎的反作用力,钻机在钻孔作业时,依靠液压油缸提供的给进压力和旋转压力的共同作用破碎岩石。一般情况下,钻进岩层存在不均质性,在这些综合作用力的影响下,钻杆就会产生振动现象。钻杆在钻进不同硬度的岩层时,钻杆振动信号会发生变化,因此对振动信号进行监测和研究可以定性识别岩层的硬度和岩层界面。
发明内容
有鉴于此,本发明的目的是提供一种基于三分量电子罗盘振动监测的随钻岩层识别方法,以解决对钻杆振动信号的监测和利用监测到的钻杆振动信号对岩层界面及岩层硬度进行识别的技术问题。
本发明基于三分量电子罗盘振动监测的随钻岩层识别方法,其特征在于:包括以下步骤:
I)将三分量电子罗盘安装在钻杆顶端的钻井液管接头上,所述三分量电子罗盘包括三向振动加速度传感器、与三向振动加速度传感器连接的控制板、与控制板连接的WIFI模块和与控制板连接的电源模块,通过三向振动加速度传感器采集钻杆在Y方向、X方向和Z方向的振动信号,Y方向为钻杆的轴向,X方向垂直于钻杆轴向,Z方向垂直于Y方向和X方向,控制板对三向振动加速度传感器采集的信号进行处理,并控制WIFI模块将处理后的振动信号发送给监测主机;
II)监测主机对接收到的振动信号进行如下处理:
a、对实时监测到的Y方向、X方向和Z方向的振动信号采用短时傅里叶变换进行时频分析:
式中,S(f,t)为进行短时傅里叶变换后的时频函数;f(τ)为振动信号的幅值;ω(τ-t)为短时分析窗口,τ为有限时间,,j∈Z,f为振幅;然后根据Y方向、X方向和Z方向的振动信号幅值f(τ)和与岩石硬度对应的参考幅值进行比较,定性分辨出当前钻遇岩石的硬度;
b、对每个采样时间段tn内的Y方向、X方向和Z方向的振动信号数据组进行方差、均方根、波峰因素和峭度指数分析,公式如下:
1)方差
2)均方根值
以上公式中N表示时间采样点的数量,xi表示一帧振动数据中第i个采样点的幅值,x表示一帧振动数据的N个采样点的振动幅值平均值,方差和均方根值反应振动信号的平均能量,将计算得到的方差σ和参考方差σref进行比较,将计算得到的均方根xr和参考均方根xref进行比较,若σ>σref且xr>xref,则判断当前采样时间段内钻遇岩层出现界面变化,否则钻遇岩层没有出现界面变化;
3)对每个采样时间段tn内的Y方向、X方向和Z方向的振动信号数据组进行波峰因素和峭度指数分析,公式如下:
波峰因素
峭度
上式中|xi|max表示一帧振动数据的N个采样点中的最大振幅,xrms一帧振动数据的波形有效值,波峰因素表示硬岩和软岩的差异程度,峭度指数表示硬岩和软岩的变化率,通过波峰因素和峭度判断当前钻遇岩层属于硬岩或软岩,以定性地对岩石性质进行分析。
4)对Y方向、X方向和Z方向的振动信号分别采用小波包分解,得到信号在每个频段的分辨率,小波包分解得到新的信号函数如下:
S=S30+S31+S32+S33+S34+S35+S36+S37
上式中S表示原始信号;S30,S31,S32,S33,S34,S35,S36和S37分别表示重构小波包分解系数。
计算小波包分解后重构信号的能量,能量计算公式如下:
上式中Sij(t)为分解后的重构信号,Xjk(t)表示重构信号的幅值,下标i表示小波包分解层数、j表示小波包分解对应的节点、k表示次数,n表示帧数,通过重构信号的能量Eij进行岩层界面识别与岩体结构特征的识别与区分。
本发明的有益效果:
本发明基于三分量电子罗盘振动监测的随钻岩层识别方法,其通过三分量电子罗盘采集钻杆在X、Y、Z三个方向上的振动信号,并对振动信号提取方差、均方根值、波峰因素和峭度等信号特征,还对振动信号进行小波包分解,并计算小波包分解后重构信号的能量,通过得到的方差、均方根值、波峰因素和峭度等信号特征以及小波包分解后重构信号的能量实现了在钻探过程中对岩层硬度及岩层界面进行定性识别。
附图说明
图1为三分量电子罗盘的结构示意图。
具体实施方式
下面结合实施例对本发明作进一步描述。
本实施例中基于三分量电子罗盘振动监测的随钻岩层识别方法包括以下步骤:
I)将三分量电子罗盘安装在钻杆顶端的钻井液管接头上,所述三分量电子罗盘包括三向振动加速度传感器1、与三向振动加速度传感器连接的控制板2、与控制板连接的WIFI模块3和与控制板连接的电源模块4,通过三向振动加速度传感器采集钻杆在Y方向、X方向和Z方向的振动信号,Y方向为钻杆的轴向,X方向垂直于钻杆轴向,Z方向垂直于Y方向和X方向,控制板对三向振动加速度传感器采集的信号进行处理,并控制WIFI模块将处理后的振动信号发送给监测主机5。
II)监测主机对接收到的振动信号进行如下处理:
a、对实时监测到的Y方向、X方向和Z方向的振动信号采用短时傅里叶变换进行时频分析:
式中,S(f,t)为进行短时傅里叶变换后的时频函数;f(τ)为振动信号的幅值;ω(τ-t)为短时分析窗口,τ为有限时间,,j∈Z,f为振幅;然后根据Y方向、X方向和Z方向的振动信号幅值f(τ)和与岩石硬度对应的参考幅值进行比较,定性分辨出当前钻遇岩石的硬度。
b、对每个采样时间段tn内的Y方向、X方向和Z方向的振动信号数据组进行方差、均方根、波峰因素和峭度指数分析,公式如下:
1)方差
2)均方根值
以上公式中N表示时间采样点的数量,xi表示一帧振动数据中第i个采样点的幅值,x表示一帧振动数据的N个采样点的振动幅值平均值,方差和均方根值反应振动信号的平均能量,将计算得到的方差σ和参考方差σref进行比较,将计算得到的均方根xr和参考均方根xref进行比较,若σ>σref且xr>xref,则判断当前采样时间段内钻遇岩层出现界面变化,否则钻遇岩层没有出现界面变化。
3)对每个采样时间段tn内的Y方向、X方向和Z方向的振动信号数据组进行波峰因素和峭度指数分析,公式如下:
波峰因素
峭度
上式中|xi|max表示一帧振动数据的N个采样点中的最大振幅,xrms一帧振动数据的波形有效值,波峰因素表示硬岩和软岩的差异程度,峭度指数表示硬岩和软岩的变化率,通过波峰因素和峭度判断当前钻遇岩层属于硬岩或软岩,以定性地对岩石性质进行分析。
4)对Y方向、X方向和Z方向的振动信号分别采用小波包分解,得到信号在每个频段的分辨率,小波包分解得到新的信号函数如下:
S=S30+S31+S32+S33+S34+S35+S36+S37
上式中S表示原始信号;S30,S31,S32,S33,S34,S35,S36和S37分别表示重构小波包分解系数。
计算小波包分解后重构信号的能量,能量计算公式如下:
上式中Sij(t)为分解后的重构信号,Xjk(t)表示重构信号的幅值,下标i表示小波包分解层数、j表示小波包分解对应的节点、k表示次数,n表示帧数,通过重构信号的能量Eij进行岩层界面识别与岩体结构特征的识别与区分。
本实施例基于三分量电子罗盘振动监测的随钻岩层识别方法,其通过三分量电子罗盘采集钻杆在X、Y、Z三个方向上的振动信号,并对振动信号提取方差、均方根值、波峰因素和峭度等信号特征,还对振动信号进行小波包分解,并计算小波包分解后重构信号的能量,通过得到的方差、均方根值、波峰因素和峭度等信号特征以及小波包分解后重构信号的能量实现了在钻探过程中对岩层硬度及岩层界面进行定性识别。在钻探过程中,可以根据对岩层硬度及岩层界面的定性识别结果,指导改变钻井参数,例如判断钻入硬质岩层时,可以减小进给速度,当判断钻入软质岩层时可以增大钻进速度,有利于提高钻探工作效率。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (1)
1.基于三分量电子罗盘振动监测的随钻岩层识别方法,其特征在于:包括以下步骤:
I)将三分量电子罗盘安装在钻杆顶端的钻井液管接头上,所述三分量电子罗盘包括三向振动加速度传感器、与三向振动加速度传感器连接的控制板、与控制板连接的WIFI模块和与控制板连接的电源模块,通过三向振动加速度传感器采集钻杆在Y方向、X方向和Z方向的振动信号,Y方向为钻杆的轴向,X方向垂直于钻杆轴向,Z方向垂直于Y方向和X方向,控制板对三向振动加速度传感器采集的信号进行处理,并控制WIFI模块将处理后的振动信号发送给监测主机;
II)监测主机对接收到的振动信号进行如下处理:
a、对实时监测到的Y方向、X方向和Z方向的振动信号采用短时傅里叶变换进行时频分析:
式中,S(f,t)为进行短时傅里叶变换后的时频函数;f(τ)为振动信号的幅值;ω(τ-t)为短时分析窗口,τ为有限时间,j∈Z,f为振幅;然后根据Y方向、X方向和Z方向的振动信号幅值f(τ)和与岩石硬度对应的参考幅值进行比较,定性分辨出当前钻遇岩石的硬度;
b、对每个采样时间段tn内的Y方向、X方向和Z方向的振动信号数据组进行方差、均方根、波峰因素和峭度指数分析,公式如下:
1)方差
2)均方根值
以上公式中N表示时间采样点的数量,xi表示一帧振动数据中第i个采样点的幅值,表示一帧振动数据的N个采样点的振动幅值平均值,方差和均方根值反应振动信号的平均能量,将计算得到的方差σ和参考方差σref进行比较,将计算得到的均方根xr和参考均方根xref进行比较,若σ>σref且xr>xref,则判断当前采样时间段内钻遇岩层出现界面变化,否则钻遇岩层没有出现界面变化,用以识别岩层界面;
3)对每个采样时间段tn内的Y方向、X方向和Z方向的振动信号数据组进行波峰因素和峭度指数分析,公式如下:
波峰因素
峭度
上式中|xi|max表示一帧振动数据的N个采样点中的最大振幅,xrms表示一帧振动数据的波形有效值,波峰因素表示硬岩和软岩的差异程度,峭度指数表示硬岩和软岩的变化率,通过波峰因素和峭度判断当前钻遇岩层属于硬岩或软岩,以定性地对岩石性质进行分析;
4)对Y方向、X方向和Z方向的振动信号分别采用小波包分解,得到信号在每个频段的分辨率,小波包分解得到新的信号函数如下:
S=S30+S31+S32+S33+S34+S35+S36+S37
上式中S表示原始信号;S30,S31,S32,S33,S34,S35,S36和S37分别表示重构小波包分解系数;
计算小波包分解后重构信号的能量,能量计算公式如下:
上式中Sij(t)为分解后的重构信号,Xjk(t)表示重构信号的幅值,下标i表示小波包分解层数、j表示小波包分解对应的节点、k表示次数,n表示帧数,通过重构信号的能量Eij进行岩层界面识别与岩体结构特征的识别与区分。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310533501.5A CN116816338A (zh) | 2023-05-12 | 2023-05-12 | 基于三分量电子罗盘振动监测的随钻岩层识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310533501.5A CN116816338A (zh) | 2023-05-12 | 2023-05-12 | 基于三分量电子罗盘振动监测的随钻岩层识别方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116816338A true CN116816338A (zh) | 2023-09-29 |
Family
ID=88121225
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310533501.5A Pending CN116816338A (zh) | 2023-05-12 | 2023-05-12 | 基于三分量电子罗盘振动监测的随钻岩层识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116816338A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117268842A (zh) * | 2023-11-23 | 2023-12-22 | 华侨大学 | 一种岩土工程勘察用取样装置及其使用方法 |
-
2023
- 2023-05-12 CN CN202310533501.5A patent/CN116816338A/zh active Pending
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117268842A (zh) * | 2023-11-23 | 2023-12-22 | 华侨大学 | 一种岩土工程勘察用取样装置及其使用方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110109895B (zh) | 适用于tbm掘进隧道的围岩分级联合预测方法及应用 | |
US10036246B2 (en) | Apparatus, computer readable medium, and program code for evaluating rock properties while drilling using downhole acoustic sensors and a downhole broadband transmitting system | |
US20140334260A1 (en) | Neural Network Signal Processing of Microseismic Events | |
EP2761335B1 (en) | Methods of evaluating rock properties while drilling using downhole acoustic sensors and telemetry system | |
US10180061B2 (en) | Methods of evaluating rock properties while drilling using downhole acoustic sensors and a downhole broadband transmitting system | |
CN101879690B (zh) | 一种基于多传感器信号的批量钻削工序质量监测方法 | |
CN116816338A (zh) | 基于三分量电子罗盘振动监测的随钻岩层识别方法 | |
CN110529087B (zh) | 一种评价地层水力压裂效果的方法及装置 | |
CN1309229A (zh) | 研究和鉴定基础的类型用的方法和装置 | |
EP2761336B1 (en) | Apparatus, computer readable medium, and program code for evaluating rock properties while drilling using downhole acoustic sensors and a downhole broadband transmitting system | |
EP2761337B1 (en) | Methods of evaluating rock properties while drilling using downhole acoustic sensors and a downhole broadband transmitting system | |
CN107678064B (zh) | 一种声波时差实时提取方法 | |
CN112781820B (zh) | 一种滚刀性能退化趋势评估方法 | |
CN107350900A (zh) | 一种基于断屑时间提取的刀具状态监测方法 | |
CN109542884B (zh) | 采用大数据耦合参数以优化破岩效率的方法和系统 | |
CN113390458A (zh) | 一种爆破区围岩损伤程度的判定方法 | |
CN116677367A (zh) | 一种基于随钻参数特征的智能实时岩层反演识别方法 | |
CN104278990B (zh) | 测井数据质量恢复方法及装置 | |
CN115950947A (zh) | 一种tbm隧洞基于振动信号实时的地质参数预测方法 | |
CN110761778B (zh) | 基于声波的压裂裂缝监测方法及系统 | |
CN112412390B (zh) | 基于深度学习模型评价固井第二界面的方法及装置 | |
CN114112651A (zh) | 一种用于人造岩心的岩石动静态力学参数转换方法及系统 | |
CN113075120B (zh) | 一种实时土体类别识别方法、系统及一种盾构机 | |
CN104571090A (zh) | 随钻测井地面机测试系统 | |
Senjoba et al. | Lithology identification during rotary percussion drilling based on acceleration waveform 1D convolutional neural network |
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 |