CN109372497B - 一种超声成像动态均衡处理的方法 - Google Patents
一种超声成像动态均衡处理的方法 Download PDFInfo
- Publication number
- CN109372497B CN109372497B CN201810948466.2A CN201810948466A CN109372497B CN 109372497 B CN109372497 B CN 109372497B CN 201810948466 A CN201810948466 A CN 201810948466A CN 109372497 B CN109372497 B CN 109372497B
- Authority
- CN
- China
- Prior art keywords
- image
- result
- frequency component
- ultrasonic
- highcom
- 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
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
- 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
Abstract
本发明公开了一种超声成像动态均衡处理的方法,包括:对输入超声成像图像进行平滑滤波,将超声输入图像与进行平滑滤波后的图像进行相减,得到高频分量,对平滑滤波图像进行椭圆拟合出椭圆中心偏移量,然后对该椭圆中心偏移量进行偏心校正,将校正结果与所得的高频分量相加,得到输入超声成像图像预处理结果;对预处理结果进行动态均衡处理,得到最终的处理结果。本发明处理效果,明显优于国外贝克休斯软件处理效果,能够显著提高裂缝等地质特征的分辨率和对比度。
Description
技术领域
本发明属于石油开发和工程领域,是一种超声成像动态均衡处理的方法。
背景技术
随着测井技术的发展,超声成像仪器在裸眼井、套管井都有广泛的应用。在裸眼井中和微电阻率扫描组合使用,在井旁构造缝洞识别等有着重要的应用。套管井中在套损检测也有着一定的应用。超声仪器测量的数据值动态范围大,分辨率低,需要对成像数据进行进一步的增强处理,提高图像的分辨率、对比度。图像动态均衡处理就是一个很重要的处理方法。但是现有的传统动态均衡处理方法,存在一些问题,只从图像均衡处理角度进行处理,没有考虑偏心的影响,处理导致处理结果的裂缝显示不清晰。本发明采用了新的流程进行处理,达到了较好的处理效果。
发明内容
本发明的目的是针对目前超声图像动态均衡处理存在的分辨率不高、裂缝显示不清晰等问题,采用了新的均衡处理流程。本发明处理流程具备创新性,处理效果优于目前在用的国外贝克休斯软件处理效果,已经应用到了油田超声图像的数据处理中,得到了较好的处理效果,得到了油田的肯定和认可。
本发明的目的是通过下述技术方案来实现的。
一种超声成像动态均衡处理方法,包括以下步骤:
步骤1,对输入超声成像图像进行平滑滤波;
步骤2,将超声输入图像与进行平滑滤波后的图像进行相减,得到高频分量;
步骤3,对平滑滤波图像进行椭圆拟合出椭圆中心偏移量,然后对该椭圆中心偏移量进行偏心校正;
步骤4,将步骤3的校正结果与步骤2所得的高频分量相加,得到输入超声成像图像预处理结果;
步骤5,对预处理结果进行动态均衡处理,得到最终的处理结果。
对于上述技术方案,本发明还有进一步优选的方案:
优选的,对超声成像图像进行11*11~25*25滤波,滤除图像细节分量,保留图像低频分量,得到平滑滤波图像。
优选的,高频分量通过以下方式得到:
超声图像IMG是由高频细节分量HighCom和低频分量LowCom组成,计算高频分量为后续处理做准备。由下式计算高频分量:
HighCom(i,j)=IMG(i,j)-LowCom(i,j)
i=0,1,…N-1,j=0,1,…M-1,
其中,IMG表示超声图像数据,HighCom表示超声图像高频细节分量,LowCom表示超声图像低频细节分量,N,M表示超声图像IMG、高频细节分量HighCom、低频分量LowCom共有N行,M列数据,i,j为整数,表示高频细节分量HighCom和低频分量LowCom的第i行,第j列。
优选的,对平滑滤波图像进行椭圆拟合偏心校正,通过下式进行校正:
θj=(2*π*j/M),
i=0,1,…N-1,j=0,1,…M-1,
其中,(x0,y0)表示椭圆中心偏移量,θj表示低频分量LowCom的(i,j)数据对应的方位。
其中,Corr(i,j)表示超声图像偏心校正结果,(x0,y0)表示椭圆中心偏移量,θj表示低频分量LowCom的(i,j)数据对应的方位,i,j为整数,表示超声图像偏心校正结果Corr和低频分量LowCom的第i行,第j列,N,M表示超声图像偏心校正结果Corr、低频分量LowCom共有N行,M列数据。
优选的,步骤4中,输入超声成像图像校正结果通过下式得到:
PreResult(i,j)=Corr(i,j)+HighCom(i,j)
i=0,1,…N-1,j=0,1,…M-1。
其中,PreResult(i,j)表示超声图像预处理结果,Corr(i,j)表示超声图像偏心校正结果,HighCom(i,j)表示高频细节分量,i,j为整数,表示偏心校正结果Corr(i,j)和高频细节分量HighCom的第i行,第j列,N,M表示超声图像预处理结果PreResult、高频分量HighCom共有N行,M列数据。
相对于现有传统图像均衡处理技术,本发明的有益效果在于:
本发明首先对超声图像进行17*17窗长滤波,原图像与滤波结果相减,得到高频细节分量,对滤波结果进行偏心校正,偏心校正结果和高频细节分量相加,得到新的预处理结果,最后对新的校正结果,进行动态均衡处理。本发明处理效果,明显优于国外贝克休斯软件处理效果,能够显著提高裂缝等地质特征的分辨率和对比度,得到了高度认可。
附图说明
图1为本发明处理流程示意图;
图2a-2f分别为实施例1中处理效果Ⅰ示意图;其中,图2a为输入超声成像图像;图2b为原静态图像处理效果示意图;图2c为原动态图像处理效果示意图;图2d为采用本发明方法处理图像的处理效果示意图;图2e为国外贝壳修斯静态图像处理效果示意图;图2f为国外贝壳修斯动态图像处理效果示意图;
图3a-3f为实施例1中处理效果Ⅱ示意图;其中,图3a为输入超声成像图像;图3b为原静态图像处理效果示意图;图3c为原动态图像处理效果示意图;图3d为采用本发明方法处理图像的处理效果示意图;图3e为国外贝壳修斯静态图像处理效果示意图;图3f为国外贝壳修斯动态图像处理效果示意图。
具体实施方式
下面结合具体相关流量测井数据,对本发明的具体的实施方式进行说明。
本发明的超声成像动态均衡处理方法,包括以下步骤:
步骤1,对输入超声成像图像进行平滑滤波,程序中采用17*17滤波。对超声成像图像进行11*11~25*25滤波,滤除图像细节分量,保留图像低频分量,得到平滑滤波图像。
步骤2,将超声输入图像与进行平滑滤波后的图像进行相减,得到高频分量。
超声图像IMG是由高频细节分量HighCom和低频分量LowCom组成,由下式计算高频分量:
HighCom(i,j)=IMG(i,j)-LowCom(i,j)
i=0,1,…N-1,j=0,1,…M-1,
其中,IMG表示超声图像数据,HighCom表示超声图像高频细节分量,LowCom表示超声图像低频细节分量,N,M表示超声图像IMG共有N行,M列数据,i,j为整数,表示高频细节分量HighCom和低频分量LowCom的第i行,第j列。
步骤3,对平滑滤波图像进行椭圆拟合出椭圆中心偏移量,然后对该椭圆中心偏移量进行偏心校正,通过下式进行校正:
θj=(2*π*j/M),
i=0,1,…N-1,j=0,1,…M-1,
其中,Corr(i,j)表示超声图像偏心校正结果,(x0,y0)表示椭圆中心偏移量,θj表示低频分量LowCom的(i,j)数据对应的方位,i,j为整数,表示超声图像偏心校正结果Corr和低频分量LowCom的第i行,第j列。
步骤4,将步骤3的校正结果与步骤2所得的高频分量相加,得到输入超声成像图像预处理结果,通过下式得到:
PreResult(i,j)=Corr(i,j)+HighCom(i,j)
i=0,1,…N-1,j=0,1,…M-1
其中,PreResult(i,j)表示超声图像预处理结果,Corr(i,j)表示超声图像偏心校正结果,HighCom(i,j)表示高频细节分量,i,j为整数,表示偏心校正结果Corr(i,j)和高频细节分量HighCom的第i行,第j列,N,M表示超声图像共有N行,M列数据。
步骤5,对预处理结果进行动态均衡处理,得到最终的处理结果。
本程序处理流程参照处理流程图1。处理效果图2a-2f、图3a-3f。从处理效果图2a-2f来看,图2a为原图,原开发的静态图2b、动态图2c处理效果和贝克休斯处理的静态图2e、动态图2f相当,本发明处理效果图2d比贝壳分辨率高、裂缝清晰。图3a-3f中,6758-6760深度段贝克休斯处理的静态图3b、动态图3c像几乎没有裂缝信息,本发明处理效果图3d可以清晰的看到相关裂缝信息,处理效果很好。
以上对本发明实施例所提供的实施方案进行了详细介绍,本文中应用了具体个例对本发明实施例的原理以及实施方式进行了阐述,以上实施例的说明只适用于帮助理解本发明实施例的原理;同时,对于本领域的一般技术人员,依据本发明实施例,在具体实施方式以及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (4)
1.一种超声成像动态均衡处理的方法,其特征在于,包括以下步骤:
步骤1,对输入超声成像图像进行平滑滤波;
步骤2,将超声输入图像与进行平滑滤波后的图像进行相减,得到高频分量;
步骤3,对平滑滤波图像进行椭圆拟合出椭圆中心偏移量,然后对该椭圆中心偏移量进行偏心校正;
对平滑滤波图像进行椭圆拟合偏心校正,通过下式进行校正:
θj=(2*π*j/M),
i=0,1,…N-1,j=0,1,…M-1,
其中,Corr(i,j)表示超声图像偏心校正结果,(x0,y0)表示椭圆中心偏移量,θj表示低频分量LowCom的(i,j)数据对应的方位,i,j为整数,表示超声图像偏心校正结果Corr和低频分量LowCom的第i行,第j列;
步骤4,将步骤3的校正结果与步骤2所得的高频分量相加,得到输入超声成像图像预处理结果;
步骤5,对预处理结果进行动态均衡处理,得到最终的处理结果。
2.根据权利要求1所述的方法,其特征在于,步骤1中,对超声成像图像进行11*11~25*25滤波,滤除图像细节分量,保留图像低频分量,得到平滑滤波图像。
3.根据权利要求1所述的方法,其特征在于,步骤2中,高频分量通过以下方式得到:
超声图像IMG是由高频分量HighCom和低频分量LowCom组成,由下式计算高频分量:
HighCom(i,j)=IMG(i,j)-LowCom(i,j)
i=0,1,…N-1,j=0,1,…M-1,
其中,IMG表示超声图像数据,HighCom表示超声图像高频分量,LowCom表示超声图像低频分量,N,M表示超声图像IMG共有N行,M列数据,i,j为整数,表示高频分量HighCom和低频分量LowCom的第i行、第j列。
4.根据权利要求1所述的方法,其特征在于,步骤4中,输入超声成像图像预处理结果通过下式得到:
PreResult(i,j)=Corr(i,j)+HighCom(i,j)
i=0,1,…N-1,j=0,1,…M-1
其中,PreResult(i,j)表示超声图像预处理结果,Corr(i,j)表示超声图像偏心校正结果,HighCom(i,j)表示高频细节分量,i,j为整数,表示偏心校正结果Corr(i,j)和高频细节分量HighCom的第i行,第j列,N,M表示超声图像共有N行,M列数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810948466.2A CN109372497B (zh) | 2018-08-20 | 2018-08-20 | 一种超声成像动态均衡处理的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810948466.2A CN109372497B (zh) | 2018-08-20 | 2018-08-20 | 一种超声成像动态均衡处理的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109372497A CN109372497A (zh) | 2019-02-22 |
CN109372497B true CN109372497B (zh) | 2022-03-29 |
Family
ID=65404062
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810948466.2A Active CN109372497B (zh) | 2018-08-20 | 2018-08-20 | 一种超声成像动态均衡处理的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109372497B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2614706A1 (en) * | 2005-07-22 | 2007-02-08 | Halliburton Energy Services, Inc. | Ultrasonic imaging in wells or tubulars |
CN102914797A (zh) * | 2012-10-16 | 2013-02-06 | 中国石油天然气股份有限公司 | 一种获得地层各向异性系数的方法及装置 |
WO2013017105A1 (en) * | 2011-08-04 | 2013-02-07 | Shenzhen Institutes Of Advanced Technology Chinese Academy Of Sciences | Ultrasonic imaging system and method for measuring elasticity of biological tissues |
CN105422089A (zh) * | 2015-11-30 | 2016-03-23 | 长江大学 | 一种利用井周超声图像生成孔隙度谱的方法 |
CN105604541A (zh) * | 2015-12-28 | 2016-05-25 | 中国石油天然气集团公司 | 一种生产测井多臂井径斜井校正处理的方法 |
CN105825487A (zh) * | 2016-04-06 | 2016-08-03 | 中国海洋石油总公司 | 一种全井周电成像图像生成方法和系统 |
CN106443802A (zh) * | 2016-09-20 | 2017-02-22 | 中石化石油工程技术服务有限公司 | 一种基于形态学滤波的电成像储层缝洞体的定量表征方法及系统 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1284975C (zh) * | 2003-01-16 | 2006-11-15 | 上海交通大学 | 遥感图像双线性插值小波变换融合的优化方法 |
CN101984346A (zh) * | 2010-10-19 | 2011-03-09 | 浙江大学 | 基于低通滤波的水果表面缺陷检测方法 |
JP5791349B2 (ja) * | 2011-04-21 | 2015-10-07 | キヤノン株式会社 | 撮像装置及びその制御方法 |
JP2013066157A (ja) * | 2011-08-31 | 2013-04-11 | Sony Corp | 画像処理装置、および画像処理方法、並びにプログラム |
WO2016194161A1 (ja) * | 2015-06-03 | 2016-12-08 | 株式会社日立製作所 | 超音波診断装置、及び画像処理方法 |
US10222501B2 (en) * | 2016-01-25 | 2019-03-05 | Baker Hughes, A Ge Company, Llc | Televiewer image wood-grain reduction techniques |
-
2018
- 2018-08-20 CN CN201810948466.2A patent/CN109372497B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2614706A1 (en) * | 2005-07-22 | 2007-02-08 | Halliburton Energy Services, Inc. | Ultrasonic imaging in wells or tubulars |
WO2013017105A1 (en) * | 2011-08-04 | 2013-02-07 | Shenzhen Institutes Of Advanced Technology Chinese Academy Of Sciences | Ultrasonic imaging system and method for measuring elasticity of biological tissues |
CN102914797A (zh) * | 2012-10-16 | 2013-02-06 | 中国石油天然气股份有限公司 | 一种获得地层各向异性系数的方法及装置 |
CN105422089A (zh) * | 2015-11-30 | 2016-03-23 | 长江大学 | 一种利用井周超声图像生成孔隙度谱的方法 |
CN105604541A (zh) * | 2015-12-28 | 2016-05-25 | 中国石油天然气集团公司 | 一种生产测井多臂井径斜井校正处理的方法 |
CN105825487A (zh) * | 2016-04-06 | 2016-08-03 | 中国海洋石油总公司 | 一种全井周电成像图像生成方法和系统 |
CN106443802A (zh) * | 2016-09-20 | 2017-02-22 | 中石化石油工程技术服务有限公司 | 一种基于形态学滤波的电成像储层缝洞体的定量表征方法及系统 |
Non-Patent Citations (9)
Title |
---|
"Differentiating Open Natural Fractures From Healed Fractures Using The New, High-Definition Oil-Based Mud Microelectrical Imager – Case Studies From Organic Rich Shales";Bingjian Li;《society of petroleum engineers》;20150930;全文 * |
"图像增强算法实现";王蓉;《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》;20150115(第1期);第9-24页 * |
"基于超声测井图像的动态直方图均衡算法研究";涂继辉;《电视技术》;20111031;第35卷(第19期);第113-114,122页 * |
"基于限定对比度直方图均衡的超声测井图像增强方法";付青青,等;《测井技术》;20170630;第41卷(第3期);全文 * |
"套管井中仪器偏心对井壁超声电视图象的影响及其校正";余厚全,等;《应用声学》;19990415(第2期);全文 * |
"成像测井中超声图像退化原因分析与复原方法研究";张健;《计算技术与自动化》;20131231;第32卷(第4期);全文 * |
"改进的超声电视幅度图像偏心校正算法";涂继辉,等;《测井技术》;20111231;第35卷(第6期);全文 * |
"超声成像测井井径的椭圆拟合及偏心时间图像的校正";余厚全,等;《测井技术》;20111031;第35卷(第5期);第418-421页 * |
J.. Sugiura, R.. Lee."Removal of Cyclic Borehole Noise From Low- and High-Resolution LWD Images and Its Impact on Image Interpretation".《society of petroleum engineers》.2012,全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN109372497A (zh) | 2019-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107561738B (zh) | 基于fcn的tft-lcd表面缺陷快速检测方法 | |
WO2020192704A1 (zh) | 图像处理模型的训练方法、图像处理方法、装置及存储介质 | |
CN110223254A (zh) | 一种基于对抗生成网络的图像去噪方法 | |
CN104821011A (zh) | 基于相机拍摄的2d户型照片生成3d户型模型的方法 | |
CN111062978B (zh) | 基于频域滤波技术的时空图像测流的纹理识别方法 | |
Aleem et al. | Random data augmentation based enhancement: a generalized enhancement approach for medical datasets | |
CN109372497B (zh) | 一种超声成像动态均衡处理的方法 | |
CN113592923A (zh) | 一种基于深度局部特征匹配的批图像配准方法 | |
CN110080754B (zh) | 一种电成像测井图像类周期性干扰处理方法 | |
Li et al. | Multipatch unbiased distance non-local adaptive means with wavelet shrinkage | |
CN112017130B (zh) | 基于自适应各向异性全变分正则化的图像复原方法 | |
CN111724326B (zh) | 图像处理方法及装置、电子设备及存储介质 | |
CN111079689B (zh) | 一种指纹图像增强方法 | |
CN110264404B (zh) | 一种超分辨图像纹理优化的方法和装置 | |
CN107610218B (zh) | 一种面向立体结构网点三维图像重建的层面数据获取方法 | |
CN107464220B (zh) | 一种基于重力叠加模型的公路面层病害图像增强方法 | |
CN108364285A (zh) | 一种基于迭代相位一致性的裂缝提取方法和应用 | |
CN115188039A (zh) | 一种基于图像频域信息的深度伪造视频技术溯源方法 | |
CN104346786B (zh) | 基于Demons算法的图像去噪算法 | |
CN111369452A (zh) | 大区域图像局部破损点优化提取方法 | |
Wang et al. | Accurate extraction method for structural features of building facades through texture fusion | |
CN110910442A (zh) | 一种基于kernel-free图像复原的高速移动对象机器视觉尺寸检测方法 | |
CN108198139A (zh) | 基于传导全变分正则化的图像去噪方法 | |
CN113392782B (zh) | 面向弱小目标检测的遥感影像深度学习训练样本增强方法 | |
CN113223021B (zh) | 一种基于UNet的肺部X光图像分割方法 |
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 |