CN114779346A - 一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法及装置 - Google Patents
一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法及装置 Download PDFInfo
- Publication number
- CN114779346A CN114779346A CN202210355648.5A CN202210355648A CN114779346A CN 114779346 A CN114779346 A CN 114779346A CN 202210355648 A CN202210355648 A CN 202210355648A CN 114779346 A CN114779346 A CN 114779346A
- Authority
- CN
- China
- Prior art keywords
- component
- well
- dipole
- phi
- data
- 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.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 81
- 238000005259 measurement Methods 0.000 title claims abstract description 30
- 238000003384 imaging method Methods 0.000 claims abstract description 73
- 238000001514 detection method Methods 0.000 claims abstract description 18
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 43
- 238000012545 processing Methods 0.000 claims description 30
- 238000005755 formation reaction Methods 0.000 claims description 27
- 238000001914 filtration Methods 0.000 claims description 26
- 238000004088 simulation Methods 0.000 claims description 16
- 238000003786 synthesis reaction Methods 0.000 claims description 16
- 238000006243 chemical reaction Methods 0.000 claims description 15
- 230000005855 radiation Effects 0.000 claims description 13
- 230000000694 effects Effects 0.000 claims description 7
- 239000007788 liquid Substances 0.000 claims description 7
- 238000001308 synthesis method Methods 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 238000000926 separation method Methods 0.000 claims description 6
- 230000003993 interaction Effects 0.000 claims description 5
- 230000009467 reduction Effects 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 4
- 230000005284 excitation Effects 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000002146 bilateral effect Effects 0.000 claims description 2
- 238000004891 communication Methods 0.000 claims description 2
- 238000009413 insulation Methods 0.000 claims description 2
- 238000006073 displacement reaction Methods 0.000 description 9
- 239000012530 fluid Substances 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 3
- 150000001875 compounds Chemical class 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 230000005405 multipole Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 235000009508 confectionery Nutrition 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 239000000178 monomer Substances 0.000 description 1
- 230000005404 monopole Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000002847 sound insulator Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000002945 steepest descent method Methods 0.000 description 1
- 238000010189 synthetic method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Environmental & Geological Engineering (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Fluid Mechanics (AREA)
- Acoustics & Sound (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法及装置,主要包括:利用正交偶极声波测井装置在井中采集四分量数据,从中分离出井外反射信号并对其偏移成像;用四分量反射信号组成φ向发射和接收的合成信号XX,变化φ得到使XX达到最大幅值时的角度φ0;用四分量反射数据合成信号Ns1和Ns2并比较其相位差异;建立Ns1和Ns2之间相位差异随反射体方位变化的规律图版;判断反射体的真实方位并转换到固定坐标系下;对最大幅值信号XX(φ0)进行偏移成像,确定反射体的离井距离。本发明同时利用反射波的幅值和相位信息,有效消除了偶极横波远探测中的方位180°不确定性。此外,本发明具有测井效率高、成像结果准确、稳定性好和适用性强的优点。
Description
技术领域
本发明属于应用地球物理声学测井领域,具体涉及一种利用井中正交偶极声源发射-考虑不同数据采集方式的四分量阵列接收的声学测量模式,以及一种基于四分量反射信号的幅值和相位特征精准探测井外地质构造位置的方法。
背景技术
对于油气资源的勘探和开发来说,井外地质构造及地质体的有效探测对于寻找储层的油气“甜点”至关重要,在这一过程中偶极横波远探测测井技术发挥了巨大的作用。偶极横波远探测成像方法具有方位灵敏度强、径向探测深度深和反射信号易提取等优势,目前,该技术已经在裂缝、断层、溶蚀孔洞等地质构造的探测中取得了良好的应用效果,具有广阔的应用前景(唐晓明,古希浩,李杨虎,苏远大,2021,井孔与弹性波的相互作用:理论、方法和应用[J]:地球物理学报,64(12),4227-4238)。
然而,现有的偶极横波远探测技术存在反射体方位的180°不确定难题,其不利于指导定向侧钻,难以满足现场应用的需求。近年来,针对井外界面的方位模糊性问题,国内外学者已经开展了一系列研究,主要提出了以下三种解决方法:
第一种是采用偶极发射-单极接收的单偶极混合测量模式来消除方位不确定性(Tang,X.M.,2004,Imaging near-borehole structure using directional acoustic-wave measurement[J]:Geophysics,69,no.6,1378-1386.)。在该模式的基础上,Xu等人(Xu,J.Q.,H.S.Hu,and Z.Wang,2019,Asymptotic solution to a 3D dipole single-well imaging system with combined monopole and dipole receivers with anapplication in elimination of azimuth ambiguity[J]:Geophysics,84,no.5,D191-D207.)提出通过比较压力和位移与参考波形的极性来唯一确定反射体的方位角,但在实际测井中,与反射体性质相关的参考波形是未知的;Li等人(Li C.,H.Chen,X.He,andX.M.Wang,2021,Identifying reflector azimuth from borehole multicomponentcross-dipole acoustic measurement[J]:Geophysics,86,no.6,P.D201-D214.)利用偶极声源辐射的P波来消除SH波确定方位时造成的不确定性,然而偶极P波反射只有在软地层情况下较为明显。
第二种是利用三分量位移传感器接收完整的反射波场(Gong,H.,H.Chen,X.He,and X.M.Wang,2015,Eliminating the azimuth ambiguity in single-well imagingusing 3C sonic data[J]:Geophysics,80,no.1,A13-A17.),采用逆时偏移技术联合反演波场的水平分量和轴向分量,增强真实成像的同时降低伪像,但这种成像结果取决于反射波三分量的接收方式,并且该方法需要新的仪器支撑,目前该方法的实际应用尚未见报道。
第三种是通过偶极发射-八方位接收的测量模式来消除方位不确定性(Li,D.,W.X.Qiao,X.H.Che,X.D.Ju,S.B.Yang,J.Q.Lu,and B.Y.Men,2020,Eliminating theazimuth ambiguity in reflected S-wave imaging logging based on the azimuthalreceiver mode[J]:Journal of Petroleum Science and Engineering,199,108295.),但这种方法对接收换能器的性能要求较高,接收数据中反射波的信噪比一般较差。
目前,从经济性和实用性上看,偶极横波远探测存在的方位180°不确定问题还未得到妥善解决,因此亟需提出一种能够有效消除方位模糊性的新方法。
发明内容
本发明的目的是提供一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法及装置。以现有的正交偶极四分量测井原理为基础,采用不同的四分量数据采集模式,根据反射波的幅值和相位特征精准探测井外地质构造,可以经济有效地解决偶极横波远探测中存在的方位180°不确定的技术问题。
为实现上述发明目的,本发明提供了一种井中正交偶极四分量测量的声学测井装置,包括:正交偶极声源发射器短节、发射和接收之间的隔声短节、四个周向等间隔布置的方位单元接收站阵列短节和处理器短节。
所述正交偶极声源发射短节包括两个正交布置的偶极声源,声源指向相互垂直;所述接收站阵列短节包括至少八个轴向上等间隔布置的方位接收站;所述方位接收站包括四个周向上等间隔布置的接收换能器;所述处理器短节包括仪器方位测量装置以及与所述声源发射短节、接收站阵列短节和地面控制系统进行通信连接的模块。
此外,为实现上述发明目的,本发明同时提出了一种根据所述声学测量装置采集数据中反射波的幅值和相位特征精准探测井外地质构造的方法,其包括以下处理步骤:
步骤一、将所述正交偶极四分量声波测井装置放于测量井中,依次激发两个正交的偶极声源,向井外地层中辐射弹性波。
步骤二、记录仪器坐标轴(如y轴)相对于固定方向(如地球北极)的方位角。
步骤三、用方位阵列接收器采集信号并按照特定模式组成四分量数据(xx,xy,yx,yy)。
步骤四、将采集的仪器方位数据和四分量数据通过电缆上传至地面控制系统。
步骤五、对四分量数据进行滤波并做波场分离处理,压制和剔除沿井孔传播的直达波信号,提取出来自井外界面的反射信号。
步骤六、将从全波中分离出的四分量井外反射信号进行偏移成像。
步骤七、在仪器坐标系(x,y)下,以相对于仪器x轴的角度φ为变量,用四分量井外反射信号组成φ向发射和接收的合成XX反射波信号XX(φ)=xx·sin2φ-(xy+yx)·sinφcosφ+yy·cos2φ,变化φ值得到使XX达到最大幅值时的角度φ0。
步骤八、利用四分量反射数据合成两个新信号Ns1和Ns2,并比较它们之间的相位差异。
步骤九、按照步骤三中的四分量采集模式和步骤八中的新信号合成方式,通过理论模拟,建立新信号Ns1和Ns2之间相位差异随井外反射体方位变化的规律图版。
步骤十、根据理论模拟图版,判断φ0是否为反射体的真实方位角:若φ0与步骤八中信号Ns1和Ns2之间相位差异的关系符合步骤九中的理论图版规律,则φ0为真实方位角;否则真实方位角为φ0+180°。
步骤十一、根据测量的仪器方位,将反射体的真实方位角(φ0或φ0+180°)从仪器坐标系转换到固定坐标系下。
步骤十二、对最大幅值信号XX(φ0)进行偏移成像。
步骤十三、利用地层横波速度对成像进行时深转换,进一步确定反射体的离井距离,得到最终反射成像结果。
所述步骤一中,两个正交偶极声源X和Y分别指向仪器坐标系的x和y轴方向,其中X(Y)向偶极源由对称分布在仪器x(y)轴两侧的两块极板组成,两块极板激发信号的能量相同、极性相反。仪器工作时,依次激发X和Y向偶极源,向井外地层中辐射弹性波。
所述步骤三中,方位阵列接收器包含轴向上等间隔布置的多组(一般至少八组)方位接收站,每一组方位接收站由对称分布在仪器x和y轴两侧的四块接收极板组成。X和Y向偶极源激发后,每一组接收站记录四个不同方位的共八个分量的声压信号,即XX1、XX2、XY1、XY2以及YX1、YX2、YY1、YY2,其中第一个字母代表声源,第二个字母及数字代表接收器。八个分量信号在处理器中按照下述方式组合,得到四分量正交偶极数据:
其中,m1和m2为调节因子;当m1和m2取不同值时,可实现不同的四分量采集模式。在本发明中,多种采集模式均可以取得良好的方位识别效果,比如:m1=1,m2=1;m1=0,m2=1;m1=-1,m2=2;等。
所述步骤五具体为:首先对四分量波形数据进行滤波处理以减小噪声干扰,然后进行直达波压振处理以减少直达波周期,最后选用中值滤波或F-K滤波等方法进行波场分离,压制和剔除直达波干扰,从全波中提取出反射波信号。
所述步骤六具体为:对步骤五得到的四分量井外反射信号,首先按照共中心点叠加、倾角叠加和径向补偿方法进行信号增强处理,然后用时深转换的偏移成像进行成像,并利用直接滤波或F-K滤波对成像结果做进一步的降噪处理。
所述步骤七中,对于合成反射信号XX(φ),调整φ值得到使XX达到最大幅值时的角度φ0,φ0即为确定的反射体方位;但此时的方位存在多解性,即可能为φ0,也可能为φ0+180°。
所述步骤八中,新信号Ns1和Ns2的合成方式取决于步骤九中理论图版的规律性,本发明中主要介绍了以下两种合成方式(但其有效合成方式不限于此):第一种是直接采用四分量中的交叉分量进行比较,即Ns1=xy,Ns2=yx;第二种是利用最佳成像两侧对称分布的四分量旋转信号进行组合,即Ns1=XX(φ0+φopt),Ns2=XX(φ0-φopt),其中φopt取10°~45°效果较好。得到新信号Ns1和Ns2后,进一步比较它们之间的相位差异(超前或滞后),可通过直接观察波形,或者利用快速傅里叶变换(Fast Fourier Transform,FFT)或动态时间规整(Dynamic Time Warping,DWT)等方法进行计算。
步骤七和步骤八可用反射数据或其时深转换后的成像数据,视实际测井中两种数据信噪比孰为优劣而定。这是因为偏移成像是一线性变换,不改变Ns1和Ns2之间的相对差异。
所述步骤九具体为:根据井孔与弹性波的相互作用理论,考虑充液井内偶极声源在井外的辐射、井外地质体的反射以及充液井孔对反射波的接收调制三个过程,建立充液井中偶极横波方位远探测接收声场的解析解;在此基础上,按照步骤三中的四分量采集模式和步骤八中的新信号合成方式,通过理论模拟,得到新信号Ns1和Ns2之间相位差异随井外反射体方位变化的规律图版。需要指出的是,建立理论图版时,只需要根据四分量反射数据合成信号Ns1和Ns2即可,而无需将理论四分量反射数据进行偏移成像。同上所述,成像运算是线性的,不会改变两个信号之间的相位差异,从反射数据和成像数据中得到的理论图版的规律是一致的。
所述步骤十二具体为:对步骤六中确定的最大幅值信号XX(φ0)首先按照共中心点叠加、倾角叠加和径向补偿方法进行信号增强处理,然后用时深转换的偏移成像进行成像,最后利用直接滤波或F-K滤波对成像结果做进一步的降噪处理。
此外,在本发明的基础上衍生出的另外一种方法是:仪器将井下接收的八个分量数据(即XX1、XX2、XY1、XY2以及YX1、YX2、YY1、YY2)都传输到地面,处理时同时考虑不同的四分量组合模式。具体来说,在步骤七中采用m1=1和m2=0模式组合的数据进行四分量旋转,而在步骤八和九中则采用m1=1和m2=1(或m1=0和m2=1;m1=-1和m2=2;等)模式组合的数据进行相位比较。但这需要记录八个分量的数据,数据的采集量翻倍,并且数据的信噪比也有所降低。
本发明的优点与积极效果如下:
1、本发明提供了一种仅利用正交偶极四分量测量就可以精准探测井外地质构造位置的声波测井方法及装置,有效解决了偶极横波远探测中存在的方位180°不确定的技术难题;
2、本发明所述方法处理结果不受地层和反射体性质、以及接收器接收信号类型(即流体声压或径向位移)的影响,具有较强的稳定性和适用性;
3、本发明所述测井装置环向上只需要接收四个方位的数据,并且在井下经过四分量组合后数据量减半,极大地减少了井下仪器与地面系统之间的数据传输时间,有效提高测井效率(与目前常用的八方位接收模式相比,本发明采集数据量只有其1/4);
4、本发明所述测井装置结构相对简单,所需接收换能器数量较少,装置研制成本低、难度小,各接收换能器性能的一致性更易保持,仪器测量结果也更加可靠;
5、本发明所述方法中的四分量采集模式可以有效消除井内单极直达波的干扰(这在偶极测井中始终存在),有效提高反射波信噪比。
综上,本发明所述声波测井装置研制成本低、难度小,测井效率高,采集数据的信噪比高;所述处理方法成像结果准确,稳定性好,适用性强;本发明整体经济性和实用性较强,具有广阔的应用前景。
附图说明
图1为本发明提供的一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法的工作流程图;
图2为本发明提供的一种正交偶极四分量声波测井装置的结构示意图;
图3为本发明所述测井装置中正交偶极声源的横截面示意图;
图4为本发明所述测井装置中方位接收站(任意一组)的横截面示意图;
图5为井外存在地质体时利用本发明中图2所示测井装置进行声波远探测的模型示意图;
图6(a)为本发明中四分量反射波形(地层a和地层b均为硬地层,假设接收换能器测量声压信号);
图6(b)为本发明中旋转四分量反射波形得到的各个方位信号的归一化幅度(地层a和地层b均为硬地层,假设接收换能器测量声压信号);
图6(c)为本发明中合成信号Ns1和Ns2的反射波形(其中Ns1=xy,Ns2=yx;地层a和地层b均为硬地层,假设接收换能器测量声压信号);
图6(d)为本发明中合成信号Ns1和Ns2的反射波形(其中Ns1=XX(φ0+20°),Ns2=XX(φ0-20°);地层a和地层b均为硬地层,假设接收换能器测量声压信号);
图7(a)、图7(b)分别为本发明中合成信号Ns1和Ns2之间相位差异随井外反射体方位变化的两种规律图版(其中Ns1=xy,Ns2=yx;地层a和地层b均为硬地层,假设接收换能器测量声压信号);
图8(a)、图8(b)分别为本发明中合成信号Ns1和Ns2之间相位差异随井外反射体方位变化的两种规律图版(其中Ns1=XX(φ0+20°),Ns2=XX(φ0-20°);地层a和地层b均为硬地层,假设接收换能器测量声压信号);
图9(a)、图9(b)分别为本发明中合成信号Ns1和Ns2之间相位差异随井外反射体方位变化的两种规律图版(其中Ns1=xy,Ns2=yx;地层a和地层b均为软地层,假设接收换能器测量声压信号);
图10(a)、图10(b)分别为本发明中合成信号Ns1和Ns2之间相位差异随井外反射体方位变化的两种规律图版(其中Ns1=xy,Ns2=yx;地层a和地层b均为硬地层,假设接收换能器测量径向位移信号);
图11为利用本发明所述测井方法及装置模拟的井段四分量远探测数据(以变密度图显示);
图12为本发明中井外四分量反射信号的偏移成像结果(以变密度图显示);
图13为利用本发明所述方法进行井外地质体探测的合成实例处理成果图(其中Ns1=xy,Ns2=yx);
图14为利用本发明所述方法进行井外地质体探测的合成实例处理成果图(其中Ns1=XX(φ0+20°),Ns2=XX(φ0-20°))。
图中编号:1、正交偶极声源;2、隔声体;3、方位阵列接收站;4、处理器;5、正交偶极声源;6、方位接收站(任意一组);7、充液井孔;8、反射界面;9、虚源。
具体实施方案
下面结合正交偶极四分量测量进行井外地质构造探测的具体原理,用理论合成实例对本发明所述方法和装置作进一步说明,以使本领域的技术人员可以更好地理解本发明内容并予以实施,但所举实例不作为对本发明的限定。
如图1和图2所示,本发明提供了一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法及装置,具体工作流程如下:
步骤一、将图2所示的正交偶极四分量声波测井仪放置于测量井中,依次激发两个正交的偶极声源X和Y,向井外地层中辐射弹性波。其中,X(Y)偶极源指向仪器坐标系的x(y)轴方向,由对称分布在仪器x(y)轴两侧的两块极板X1和X2(Y1和Y2)组成,两块极板激发信号的能量相同、极性相反,如图3中所示。
步骤二、记录仪器坐标轴(如y轴)相对于固定方向(如地球北极)的方位角。
步骤三、用方位阵列接收器采集信号并按照特定模式组成四分量数据(xx,xy,yx,yy)。其中,方位阵列接收器包含轴向上等间隔布置的多组(一般至少八组)方位接收站,每一组方位接收站由对称分布在仪器x和y轴两侧的四块接收极板组成,如图4所示。X和Y向偶极源激发后,每一组接收站记录四个不同方位的共八个分量的声压信号,即XX1、XX2、XY1、XY2以及YX1、YX2、YY1、YY2,其中第一个字母代表声源,第二个字母及数字代表接收器。八个分量信号在仪器上的处理器中按照下述方式组合,得到四分量正交偶极数据:
其中,m1和m2为调节因子;当m1和m2取不同值时,可实现不同的四分量采集模式。在本发明中,多种采集模式均可以取得良好的方位识别效果,比如:m1=1,m2=1;m1=0,m2=1;m1=-1,m2=2;等。本发明中以m1=0,m2=1这一种采集模式为例进行介绍,其他采集模式也可得到相似的结论,本发明中不再一一列举。
为了利用理论实例对本发明所述方法进行更好地阐述,首先根据井孔与弹性波的相互作用理论,建立充液井中偶极横波远探测四分量接收声场的解析解。图5给出了利用图2中所示的声波测井仪对井旁地质体进行探测的示意图,仪器上的偶极声源向井外辐射弹性波,经反射后入射到井中,被仪器上的探测器接收。在辐射距离远大于波长的远场辐射条件下,图5中X向偶极声源辐射到地层中的SH和SV横波位移势函数χ和Γ的远场渐近解为(Tang,X.M.,J.J.Cao,and Z.T.Wei,2014,Shear-wave radiation,reception,andreciprocity of a borehole dipole source:With application to modeling ofshear-wave reflection survey[J]:Geophysics,79,no.2,T43-T50.)
其中,为辐射横波波数的最速下降解,ω为圆频率,β为地层横波速度,θt为波辐射方向与z轴正向的夹角;S为声源函数谱;R是井中声源到场点之间的辐射距离;D1和F1分别为SH和SV横波的远场振幅系数,由井壁处的边界条件求得。将(2)式中的和对调,便可得到Y向偶极声源辐射的情况。当辐射横波与地层中的反射体作用后,可将反射回井孔的入射横波波场表示为
式中,RD为声源在井外远处的辐射函数,由(2)式的中括号里的部分组成;RF为声波在反射体处的反射系数,可由佐普利兹方程求得;S、RD、RF都是关于圆频率ω的函数;D为声波在地层中总的传播距离,其传播路径为一条折线。将远探测声场看作是声源在反射体向外一侧的镜像点的辐射声场,此时声波的传播路径D转化为一条直线,如图5中所示,从而(3)式中的球面波传播因子可以展开为多极子柱面波叠加的形式(Li,Y.H.,X.M.Tang,H.R.Li,and S.Q.Lee,2021,Characterizing the borehole response for single-wellshear-wave reflection imaging[J]:Geophysics,86,no.1,D15-D26;Tang,X.M.,andC.H.Cheng,2004,Quantitative borehole acoustic methods[M]:Elsevier SciencePublishing.)
式中,和分别为虚源和场点的位置;同时也是反射体所在方位;k和s分别为入射横波的轴向和径向波数,s=(k2-ω2/β2)1/2;In和Kn分别表示第一类和第二类n阶变型贝塞尔函数;n为多极子的阶数;εn为Neumann因子,n=0时εn=1,n>0时εn=2;将(4)式代入(3)式,得到入射SH和SV横波的位移势函数由下式给出
其中,上标i表示入射波场;下标SH和SV表示计算声源的辐射函数RD和反射系数RF时应分别对应SH和SV横波。
(4)式对(3)式的应用将界面反射波转化为来自虚源的球面入射波,进一步,考虑入射横波与井眼的相互作用,通过衔接边界条件,求解矩阵方程,并将波数积分运用最速下降法进行计算,得到入射横波引起的井内流体的声压和径向位移场,如下(具体求解过程见Li,Y.H.,X.M.Tang,H.R.Li,and S.Q.Lee,2021,Characterizing the borehole responsefor single-well shear-wave reflection imaging[J]:Geophysics,86,no.1,D15-D26.)
式中,为流体纵波的径向波数,k0=(ω/β)cos(π-θr)为入射横波波数的最速下降解,θr为入射波与z轴正向之间的夹角;αf为流体纵波声速,ρf为流体密度;此处的D与(3)式和(4)式中等同,只是与井内场点的径向距离r有关的项在r<<r0的条件下被略去了;An和A′n为井中声场的振幅系数。
结合图5所示的远探测模型,在(2)~(6)式中,对于给定的反射体位置依次激发X和Y向偶极声源,并将(6)式中的分别取为0°、90°、180°和270°,便可以模拟XX1、XY1、XX2、XY2以及YX1、YY1、YX2、YY2八个分量数据,再按照(1)式中的采集模式进行组合,最后便得到理论模拟的四分量反射数据。以反射体方位采集模式为m1=0和m2=1为例,其四分量反射波形在图6(a)中给出(其中,反射体设为地层a和地层b之间倾角为60°的界面;地层a的纵、横波速度分别为4000m/s和2300m/s,密度为2500kg/m3;地层b的纵、横波速度分别为4500m/s和3000m/s,密度为2600kg/m3)。
步骤四、将采集的仪器方位数据和四分量数据通过电缆上传至地面控制系统。
步骤五、对四分量数据进行滤波并做波场分离处理,具体包括:首先对四分量波形数据进行滤波处理以减小噪声干扰,然后进行直达波压振处理以减少直达波周期,最后选用中值滤波或F-K滤波等方法进行波场分离,压制和剔除直达波干扰,从全波中提取出来自井外界面的反射信号。
步骤六、对步骤五得到的四分量井外反射信号,首先按照共中心点叠加、倾角叠加和径向补偿方法进行信号增强处理,然后用时深转换的偏移成像进行成像,并利用直接滤波或F-K滤波对成像结果做进一步的降噪处理。
步骤七、在仪器坐标系(x,y)下,以相对于仪器x轴的角度φ为变量,用四分量井外反射成像信号组成φ向发射和接收的合成XX反射波信号XX(φ)=xx·sin2φ-(xy+yx)·sinφcosφ+yy·cos2φ。将图6(a)中的四分量数据,按照该公式进行四分量旋转,图6(b)的极坐标图给出了归一化幅值|XX(φ)|随旋转角φ的变化曲线。图中显示,当旋转角为117°时,曲线取得极大值,正好对应着反射体的方位(真实方位为120°,在实际远探测测井中,3°的微小误差可以忽略)。然而,由四分量旋转确定的方位存在多解性,即可能为117°或者297°。
步骤八、利用四分量反射数据合成两个新信号Ns1和Ns2,并比较它们之间的相位差异。本发明中主要介绍了以下两种合成方式的分析结果(但其有效合成方式不限于此):第一种是直接采用四分量中的交叉分量进行比较,即Ns1=xy,Ns2=yx;第二种是利用最佳成像两侧对称分布的四分量旋转信号进行比较,即Ns1=XX(φ0+φopt),Ns2=XX(φ0-φopt),其中φopt取10°~45°效果较好。得到新信号Ns1和Ns2后,进一步比较它们之间的相位差异(超前或滞后),可通过直接观察波形,或者利用快速傅里叶变换(Fast FourierTransform,FFT)或动态时间规整(Dynamic Time Warping,DWT)等方法进行计算。
步骤七和步骤八可用反射数据或其时深转换后的成像数据,视实际测井中两种数据信噪比孰为优劣而定。这是因为偏移成像是一线性变换,不改变Ns1和Ns2之间的相对差异。对于无噪声干扰的理论数据,可直接比较信号Ns1和Ns2的反射波形即可。图6(c)给出了第一种合成方式(即Ns1=xy,Ns2=yx)的波形,显示信号Ns1超前于Ns2;图6(c)中横坐标为time(ms),纵坐标为Normalized amplitude,实线为Ns1=xy,虚线为Ns2=yx。
图6(d)给出了第二种合成方式(即Ns1=XX(φ0+φopt),Ns2=XX(φ0-φopt),其中φopt取20°)的波形,显示信号Ns1超前于Ns2。图6(d)中横坐标为time(ms),纵坐标为Normalized amplitude,实线为Ns1=XX(137°),虚线为Ns2=XX(97°)。
步骤九、按照步骤三中的四分量采集模式和步骤八中的信号合成方式,通过理论模拟,建立新信号Ns1和Ns2之间相位差异随井外反射体方位变化的规律图版。同上所述,成像运算是线性的,不会改变两个信号之间的相位差异,从反射数据和成像数据中得到的理论图版的规律是一致的。因此建立理论图版时,只需要根据四分量反射数据合成信号Ns1和Ns2即可,而无需将理论四分量反射数据进行偏移成像。采用图5中所示的远探测模型和步骤三中的理论声场分析,模拟当仪器在井中固定、反射体方位从0°变化到360°时,信号Ns1和Ns2之间相位差异的变化规律。同样考虑步骤八中的两种合成方式,第一种合成方式(即Ns1=xy,Ns2=yx)的结果在图7(a)中给出,其他参数和图6(a)中相同。图中空心圆标记表示信号Ns1超前于Ns2,实心圆标记表示信号Ns2超前于Ns1。环向坐标表示反射体方位,径向刻度表示Ns1和Ns2之间相位差的大小(单位:弧度)。为了更加直观地显示,将图7(a)中Ns1和Ns2之间的超前或滞后关系(不考虑相位差的大小)在图7(b)中给出。结果显示,当反射体位于45°~225°之间时,信号Ns1超前于Ns2;而当反射体位于0°~45°或者225°~360°之间时,信号Ns2超前于Ns1。重复图7(a)、图7(b)的模拟,得到第二种合成方式(即Ns1=XX(φ0+20°),Ns2=XX(φ0-20°))的方位响应结果,在图8(a)、图8(b)中给出。结果显示,当反射体位于0°~45°或90°~180°或225°~270°之间时,信号Ns1超前于Ns2;而当反射体位于其他方位时,信号Ns2超前于Ns1。
步骤十、根据理论模拟图版,确定反射体的真实方位角。具体来说,图6(b)的四分量旋转结果显示反射体的方位为117°或者297°,对于第一种合成方式,图6(c)中显示信号Ns1超前于Ns2,结合图7(a)、图7(b)中的理论模拟图版可确定反射体的真实方位应该为117°;对于第二种合成方式,图6(d)中显示信号Ns1超前于Ns2,结合图8(a)、图8(b)中的理论模拟图版可确定反射体的真实方位应该为117°。两种方式均有效消除了步骤七结果中的方位180°不确定性。判断结果与理论设定的120°基本吻合,证明了上述两种方法的正确性。
然而,实际的测井环境是非常复杂的,图7(a)、图7(b)和图8(a)、图8(b)中所示的信号Ns1和Ns2之间相位差异随反射体方位变化的规律图版的适用性需要进一步考察。一方面,在同一口井中,随着地层埋深增加,地层性质是不断变化的,不同井之间的地层性质也有所不同,而反射体性质通常是未知的。我们将图5中的地层a和b都设为软地层来模拟这种情况(其中,地层a的纵、横波速度分别为2074m/s和869m/s,密度为2250kg/m3;地层b的纵、横波速度分别为1100m/s和600m/s,密度为1400kg/m3),其他参数和图7(a)、图7(b)中相同,由此模型计算的结果如图9(a)、图9(b)所示。另一方面,考虑到不同类型的换能器在井孔流体中接收信号类型的区别(即接收声压或径向位移信号),我们假设换能器在井中接收的是流体径向位移,重复图7(a)、图7(b)的模拟,其结果在图10(a)、图10(b)中给出。图9(a)、图9(b)和图10(a)、图10(b)的结果显示出与图7(a)、图7(b)中相同的规律(对于图8(a)、图8(b)所示的第二种合成方式的结果,重复模拟同样具有相同的结论,此处未作展示)。需要说明的是,图7(a)、图7(b)、图9(a)、图9(b)和图10(a)、图10(b)中有个别点不满足这一规律,但实际上,由于测量时仪器的旋转,恰好消除了这些特殊点对整体结果判断的影响。上述分析表明利用本方明所述方法识别井外反射体方位时不受地层和反射体性质,以及接收信号类型的影响,进一步验证了本发明所述方法的普适性和有效性。
步骤十一、根据测量的仪器方位,将反射体的真实方位角(φ0或φ0+180°)从仪器坐标系转换到固定坐标系下。
步骤十二、对步骤六中确定的最大幅值信号XX(φ0)首先按照共中心点叠加、倾角叠加和径向补偿方法进行信号增强处理,然后用时深转换的偏移成像进行成像,最后利用直接滤波或F-K滤波对成像结果做进一步的降噪处理。
步骤十三、获得了反射体所在方位的成像结果后,利用地层横波速度对成像进行时深转换,进一步确定反射体的离井距离,得到最终反射成像结果。
以下,结合具体的合成井段测井实例,进一步说明本发明所述基于正交偶极四分量测量精准探测井外地质构造的声波测井方法及装置的应用效果。模拟井段深度区间为50m,设置反射界面倾角为75°,位于60°方位(相对于北极)。对此模型模拟由图2中的正交偶极测井装置按照m1=0和m2=1的组合模式采集四分量偶极横波远探测阵列数据。偶极声源位于阵列下方3m处,测井声场的深度和时间采样间隔分别为0.1524m和36μs,采样点数为512。模拟所用声源的中心频率为3kHz。图11给出了模拟得到的第一组接收站上的四分量声压数据(xx,xy,yx,yy),第1道中仪器相对于北极的方位曲线显示仪器在测量时不断旋转,其后的4道以变密度图形式显示了四分量数据。从图中可以看出,来自井外地层界面的反射波在强幅度的井孔直达波之后。利用本发明所述方法对该井段四分量数据进行处理,首先对波形进行滤波和波场分离,消除沿井孔传播的直达波信号,提取来自井外界面的反射信号。然后将四分量井外反射信号按照共中心点叠加、倾角叠加和径向补偿方法进行信号增强处理,并用时深转换的偏移成像进行成像,其反射成像结果在图12中给出(以下行波为例)。
利用图12中的四分量成像数据进行反射体方位识别,考虑步骤八中信号Ns1和Ns2的两种合成方式。对于第一种合成方式(即Ns1=xy,Ns2=yx),其方位处理结果如图13所示。首先将反射成像信号按照步骤七中公式进行四分量旋转得到反射体的两条方位曲线,在第2道中分别以实线和虚线给出。两条方位曲线之间相差180°,这正是偶极横波远探测存在的方位180°不确定性。第3道给出了信号Ns1和Ns2之间的相位关系,其中,空心圆标记表示该深度点的信号Ns1超前于Ns2,实心圆标记表示该深度点的信号Ns2超前于Ns1。根据Ns1和Ns2之间相位差异随反射体方位变化的理论图版(即图7(a)、图7(b)、图9(a)、图9(b)和图10(a)、图10(b)),可由第3道中的相位结果判断第2道中两条方位曲线的“真伪”。具体而言,第3道中空心圆标记的深度点对应反射体方位应该在45°~225°之间,实心圆标记的深度点对应反射体方位应该在0°~45°或者225°~360°之间,由此得到反射体在仪器坐标系下的真实方位曲线,在第4道中给出。根据第1道中的仪器方位曲线将反射体的方位校正到大地坐标系下(即相对于北极),在第5道中给出。图中显示,在整个井段区间,反射体的真实方位约为60°,与理论模型相一致,验证了反演结果的正确性。确定了反射体的真实方位后,对该方位上的信号按照共中心点叠加、倾角叠加和径向补偿方法进行信号增强处理,然后用时深转换的偏移成像进行成像,并利用直接滤波或F-K滤波对成像结果做进一步的降噪处理。最后,利用地层横波速度对成像进行时深转换,进一步确定反射体的离井距离,得到第6道中所示的最终反射成像结果。
对于第二种合成方式(即Ns1=XX(φ0+20°),Ns2=XX(φ0-20°)),其方位处理结果如图14所示。首先按照步骤七中的四分量旋转公式得到第2道中的两条方位曲线,即为最强成像所在方位,然后取该方位正负各20°,再次利用四分量旋转公式合成信号Ns1和Ns2。Ns1和Ns2之间的相位关系在第3道中给出,图中标记含义与图13中相同。根据Ns1和Ns2之间相位差异随反射体方位变化的理论图版(即图8(a)、图8(b)),可由第3道中的相位结果判断第2道中两条方位曲线的“真伪”。具体而言,第3道中空心圆标记的深度点对应反射体方位应该在0°~45°或90°~180°或225°~270°之间,实心圆标记的深度点对应反射体方位应该在45°~90°或180°~225°或270°~360°之间,由此得到反射体在仪器坐标系下的真实方位曲线,在第4道中给出。按照图13中的处理方法可得到反射体的真实方位及其成像结果,分别在第5道中和第6道中给出。图中显示,在整个井段区间,反射体的真实方位约为60°,与理论模型相一致,验证了反演结果的正确性。
需要注意的是,图13和图14的第5道中同时显示,极少数深度点处的反射体方位判断有误(约为240°)。这是由于仪器的旋转,使得在这些深度点处时合成信号Ns1和Ns2的幅度和信噪比减少,信噪比降低使得反演结果的准确度下降,这从图7(a)、8(a)、9(a)和10(a)中也可以看出。实际上,这些位置处的结果已经不可靠,实际数据处理中,在对整个深度区间进行反射体方位识别时,应该剔除信噪比低的数据或数据点。
以上所述实例仅是为充分说明本发明而所举的一些优选实例,本发明的保护范围不限于此。本技术领域的技术人员在本发明基础上所作的等同替换或变换,均在本发明的保护范围之内。本发明的保护范围以权利要求书为准。
Claims (12)
1.一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法,采用如下处理步骤:
步骤一、将正交偶极四分量声波测井装置放于测量井中,依次激发两个正交的偶极声源,向井外地层中辐射弹性波;
步骤二、记录仪器坐标轴相对于固定方向的方位角;
步骤三、用方位阵列接收器采集信号并按照特定模式组成四分量数据xx,xy,yx,yy;
步骤四、将采集的仪器方位数据和四分量数据通过电缆上传至地面控制系统;
步骤五、对四分量数据进行滤波并做波场分离处理,压制和剔除沿井孔传播的直达波信号,提取出来自井外界面的反射信号;
步骤六、将从全波中分离出的四分量井外反射信号进行偏移成像;
步骤七、在仪器坐标系下,以相对于仪器x轴的角度φ为变量,用四分量井外反射信号组成φ向发射和接收的合成XX反射波信号XX(φ)=xx·sin2φ-(xy+yx)·sinφcosφ+yy·cos2φ,变化φ值得到使XX达到最大幅值时的角度φ0;
步骤八、利用四分量反射数据合成两个新信号Ns1和Ns2,并比较它们之间的相位差异;
步骤九、按照步骤三中的四分量采集模式和步骤八中的新信号合成方式,通过理论模拟,建立新信号Ns1和Ns2之间相位差异随井外反射体方位变化的规律图版;
步骤十、根据理论模拟图版,判断φ0是否为反射体的真实方位角:若φ0与步骤八中信号Ns1和Ns2之间相位差异的关系符合步骤九中的理论图版规律,则φ0为真实方位角;否则真实方位角为φ0+180°;
步骤十一、根据测量的仪器方位,将反射体的真实方位角φ0或φ0+180°从仪器坐标系转换到固定坐标系下;
步骤十二、对最大幅值信号XX(φ0)进行偏移成像;
步骤十三、利用地层横波速度对成像进行时深转换,进一步确定反射体的离井距离,得到最终反射成像结果。
2.根据权利要求1所述的一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法,所述步骤一中,两个正交偶极声源X和Y分别指向仪器坐标系的x和y轴方向,其中X和Y向偶极源由对称分布在仪器x和y轴两侧的两块极板组成,两块极板激发信号的能量相同、极性相反;仪器工作时,依次激发X和Y向偶极源,向井外地层中辐射弹性波。
4.根据权利要求1所述的一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法,所述步骤五具体为:首先对四分量波形数据进行滤波处理以减小噪声干扰,然后进行直达波压振处理以减少直达波周期,最后选用中值滤波或F-K滤波等方法进行波场分离,压制和剔除直达波干扰,从全波中提取出反射波信号。
5.根据权利要求1所述的一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法,所述步骤六具体为:对步骤五得到的四分量井外反射信号,首先按照共中心点叠加、倾角叠加和径向补偿方法进行信号增强处理,然后用时深转换的偏移成像进行成像,并利用直接滤波或F-K滤波对成像结果做进一步的降噪处理。
6.根据权利要求1所述的一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法,所述步骤七中,对于合成反射信号XX(φ),调整φ值得到使XX达到最大幅值时的角度φ0,φ0即为确定的反射体方位;但此时的方位存在多解性,即可能为φ0,也可能为φ0+180°。
7.根据权利要求1所述的一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法,所述步骤八中,新信号Ns1和Ns2的合成方式取决于步骤九中理论图版的规律性,包括以下两种合成方式:第一种是直接采用四分量中的交叉分量进行比较,即Ns1=xy,Ns2=yx;第二种是利用最佳成像两侧对称分布的四分量旋转信号进行组合,即Ns1=XX(φ0+φopt),Ns2=XX(φ0-φopt),其中φopt取10°~45°效果较好;得到新信号Ns1和Ns2后,进一步比较它们之间的相位差异,即超前或滞后,通过直接观察波形,或者利用快速傅里叶变换或动态时间规整方法进行计算;
步骤七和步骤八用反射数据或其时深转换后的成像数据,视实际测井中两种数据信噪比孰为优劣而定;因为偏移成像是一线性变换,不改变Ns1和Ns2之间的相对差异。
8.根据权利要求1所述的一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法,所述步骤九具体为:根据井孔与弹性波的相互作用理论,考虑充液井内偶极声源在井外的辐射、井外地质体的反射以及充液井孔对反射波的接收调制三个过程,建立充液井中偶极横波方位远探测接收声场的解析解;在此基础上,按照步骤三中的四分量采集模式和步骤八中的新信号合成方式,通过理论模拟,得到新信号Ns1和Ns2之间相位差异随井外反射体方位变化的规律图版;建立理论图版时,只需要根据四分量反射数据合成信号Ns1和Ns2即可,而无需将理论四分量反射数据进行偏移成像;因为成像运算是线性的,不会改变两个信号之间的相位差异,从反射数据和成像数据中得到的理论图版的规律是一致的。
9.根据权利要求1所述的一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法,所述步骤十二具体为:对步骤六中确定的最大幅值信号XX(φ0)首先按照共中心点叠加、倾角叠加和径向补偿方法进行信号增强处理,然后用时深转换的偏移成像进行成像,最后利用直接滤波或F-K滤波对成像结果做进一步的降噪处理。
10.根据权利要求1所述的一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法,在本发明的基础上衍生出的另外一种方法是:仪器将井下接收的八个分量数据XX1、XX2、XY1、XY2以及YX1、YX2、YY1、YY2都传输到地面,处理时同时考虑不同的四分量组合模式;具体来说,在步骤七中采用m1=1和m2=0模式组合的数据进行四分量旋转,而在步骤八和九中则采用m1=1和m2=1,或m1=0和m2=1,或m1=-1和m2=2模式组合的数据进行相位比较;但这需要记录八个分量的数据,数据的采集量翻倍,并且数据的信噪比也有所降低。
11.根据权利要求1所述的一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法,实现该方法的一种基于正交偶极四分量测量精准探测井外地质构造的声波测井装置,其特征在于,包括以下部分:正交偶极声源发射器短节、发射和接收之间的隔声短节、四个周向等间隔布置的方位单元接收站阵列短节和处理器短节。
12.根据权利要求11所述的一种基于正交偶极四分量测量精准探测井外地质构造的声波测井装置,其特征在于,所述正交偶极声源发射短节包括两个正交布置的偶极声源,声源指向相互垂直;所述接收站阵列短节包括至少八个轴向上等间隔布置的方位接收站;所述方位接收站包括四个周向上等间隔布置的接收换能器;所述处理器短节包括仪器方位测量装置以及与所述声源发射短节、接收站阵列短节和地面控制系统进行通信连接的模块。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210355648.5A CN114779346A (zh) | 2022-04-06 | 2022-04-06 | 一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210355648.5A CN114779346A (zh) | 2022-04-06 | 2022-04-06 | 一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114779346A true CN114779346A (zh) | 2022-07-22 |
Family
ID=82428131
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210355648.5A Withdrawn CN114779346A (zh) | 2022-04-06 | 2022-04-06 | 一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114779346A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115234225A (zh) * | 2022-07-29 | 2022-10-25 | 中国石油大学(华东) | 一种声波远探测数据质量检测方法 |
-
2022
- 2022-04-06 CN CN202210355648.5A patent/CN114779346A/zh not_active Withdrawn
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115234225A (zh) * | 2022-07-29 | 2022-10-25 | 中国石油大学(华东) | 一种声波远探测数据质量检测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2011326567B2 (en) | System and method for investigating sub-surface features of a rock formation | |
US7035165B2 (en) | Imaging near-borehole structure using directional acoustic-wave measurement | |
US8154951B2 (en) | Model-based relative bearing estimation of three-component receivers | |
US9772419B2 (en) | Decomposing full-waveform sonic data into propagating waves for characterizing a wellbore and its immediate surroundings | |
US8797825B2 (en) | Method and apparatus for measuring formation anisotropy while drilling | |
CN112068206B (zh) | 一种基于井孔与弹性波相互作用理论的邻井探测方法 | |
EP1062526A1 (en) | Single well system for mapping sources of acoustic energy | |
CN109143374B (zh) | 一种井周散射体成像方法及系统 | |
US4713968A (en) | Method and apparatus for measuring the mechanical anisotropy of a material | |
Daley et al. | Orbital vibrator seismic source for simultaneous P-and S-wave crosswell acquisition | |
CN109407144A (zh) | 一种基于多波的单孔孤石立体探测方法 | |
Yang et al. | Numerical simulation of acoustic reflection logging while drilling based on a cylindrical phased array acoustic receiver station | |
CN108957548A (zh) | 一种多波多分量联合观测地震页岩气富集区预测技术 | |
Li et al. | Eliminating the azimuth ambiguity in reflected S-wave imaging logging based on the azimuthal receiver mode | |
Gu et al. | Delineating a cased borehole in unconsolidated formations using dipole acoustic data from a nearby well | |
CN105044779B (zh) | 基于相控接收指向性的反射界面方位定量判定方法及装置 | |
CN112647932B (zh) | 一种井中三维探测的多分量声波测井装置 | |
CN103558637A (zh) | 基于三分量传感器的远探测方法 | |
Che et al. | Experimental study on the performance of an azimuthal acoustic receiver sonde for a downhole tool | |
US6176344B1 (en) | Method and system for determining the azimuth position and distance of a reflecting subsurface formation | |
CN207829030U (zh) | 桩基浇注前桩底隐患探查装置 | |
CN114779346A (zh) | 一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法及装置 | |
Zhao et al. | Borehole azimuthal acoustic imaging using 3D spatial scanning: Application in acoustic detection of nearby wells | |
CN114779345A (zh) | 一种基于偶极声波测井仪偏心测量消除井外地质体方位不确定性的方法 | |
Zhao et al. | Numerical simulation of borehole 3D scanning acoustic imaging using scattered waves |
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 | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20220722 |