CN110656933B - 一种破碎带裂隙密度确定方法及系统 - Google Patents
一种破碎带裂隙密度确定方法及系统 Download PDFInfo
- Publication number
- CN110656933B CN110656933B CN201910948240.7A CN201910948240A CN110656933B CN 110656933 B CN110656933 B CN 110656933B CN 201910948240 A CN201910948240 A CN 201910948240A CN 110656933 B CN110656933 B CN 110656933B
- Authority
- CN
- China
- Prior art keywords
- curve
- determining
- logging
- fracture
- fractured zone
- 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
- 238000000034 method Methods 0.000 title claims abstract description 48
- 230000004927 fusion Effects 0.000 claims abstract description 69
- 238000005516 engineering process Methods 0.000 claims abstract description 14
- 238000003384 imaging method Methods 0.000 claims abstract description 13
- 238000005553 drilling Methods 0.000 claims abstract description 10
- 238000000354 decomposition reaction Methods 0.000 claims description 46
- 238000007781 pre-processing Methods 0.000 claims description 15
- 238000010606 normalization Methods 0.000 claims description 13
- 230000035945 sensitivity Effects 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 4
- 230000007613 environmental effect Effects 0.000 claims description 4
- 238000004519 manufacturing process Methods 0.000 abstract description 5
- 239000003245 coal Substances 0.000 abstract description 4
- 239000002184 metal Substances 0.000 abstract description 4
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009533 lab test Methods 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
Images
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- Dispersion Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开一种破碎带裂隙密度确定方法及系统,该方法首先获取野外的测井曲线;根据测井曲线确定敏感测井曲线;其次利用测井曲线融合技术根据敏感测井曲线确定裂隙融合曲线;再次采用离散小波变换方法根据裂隙融合曲线确定破碎带层段;然后在破碎带层段,依据标准钻孔的成像测井资料和敏感测井曲线,建立裂隙密度关系式;最后根据裂隙密度关系式确定待检测钻孔破碎带的裂隙密度。本发明综合利用测井曲线融合技术和离散小波变换方法确定破碎带层段,依据标准钻孔的成像测井资料,建立测井曲线与裂隙密度之间关系,以此确定待检测钻孔破碎带的裂隙密度,不仅能节省成本和时间,还能满足油田、煤田以及金属矿等野外实际生产的需要。
Description
技术领域
本发明涉及地球物理测井技术领域,特别是涉及一种破碎带裂隙密度确定方法及系统。
背景技术
目前,地层破碎带估算裂隙密度的常用方法是根据岩心实验室测试和特殊测井方法(如成像测井、地层倾角测井等),综合利用岩心实验室测试和特殊测井资料能够有效的检测破碎带层段以及估算裂隙密度。但是这种方法不仅成本高和耗费时间,而且不能满足野外实际生产的需要。
发明内容
本发明的目的是提供一种破碎带裂隙密度确定方法及系统,以实现节省成本和时间,满足油田、煤田以及金属矿等野外实际生产的需要。
为实现上述目的,本发明提供了一种破碎带裂隙密度确定方法,所述方法包括:
获取野外的测井曲线;
根据所述测井曲线确定敏感测井曲线;
利用测井曲线融合技术根据所述敏感测井曲线确定裂隙融合曲线;
采用离散小波变换方法根据所述裂隙融合曲线确定破碎带层段;
在破碎带层段,依据标准钻孔的成像测井资料和所述敏感测井曲线,建立裂隙密度关系式;
根据所述裂隙密度关系式确定待检测钻孔破碎带的裂隙密度。
可选的,所述根据所述测井曲线确定敏感测井曲线,具体包括:
对所述测井曲线进行预处理;
绘制破碎带和非破碎带预处理后的所述测井曲线的频率直方图;
根据所述频率直方图确定敏感测井曲线。
可选的,所述利用测井曲线融合技术根据所述敏感测井曲线确定裂隙融合曲线;具体包括:
从所述敏感测井曲线中选取一条敏感度最大的特征曲线作为基值曲线,并令j=1;
对所述敏感测井曲线中剩余的特征曲线进行归一化处理,获得归一化测井曲线;
对所述归一化测井曲线中的第j条特征曲线进行加权处理,获得第一加权曲线;
根据所述第一加权曲线和所述基值曲线确定第二加权曲线;
判断j是否小于所述归一化测井曲线中特征曲线的个数;如果j小于所述归一化测井曲线中特征曲线的个数,则将所述第二加权曲线作为基值曲线,并令j=j+1,返回步骤“对所述归一化测井曲线中的第j条特征曲线进行加权处理,获得第一加权曲线”;如果j大于或等于所述归一化测井曲线中特征曲线的个数,则将所述第二加权曲线作为裂隙融合曲线。
可选的,所述采用离散小波变换方法根据所述裂隙融合曲线确定破碎带层段,具体包括:
选择不同类型的小波函数对裂隙融合曲线进行离散小波变换,获得不同小波函数对应的不同分解尺度的高频小波系数;
根据所述不同分解尺度的高频小波系数确定最优分解尺度和最优小波函数;
利用所述最优小波函数在最优分解尺度下对所述裂隙融合曲线进行离散小波分解,获得最优高频小波系数;
利用所述最优高频小波系数检测破碎带层段。
可选的,所述预处理包括AD值转化、非点剔除、深度归位和环境校正。
本发明还提供一种破碎带裂隙密度确定系统,所述系统包括:
获取模块,用于获取野外的测井曲线;
敏感测井曲线确定模块,用于根据所述测井曲线确定敏感测井曲线;
裂隙融合曲线确定模块,用于利用测井曲线融合技术根据所述敏感测井曲线确定裂隙融合曲线;
破碎带层段确定模块,用于采用离散小波变换方法根据所述裂隙融合曲线确定破碎带层段;
裂隙密度关系式建立模块,用于在破碎带层段,依据标准钻孔的成像测井资料和所述敏感测井曲线,建立裂隙密度关系式;
裂隙密度确定模块,用于根据所述裂隙密度关系式确定待检测钻孔破碎带的裂隙密度。
可选的,所述敏感测井曲线确定模块,具体包括:
预处理单元,用于对所述测井曲线进行预处理;
频率直方图绘制单元,用于绘制破碎带和非破碎带预处理后的所述测井曲线的频率直方图;
敏感测井曲线确定单元,用于根据所述频率直方图确定敏感测井曲线。
可选的,所述裂隙融合曲线确定模块,具体包括:
基值曲线确定单元,用于从所述敏感测井曲线中选取一条敏感度最大的特征曲线作为基值曲线,并令j=1;
归一化处理单元,用于对所述敏感测井曲线中剩余的特征曲线进行归一化处理,获得归一化测井曲线;
加权处理单元,用于对所述归一化测井曲线中的第j条特征曲线进行加权处理,获得第一加权曲线;
第二加权曲线确定单元,用于根据所述第一加权曲线和所述基值曲线确定第二加权曲线;
判断单元,用于判断j是否小于所述归一化测井曲线中特征曲线的个数;如果j小于所述归一化测井曲线中特征曲线的个数,则将所述第二加权曲线作为基值曲线,并令j=j+1,返回步骤“对所述归一化测井曲线中的第j条特征曲线进行加权处理,获得第一加权曲线”;如果j大于或等于所述归一化测井曲线中特征曲线的个数,则将所述第二加权曲线作为裂隙融合曲线。
可选的,所述破碎带层段确定模块,具体包括:
高频小波系数确定单元,用于选择不同类型的小波函数对裂隙融合曲线进行离散小波变换,获得不同小波函数对应的不同分解尺度的高频小波系数;
最优参数确定单元,用于根据所述不同分解尺度的高频小波系数确定最优分解尺度和最优小波函数;
最优高频小波系数确定单元,用于利用所述最优小波函数在最优分解尺度下对所述裂隙融合曲线进行离散小波分解,获得最优高频小波系数;
破碎带层段检测单元,用于利用所述最优高频小波系数检测破碎带层段。
可选的,所述预处理包括AD值转化、非点剔除、深度归位和环境校正。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明提供一种破碎带裂隙密度确定方法及系统,所述方法综合利用测井曲线融合技术和离散小波变换方法确定破碎带层段,依据标准钻孔的成像测井资料,建立测井曲线与裂隙密度之间关系,以此确定待检测钻孔破碎带的裂隙密度,不仅能节省成本和时间,还能满足油田、煤田以及金属矿等野外实际生产的需要。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例破碎带裂隙密度确定方法流程图;
图2为本发明实施例破碎带裂隙密度确定系统结构图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种破碎带裂隙密度确定方法及系统,以实现节省成本和时间,满足油田、煤田以及金属矿等野外实际生产的需要。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明实施例破碎带裂隙密度确定方法流程图,如图1所示,本发明提供一种破碎带裂隙密度确定方法,所述方法包括:
步骤S1:获取野外测井曲线;
步骤S2:根据所述测井曲线确定敏感测井曲线;
步骤S3:利用测井曲线融合技术根据所述敏感测井曲线确定裂隙融合曲线;
步骤S4:采用离散小波变换方法根据所述裂隙融合曲线确定破碎带层段;
步骤S5:在破碎带层段,依据标准钻孔的成像测井资料和所述敏感测井曲线,建立裂隙密度关系式;
步骤S6:根据所述裂隙密度关系式确定待检测钻孔地层破碎带的裂隙密度。
下面对各个步骤进行详细论述:
步骤S1:获取野外测井曲线;所述测井曲线包括井径特征曲线、自然伽马特征曲线、三侧向电阻率特征曲线、声波时差特征曲线、补偿密度特征曲线、补偿中子特征曲线和岩性密度特征曲线。
步骤S2:根据所述测井曲线确定敏感测井曲线,具体包括:
步骤S21:对所述测井曲线进行预处理;所述预处理包括AD值转化、非点剔除、深度归位和环境校正;
步骤S22:绘制破碎带和非破碎带预处理后的所述测井曲线的频率直方图;
步骤S23:根据所述频率直方图确定对破碎带反映敏感的敏感测井曲线;所述敏感测井曲线包括井径特征曲线、自然伽马特征曲线、三侧向电阻率特征曲线、声波时差特征曲线、补偿密度特征曲线、补偿中子特征曲线和岩性密度特征曲线;所述敏感测井曲线在所述频率直方图中表现为单条测井曲线在破碎带和非破碎带具有不同的分布区间,且没有或有小于设定值的重叠区间。
步骤S3:所述利用测井曲线融合技术根据所述敏感测井曲线确定裂隙融合曲线;具体包括:
步骤S31:从所述敏感测井曲线中选取一条敏感度最大的特征曲线作为基值曲线,并令j=1;通常所述基值曲线为补偿密度曲线;所述敏感度最大的曲线为重叠区域最大的特征曲线;
步骤S32:对所述敏感测井曲线中剩余的特征曲线进行归一化处理,获得归一化测井曲线;
为了保证所述测井曲线中各条特征曲线贡献一致,对所述敏感测井曲线中剩余的特征曲线进行归一化处理,使各条特征曲线数值范围分布在[0,1]之间,具体归一化处理公式如下:
式中,Cij为第i个样品第j个特征曲线归一化后的特征曲线,cij为第i个样品第j个特征曲线归一化前的特征曲线,cmin为该特征曲线所有数据点中的最小值,cmax为该特征曲线所有数据点中的最大值,n为参与整合测井曲线中特征曲线的个数。
步骤S33:对所述归一化测井曲线中的第j条特征曲线进行加权处理,获得第一加权曲线;具体公式为:
其中,ωj为所述归一化测井曲线中的第j条特征曲线的加权系数,Cij为归一化后的第i个样品第j个特征曲线,为Cij的平均值,Hj为中间参数变量,m为样品数量,n为参与整合测井曲线中特征曲线的个数,cij为归一化前的第i个样品第j个特征曲线。
步骤S34:根据所述第一加权曲线和所述基值曲线确定第二加权曲线;具体公式为:
步骤S35:判断j是否小于所述归一化测井曲线中特征曲线的个数;如果j小于所述归一化测井曲线中特征曲线的个数,则将所述第二加权曲线作为基值曲线,并令j=j+1,返回步骤“对所述归一化测井曲线中的第j条特征曲线进行加权处理,获得第一加权曲线”;如果j大于或等于所述归一化测井曲线中特征曲线的个数,则将所述第二加权曲线作为裂隙融合曲线。
步骤S4:采用离散小波变换方法根据所述裂隙融合曲线确定破碎带层段,具体包括:
步骤S41:选择不同类型的小波函数对裂隙融合曲线进行离散小波变换,获得不同小波函数对应的不同分解尺度的高频小波系数;具体小波函数包括Haar小波、sym小波、coif小波、bior小波、rbio小波等。
步骤S42:根据所述不同分解尺度的高频小波系数确定最优分解尺度和最优小波函数;具体包括:
根据所述不同分解尺度的高频小波系数确定高频小波系数的波动率;
利用小波系数的波动率确定最优分解尺度和最优小波函数。
步骤S43:利用所述最优小波函数在最优分解尺度下对所述裂隙融合曲线进行离散小波分解,获得高频小波系数;具体离散小波分解过程如下:
其中,f(t)为裂隙融合曲线,t为深度变量,与测井曲线的测量深度有关,k为位移因子,φ()为尺度函数, h1(t-2k)为高通滤波函数,h0(t-2k)为低通滤波函数,为小波函数,ψ()为尺度函数的正交函数,为信号第(j-1)层分解的低频小波系数,信号第(j-1)层分解的高频小波系数;为信号第(j-2)分解的低频部分。
步骤S5:在破碎带层段,依据标准钻孔的成像测井资料和所述敏感测井曲线,建立裂隙密度关系式;具体公式为:
其中,log()为对数函数,FD为裂缝密度,A和B分别为区域特征常数,MEL为测井能量平均值,所述测井能量平均值为敏感测井曲线中多条特征曲线平均值的均值,ELp为第p个破碎带内单条特征曲线的平均值,c为每个钻孔内破碎带的总个数,np为第p个破碎带内单条特征曲线对应的总个数,LOGpq为第p个破碎带内第q条特征曲线的数值,d为一个破碎带内单条特征曲线对应的总个数。
图2为本发明实施例破碎带裂隙密度确定系统结构图,如图2所示,本发明还提供一种破碎带裂隙密度确定系统,所述系统包括:
获取模块1,用于获取野外的测井曲线;
敏感测井曲线确定模块2,用于根据所述测井曲线确定敏感测井曲线;
裂隙融合曲线确定模块3,用于利用测井曲线融合技术根据所述敏感测井曲线确定裂隙融合曲线;
破碎带层段确定模块4,用于采用离散小波变换方法根据所述裂隙融合曲线确定破碎带层段;
裂隙密度关系式建立模块5,用于在破碎带层段,依据标准钻孔的成像测井资料和所述敏感测井曲线,建立裂隙密度关系式;
裂隙密度确定模块6,用于根据所述裂隙密度关系式确定待检测钻孔破碎带的裂隙密度。
下面对各个模块进行详细论述:
所述敏感测井曲线确定模块2,具体包括:
预处理单元,用于对所述测井曲线进行预处理;所述预处理包括AD值转化、非点剔除、深度归位和环境校正;
频率直方图绘制单元,用于绘制破碎带和非破碎带预处理后的所述测井曲线的频率直方图;
敏感测井曲线确定单元,用于根据所述频率直方图确定敏感测井曲线。
所述裂隙融合曲线确定模块3,具体包括:
基值曲线确定单元,用于从所述敏感测井曲线中选取一条敏感度最大的特征曲线作为基值曲线,并令j=1;
归一化处理单元,用于对所述敏感测井曲线中剩余的特征曲线进行归一化处理,获得归一化测井曲线;
加权处理单元,用于对所述归一化测井曲线中的第j条特征曲线进行加权处理,获得第一加权曲线;
第二加权曲线确定单元,用于根据所述第一加权曲线和所述基值曲线确定第二加权曲线;
判断单元,用于判断j是否小于所述归一化测井曲线中特征曲线的个数;如果j小于所述归一化测井曲线中特征曲线的个数,则将所述第二加权曲线作为基值曲线,并令j=j+1,返回步骤“对所述归一化测井曲线中的第j条特征曲线进行加权处理,获得第一加权曲线”;如果j大于或等于所述归一化测井曲线中特征曲线的个数,则将所述第二加权曲线作为裂隙融合曲线。
所述破碎带层段确定模块4,具体包括:
高频小波系数确定单元,用于选择不同类型的小波函数对裂隙融合曲线进行离散小波变换,获得不同小波函数对应的不同分解尺度的高频小波系数;
最优参数确定单元,用于根据所述不同分解尺度的高频小波系数确定最优分解尺度和最优小波函数;
最优高频小波系数确定单元,用于利用所述最优小波函数在最优分解尺度下对所述裂隙融合曲线进行离散小波分解,获得最优高频小波系数;
破碎带层段检测单元,用于利用所述最优高频小波系数检测破碎带层段。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种破碎带裂隙密度确定方法,其特征在于,所述方法包括:
获取野外的测井曲线;
根据所述测井曲线确定敏感测井曲线;
利用测井曲线融合技术根据所述敏感测井曲线确定裂隙融合曲线;
采用离散小波变换方法根据所述裂隙融合曲线确定破碎带层段;
在破碎带层段,依据标准钻孔的成像测井资料和所述敏感测井曲线,建立裂隙密度关系式;
根据所述裂隙密度关系式确定待检测钻孔破碎带的裂隙密度;
所述利用测井曲线融合技术根据所述敏感测井曲线确定裂隙融合曲线,具体包括:
从所述敏感测井曲线中选取一条敏感度最大的特征曲线作为基值曲线,并令j=1;
对所述敏感测井曲线中剩余的特征曲线进行归一化处理,获得归一化测井曲线;
对所述归一化测井曲线中的第j条特征曲线进行加权处理,获得第一加权曲线;
根据所述第一加权曲线和所述基值曲线确定第二加权曲线;
判断j是否小于所述归一化测井曲线中特征曲线的个数;如果j小于所述归一化测井曲线中特征曲线的个数,则将所述第二加权曲线作为基值曲线,并令j=j+1,返回步骤“对所述归一化测井曲线中的第j条特征曲线进行加权处理,获得第一加权曲线”;如果j大于或等于所述归一化测井曲线中特征曲线的个数,则将所述第二加权曲线作为裂隙融合曲线。
2.根据权利要求1所述的破碎带裂隙密度确定方法,其特征在于,所述根据所述测井曲线确定敏感测井曲线,具体包括:
对所述测井曲线进行预处理;
绘制破碎带和非破碎带预处理后的所述测井曲线的频率直方图;
根据所述频率直方图确定敏感测井曲线。
3.根据权利要求1所述的破碎带裂隙密度确定方法,其特征在于,所述采用离散小波变换方法根据所述裂隙融合曲线确定破碎带层段,具体包括:
选择不同类型的小波函数对裂隙融合曲线进行离散小波变换,获得不同小波函数对应的不同分解尺度的高频小波系数;
根据所述不同分解尺度的高频小波系数确定最优分解尺度和最优小波函数;
利用所述最优小波函数在最优分解尺度下对所述裂隙融合曲线进行离散小波分解,获得最优高频小波系数;
利用所述最优高频小波系数检测破碎带层段。
4.根据权利要求2所述的破碎带裂隙密度确定方法,其特征在于,所述预处理包括AD值转化、非点剔除、深度归位和环境校正。
5.一种破碎带裂隙密度确定方法,其特征在于,所述方法包括:
获取野外的测井曲线;
根据所述测井曲线确定敏感测井曲线;
利用测井曲线融合技术根据所述敏感测井曲线确定裂隙融合曲线;
采用离散小波变换方法根据所述裂隙融合曲线确定破碎带层段;
在破碎带层段,依据标准钻孔的成像测井资料和所述敏感测井曲线,建立裂隙密度关系式;
根据所述裂隙密度关系式确定待检测钻孔破碎带的裂隙密度;
所述采用离散小波变换方法根据所述裂隙融合曲线确定破碎带层段,具体包括:
选择不同类型的小波函数对裂隙融合曲线进行离散小波变换,获得不同小波函数对应的不同分解尺度的高频小波系数;
根据所述不同分解尺度的高频小波系数确定最优分解尺度和最优小波函数;
利用所述最优小波函数在最优分解尺度下对所述裂隙融合曲线进行离散小波分解,获得最优高频小波系数;
利用所述最优高频小波系数检测破碎带层段。
6.一种破碎带裂隙密度确定系统,其特征在于,所述系统包括:
获取模块,用于获取野外的测井曲线;
敏感测井曲线确定模块,用于根据所述测井曲线确定敏感测井曲线;
裂隙融合曲线确定模块,用于利用测井曲线融合技术根据所述敏感测井曲线确定裂隙融合曲线;
破碎带层段确定模块,用于采用离散小波变换方法根据所述裂隙融合曲线确定破碎带层段;
裂隙密度关系式建立模块,用于在破碎带层段,依据标准钻孔的成像测井资料和所述敏感测井曲线,建立裂隙密度关系式;
裂隙密度确定模块,用于根据所述裂隙密度关系式确定待检测钻孔破碎带的裂隙密度;
所述裂隙融合曲线确定模块,具体包括:
基值曲线确定单元,用于从所述敏感测井曲线中选取一条敏感度最大的特征曲线作为基值曲线,并令j=1;
归一化处理单元,用于对所述敏感测井曲线中剩余的特征曲线进行归一化处理,获得归一化测井曲线;
加权处理单元,用于对所述归一化测井曲线中的第j条特征曲线进行加权处理,获得第一加权曲线;
第二加权曲线确定单元,用于根据所述第一加权曲线和所述基值曲线确定第二加权曲线;
判断单元,用于判断j是否小于所述归一化测井曲线中特征曲线的个数;如果j小于所述归一化测井曲线中特征曲线的个数,则将所述第二加权曲线作为基值曲线,并令j=j+1,返回步骤“对所述归一化测井曲线中的第j条特征曲线进行加权处理,获得第一加权曲线”;如果j大于或等于所述归一化测井曲线中特征曲线的个数,则将所述第二加权曲线作为裂隙融合曲线。
7.根据权利要求6所述的破碎带裂隙密度确定系统,其特征在于,所述敏感测井曲线确定模块,具体包括:
预处理单元,用于对所述测井曲线进行预处理;
频率直方图绘制单元,用于绘制破碎带和非破碎带预处理后的所述测井曲线的频率直方图;
敏感测井曲线确定单元,用于根据所述频率直方图确定敏感测井曲线。
8.根据权利要求6所述的破碎带裂隙密度确定系统,其特征在于,所述破碎带层段确定模块,具体包括:
高频小波系数确定单元,用于选择不同类型的小波函数对裂隙融合曲线进行离散小波变换,获得不同小波函数对应的不同分解尺度的高频小波系数;
最优参数确定单元,用于根据所述不同分解尺度的高频小波系数确定最优分解尺度和最优小波函数;
最优高频小波系数确定单元,用于利用所述最优小波函数在最优分解尺度下对所述裂隙融合曲线进行离散小波分解,获得最优高频小波系数;
破碎带层段检测单元,用于利用所述最优高频小波系数检测破碎带层段。
9.根据权利要求7所述的破碎带裂隙密度确定系统,其特征在于,所述预处理包括AD值转化、非点剔除、深度归位和环境校正。
10.一种破碎带裂隙密度确定系统,其特征在于,所述系统包括:
获取模块,用于获取野外的测井曲线;
敏感测井曲线确定模块,用于根据所述测井曲线确定敏感测井曲线;
裂隙融合曲线确定模块,用于利用测井曲线融合技术根据所述敏感测井曲线确定裂隙融合曲线;
破碎带层段确定模块,用于采用离散小波变换方法根据所述裂隙融合曲线确定破碎带层段;
裂隙密度关系式建立模块,用于在破碎带层段,依据标准钻孔的成像测井资料和所述敏感测井曲线,建立裂隙密度关系式;
裂隙密度确定模块,用于根据所述裂隙密度关系式确定待检测钻孔破碎带的裂隙密度;
所述破碎带层段确定模块,具体包括:
高频小波系数确定单元,用于选择不同类型的小波函数对裂隙融合曲线进行离散小波变换,获得不同小波函数对应的不同分解尺度的高频小波系数;
最优参数确定单元,用于根据所述不同分解尺度的高频小波系数确定最优分解尺度和最优小波函数;
最优高频小波系数确定单元,用于利用所述最优小波函数在最优分解尺度下对所述裂隙融合曲线进行离散小波分解,获得最优高频小波系数;
破碎带层段检测单元,用于利用所述最优高频小波系数检测破碎带层段。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910948240.7A CN110656933B (zh) | 2019-10-08 | 2019-10-08 | 一种破碎带裂隙密度确定方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910948240.7A CN110656933B (zh) | 2019-10-08 | 2019-10-08 | 一种破碎带裂隙密度确定方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110656933A CN110656933A (zh) | 2020-01-07 |
CN110656933B true CN110656933B (zh) | 2022-12-13 |
Family
ID=69038503
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910948240.7A Active CN110656933B (zh) | 2019-10-08 | 2019-10-08 | 一种破碎带裂隙密度确定方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110656933B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103293563A (zh) * | 2013-05-13 | 2013-09-11 | 中国石油大学(华东) | 一种确定油气储层岩石裂隙发育度和流体性质的方法 |
CN103412336A (zh) * | 2013-07-22 | 2013-11-27 | 中国石油天然气股份有限公司 | 一种非均质油藏中岩石系统的纵波速度预测方法 |
CN104155701A (zh) * | 2014-07-03 | 2014-11-19 | 中国石油天然气集团公司 | 一种利用叠前地震资料及井信息的多尺度裂缝预测方法 |
WO2016041189A1 (zh) * | 2014-09-19 | 2016-03-24 | 杨顺伟 | 一种评价页岩气储层及寻找甜点区的方法 |
CN109407173A (zh) * | 2018-09-29 | 2019-03-01 | 核工业北京地质研究院 | 基于常规测井曲线的岩性精细化和自动化识别方法 |
CN109471168A (zh) * | 2018-11-06 | 2019-03-15 | 河海大学 | 一种孔裂隙介质中纵波速度与衰减的预测方法 |
-
2019
- 2019-10-08 CN CN201910948240.7A patent/CN110656933B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103293563A (zh) * | 2013-05-13 | 2013-09-11 | 中国石油大学(华东) | 一种确定油气储层岩石裂隙发育度和流体性质的方法 |
CN103412336A (zh) * | 2013-07-22 | 2013-11-27 | 中国石油天然气股份有限公司 | 一种非均质油藏中岩石系统的纵波速度预测方法 |
CN104155701A (zh) * | 2014-07-03 | 2014-11-19 | 中国石油天然气集团公司 | 一种利用叠前地震资料及井信息的多尺度裂缝预测方法 |
WO2016041189A1 (zh) * | 2014-09-19 | 2016-03-24 | 杨顺伟 | 一种评价页岩气储层及寻找甜点区的方法 |
CN109407173A (zh) * | 2018-09-29 | 2019-03-01 | 核工业北京地质研究院 | 基于常规测井曲线的岩性精细化和自动化识别方法 |
CN109471168A (zh) * | 2018-11-06 | 2019-03-15 | 河海大学 | 一种孔裂隙介质中纵波速度与衰减的预测方法 |
Non-Patent Citations (2)
Title |
---|
孔、裂隙并存地层中的声波测井理论及多极子声场特征;陈雪莲等;《地球物理学报》;20120630;第55卷(第6期);第2129-2140页 * |
辽河油田兴隆台太古界古潜山储层评价标准研究;王博等;《高校地质学报》;20100930;第16卷(第3期);第365-363页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110656933A (zh) | 2020-01-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11598203B2 (en) | Reducing resonant noise in seismic data acquired using a distributed acoustic sensing system | |
US10180061B2 (en) | Methods of evaluating rock properties while drilling using downhole acoustic sensors and a downhole broadband transmitting system | |
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 | |
CN105116442A (zh) | 岩性油气藏弱反射地震信号的重构方法 | |
US20120044784A1 (en) | Determining a position of a geological layer relative to a wavelet response in seismic data | |
EP2761337B1 (en) | Methods of evaluating rock properties while drilling using downhole acoustic sensors and a downhole broadband transmitting system | |
CN108375785B (zh) | 裂缝带位置校正方法及装置 | |
WO2020106287A1 (en) | Enhanced anisotropy analysis with multi-component dipole sonic data | |
GB2472300A (en) | Using dipole compressional data obtained in a borehole to determined properties of a subterranean formation | |
CN106154328B (zh) | 一种基于互相关的微地震事件识别拾取方法及其系统 | |
CN105277979B (zh) | 一种地震属性的优化方法和装置 | |
CN111046328A (zh) | 基于测井曲线小波Mallet算法的岩石类型识别方法 | |
WO2017127045A1 (en) | Method of minimizing tool response for downhole logging operations | |
CN117031566A (zh) | 一种地磁异常数据集的构建方法及系统 | |
CN110656933B (zh) | 一种破碎带裂隙密度确定方法及系统 | |
CN112112618A (zh) | 一种基于流程化的压裂选井选层方法及系统 | |
CN113327070B (zh) | 智能勘查煤系气的方法、装置及电子设备 | |
US20220229202A1 (en) | Method and system for automatic picking of borehole acoustic events based on new objective function | |
CN115879020A (zh) | 储集体的分类刻画方法、装置、电子设备及存储介质 | |
CN116027431A (zh) | 储层含气性预测方法、装置、存储介质及电子设备 | |
Muhamedyev et al. | Geophysical research of boreholes: Artificial neural networks data analysis | |
Paap et al. | A neural network approach for improved seismic event detection in the Groningen gas field, the Netherlands | |
CN112558147B (zh) | 一种井中微地震数据的偏振分析方法及系统 | |
CN117130066A (zh) | 区域裂缝识别方法、装置、存储介质及电子设备 | |
CN114076990B (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 |