1、测 井 技 术WELL LOGGING TECHNOLOGYVol.46 No.6 Dec 2022第46卷 第6期 2022年12月文章编号:10041338(2022)06069605相控雷达成像测井仪器中收发天线系统设计李国英1,嵇成高1,于刚刚2,关浩3(1中国石油集团测井有限公司物资装备公司,天津 300280;2中国石油集团测井有限公司天津分公司,天津 300280;3.中国石油集团测井有限公司长庆分公司,陕西 西安 723507)摘要:复杂储层精细勘探需要了解多维几何空间更多的地质信息。相控雷达成像测井技术可以利用雷达反射波成像来展现井周较远范围内裂缝带的空间展布,可通过提取地层
2、介电常数、电导率信息识别复杂油气层。通过分析超宽带雷达波的瞬态脉冲在不同地层、不同流体中的响应机理,确定相控雷达成像测井仪器的基本参数。应用仿真优化设计出一种含有全向偶极子发射天线和定向方位阵列接收天线的优化天线系统,实现径向探测深度8 m、方位分辨率120。相控雷达成像测井技术的研发可为复杂储层油气藏勘探开发提供技术支持。关键词:相控雷达;成像测井;仿真设计;介电常数;天线设计中图分类号:P631.84 文献标识码:ADoi:10.16489/j.issn.10041338.2022.06.009Design of Transceiver Antenna System in Phased R
3、adar Imaging logging ToolLI Guoying1,JI Chenggao1,YU Ganggang2,GUAN Hao3(1.Material Equipment Company,China National Logging Corporation,Tianjin 300280,China;2.Tianjin Branch,China National Logging Corporation,Tianjin 300280,China;3.Changqing Branch,China National Logging Corporation,Xian,Shaanxi 72
4、3507,China)Abstract:Detailed exploration of complex reservoirs requires more information from multidimensional geometry.Phased radar imaging logging technology can show the spatial distribution of fracture zone in the far range around the well by using radar reflection wave imaging.The complex oil a
5、nd gas layers can be identified by extracting the information of layer permittivity and conductivity.In this paper,the basic parameters of phased radar imaging logging instrument are determined by analyzing the response mechanism of ultra bandwidth(UWB)radar wave transient pulse in different layer a
6、nd different fluids.An optimized antenna system with omnidirectional dipole transmitting antenna and directional azimuth array receiving antenna is designed by simulation optimization.The radial detection depth is 8 m and the azimuth resolution is 120.Phased radar imaging logging technology can prov
7、ide technical support for the exploration and development of complex reservoirs.Keywords:phased radar;imaging logging;simulation design;dielectric constant;antenna design0 前 言随着石油开采的深入,碳酸盐岩勘探力度也在逐步加大1。利用现有测井技术可对碳酸盐岩储层类型、非均质性、展布情况等进行综合评价,但都很难给出井周远距离范围内精细的、定量的、多维几何空间的地层信息2。中国的碳酸盐岩勘探覆盖面广,西部塔里木碳酸盐岩油气藏、西
8、南的滩相碳酸盐岩油气藏、北部鄂尔多斯的奥陶系碳酸盐岩气藏及东部大港油田千米桥潜山凝析气藏都有很不错的勘探前景3。随着中国油气资源开采难度越来越大,勘探对象已经逐步由传统构造性油气藏过渡到岩性油气藏、非常规油气藏,因此,测井技术的发展必须适应勘探开发需求。相控雷达测井成像技术可以利用雷达反射波成像来探测裂缝带的空间展布,可以从雷达波中提取地层介电常数、电导率信息来识别复杂油气层4。雷达测井在应用初期工作频率不超过1 000 MHz,频率较低。它主要是对井周较远范围内的地质体进行探测,分辨率较高,通过数据成像实现对井周地质体结构的评价分析5。2002年,斯伦贝谢公司的首席顾问Chen等验证了井中雷
9、达可用于油气资源勘探的可行性,初步估计单孔井中雷达在油井中的径向探测深度为0.510.0 m1-2。相应地也开发了一些商业化的雷达测井仪器,但是这些仪器只适用于浅地层探测,第一作者:李国英,女,1971 年生,高级工程师,从事测井井下仪器的维修保养工作。E-mail:李国英,等:相控雷达成像测井仪器中收发天线系统设计第46卷 第6期697 井下探测器部分尚不能承受油井下的高温高压环境。2012年,中国研发出瞬态脉冲雷达成像测井原理样机,并在作业井中进行了首次试验测试,测井深度达1 596 m,径向探测深度达5.0 m,初步实现了对井周异常目标体的定位5。但是,由于钻井液和井周地层的衰减作用,限
10、制了该原理样机的径向探测深度。本文设计的相控雷达成像测井仪器是基于时域脉冲技术,通过控制各天线元件的相位,保证天线的非色散特性,有效减小时域脉冲信号通过天线接收或发射的失真,提高系统的分辨率和瞬时功率。该仪器能够获取更丰富的井周原状地层信息,为复杂储层地质评价提供更精细的数据,具有广阔的市场空间。1 瞬态脉冲波传输特性理论分析1.1 工作频率的选取雷达波在损耗型介质中传播时,会发生幅度和相位的衰减。其幅度和相位的衰减与雷达波在损耗介质中传播时衰减常数、介质的磁导率、电导率、介电常数有关。瞬态脉冲是一种时域正演,具有快速、准确的特点,可满足测井快速分析需求。岩层中,瞬态脉冲的速度表达式为=+-2
11、1121 2/(1)式中,为传播速度,m/s;为介质磁导率,H/m;为点电导率,S/m;为介电常数,F/m;为角频率,rad*/s。衰减常数和速度的简化表达式为 20r(2)c0r (3)式中,0为自由空间的波阻抗,通常取377;r为相对介电常数;c0为自由空间中的波速,通常取3108 5。因此,对于大多数地质材料而言,在高频或低电导率的情况下,瞬态脉冲的波速仅取决于介质的相对介电常数。把常见地质材料的相对介电常数和电导率值代入雷达波衰减公式分析发现,雷达波频率在100 MHz时,其相位和幅度的衰减随着频率的增加可以减小到一定数值。而在较低频段,相位离散却非常严重6。如果工作频率高于100 M
12、Hz,在/1的岩层中,瞬态脉冲的衰减量(ATT)和相位离散值*非法定计量单位,1 rad=(180/)()有效电导率/(Sm-1)花岗岩石灰岩花岗岩石灰岩频率/GHz 频率/GHz 有效相对介电常数100908070605040302010016141210864201001010.10.010.0010.000 10.000 010.000 0010.000 000 1图1 频率与介电常数关系图图2 频率与电导率关系图(PHerr)与频率无关(ATT=343.4/,PHerr=1),雷达脉冲的传输没有离散,而且其衰减量基本达到一个固定值。因此,/1是雷达测井正常工作的前提条件5。本文利用雷达
13、波探测高速岩层,分析石灰石和花岗岩等典型岩石的介电常数、电导率、波速和衰减系数与工作频率的关系特性7。从图1图4可以看出,在频率低于1 GHz时,这2种介质的有效相对介电常数近似为常数。工作频率在500 MHz以下时,相对导电率近似为常数。由于弛豫效应,当工作频率大于1 GHz时,电磁波在花岗岩和石灰石中的衰减急剧增加;但工作频率小于1 GHz时,电磁波在花岗岩和石灰石中的衰减很小。工作频率在100 MHz1 GHz时,波速近似为常数;而工作频率低于100 MHz时,波速的变化会导致脉冲轻微色散。综合考虑系统的中心频率、工作带宽与地质特征等,地层对电磁波的衰减与色散特性决定了工作频率适合100
14、 MHz500 MHz。为尽量增大径向探测距离,应选择较低的中心频率,但探头的尺寸与所辐射电磁信号的波长相关,受到井下空间限制,同时兼顾系统分辨率,并考虑在高温高压环境下过高的发射脉冲及接收频率会对系统的散热带来较大的负担,因此,按照中心频率230 MHz对雷达测井仪器进行设计。2022年测 井 技 术698 1.2 超宽带电磁波与油气水的响应机理分析为了解超宽带电磁波在油、气、水中的传播衰减情况,建立一种数值模型,分析其在碳酸盐岩裂缝结构中的响应机理。图5分别显示了在碳酸盐岩层裂缝结构中填充空气、油和淡水时的测量信号,图5中可以看到模拟测量处不同介质的成像,其中淡水最明显,其次是空气,最后是
15、油。(a)淡水(b)空气(c)油1.51.00.50.0-0.5-1.0-1.5-2.0-2.50.240.160.080.00-0.08-0.16-0.240.200.150.100.050.00-0.05-0.10-0.15图5 碳酸盐岩裂缝充填淡水、空气、油的回波信号2 相控雷达测井仪器天线设计天线是相控雷达测井仪器中最重要的组成部分之一,其在雷达设备中主要有能量转换功能和定向辐射(或接收)功能8。雷达波在地层中传播、衰减后被接收探头所接收,通过对接收信号的处理分析,可获得雷达波所经过储层的信息,进而实现对储层的精细评价。2.1 发射天线设计根据相控雷达成像测井仪器的径向探测深度及方位探
16、测性能,设计采用全向偶极子天线作为发射天线。设计的全向偶极子发射天线具有2个臂,同时臂间有加载电阻;外围设计介质套层以实现与钻孔环境介质的匹配(见图6)。综合分析探头的内外径变化,以及不同馈电点缝隙宽度对反射系数S11和辐射电磁脉冲幅度的影响情况(见图7、图8),设计探头半径为16 mm,探头长度为70 cm,馈电点间距为18 mm。设计发射天线的电路控制采用高幅值、低抖动、小拖尾的纳秒级瞬时脉冲激励信号。该激励信号的产生采用雪崩管Marx电路的方式,形成多路双极性脉冲子源,经过宽带的相干功率合成器将多路脉冲子源进行功率叠加,叠加后的瞬时功率增大。大功率瞬时脉冲最终经匹配网络馈入天线。铝质金属
17、管加载电阻图6 发射天线模型图图3 频率与衰减系数关系图图4 频率与波速关系图波速/(mns-1)衰减系数/(dBm-1)花岗岩石灰岩花岗岩石灰岩频率/GHz 频率/GHz 1001010.10.011001010.10.010.0010.11100.1900.1800.1500.1450.1400.1300.1100.1050.1000.0500李国英,等:相控雷达成像测井仪器中收发天线系统设计第46卷 第6期699 2.2 接收天线设计接收天线设计为定向的单极子天线阵。每个单极子天线的辐射片为1个水滴型金属片,反射板底端与金属地相连,上端通过一定数值的电容和电阻与辐射片相连。分析230 M
18、Hz单极子接收天线在模型井中频率与S11参数的变化情况。模型井内的井液设置为泥浆,相对介电常数为78,直径为200 mm;井周岩层的相对介电常数为8、相对磁导率为1。其衰减变化特性见图9,S11参数 6 dB的频率范围为100261 MHz,与设计主频230 MHz相符合。接收天线共设计3个,均匀分布在不同的轴向上,角度相差120。1.00.80.60.40.200-10-20-30频率/GHzS11/dBr2=14 mm r2=16 mm r2=18 mm r2=20 mm 图8 馈电点缝隙宽度(L3)对S11参数影响图7 不同外径(r2)对探头S11参数的影响L3=15 mmL3=12 m
19、mL3=18 mmL3=21 mm频率/GHz0.80.70.60.50.40.30.20.100-2-4-6-8-10-12-14-16-18-20-22S11/dB8007006005004003002001001200-1-2-3-4-5-6-7-8-9-10-11频率/GHzS11/dB1 (100.57,-5.9635)2 (261.94,-5.9947)图9 天线在模型井中S11参数随频率的衰减变化关系2.3 收发天线组合仿真分析根据设计的收发天线特性,仿真模拟单发三收组合天线在模型井中的信号收发情况。设计二号天线正对目标体,分析每个不同方位单极子信号波形。二号天线接收到波形幅度最
20、大,因其正对目标体,信号最强。一号天线与三号天线接收到信号强度相等,仿真结果与预期结果一致。对单发双收和单发三收的接收波形进行对比,得到结果大致相同,表明单极子定向接收天线之间间距为200 mm(一号天线与二号天线、二号天线与三号天线之间间距均为200 mm)时能够实现最小的耦合。确定单极子定向接收天线的最佳组阵间距为200 mm。接收天线阵列中一号天线主辐射方向为0,二号天线主辐射方向为120,三号天线主辐射方向为240。采用算法矫正的3个不同平面的天线可等效于同一平面主辐射方向不同的3个天线。因此,可根据回波的相位信息和幅度信息,利用相控技术反推目标的角度信息,进而实现目标的三维定位。3
21、应用效果 相控雷达成像测井仪器研制完成后,在大港油田的双X井中进行下井试验。该井目的层为潜山奥陶系,完钻井深3 318.0 m。从雷达成像图和能量曲线分析可以看出在2 950.02 958.0 m处存在井旁裂缝,径向距离在4.05.0 m。从介电常数曲线看,在均匀碳酸盐岩层段其平均值约为9,在裂缝段增大,而裂缝段的相对电导率也增大,反映裂缝内存在介电常数大而电阻率低的填充物。在同深度段常规自然伽马值较高,而无铀伽马值较低(与该井段上下地层基本一致),说明铀含量相对较高,在这段厚度较大且岩性较均质的碳酸盐岩地层中出现铀含量较高的地层,反映该地层发育裂缝或孔洞,曾有流体活动(见图10)。2022年
22、测 井 技 术700 图10 雷达成像测井与常规资料对比图*非法定计量单位,1 b/eV=6.241 46 1010 m2/J测井解释结论深度/m垂直深度/m100气测全烃/%6010井径/cm1500自然伽马/API1500无铀伽马/API1002 9302 9402 9502 9602 9702 9301091101112 9402 9502 9602 9702200 0002200 0002200 0002200 0002200 0000钻时/min录井剖面阵列侧向测井(MLR1C)/(m)阵列侧向测井(MLR2C)/(m)阵列侧向测井(MLR3C)/(m)阵列侧向测井(MLR4C)/(
23、m)球形聚焦电阻率/(m)-1545补偿中子/%32补偿密度/(gcm-3)100600补偿声波/(sm-1)200光电吸收截面指数/(beV-1)*01000100含水饱和度/%050残余油可动油砂岩白云岩石灰岩泥岩孔隙度/%孔隙度/%(b)常规测井资料解释图141210864200.207.0 7.75 8.5 9.25 10 10.750.42 9202 9302 9402 9502 9602 9702 9802 9202 9302 9402 9502 9602 9702 9802 9202 9302 9402 9502 9602 9702 980归一化回波能量距离/m相对介电常数深度/
24、m深度/m深度/m(a)雷达成像测井图(下转第 706 页)2022年测 井 技 术706 5 结 论(1)声波测井换能器的结构比较复杂,难以通过解析求解的方法进行描述,有限元数值模拟已经成为换能器优化设计的基本手段。通过有限元分析软件数值模拟边界条件、结构尺寸和压电材料对换能器的影响,得到的结果可服务于声波测井换能器的设计和研制。(2)通过声波测井声场传播理论研究,开展不同地质条件下声波测井正演模拟,可明确声波测井换能器频率与激发能量、换能器灵敏度等关键技术参数。(3)通过基于通用软件的换能器有限元数值模拟,建立符合实际仪器结构的有限元模型,进行声波测井换能器多物理场耦合数值模拟,提升关键参
25、数模拟计算精度,优化计算模型,开展对比计算分析,可得到换能器详细结构设计参数。随着勘探难度的增大,井深增加,大斜度、水平井增多,高性能换能器的研究与开发势在必行。参考文献:1 汤天知,岳爱忠,朱益华,等.测井数值模拟平台设计与实现J.测井技术,2010,34(6):585-588.2 中仿科技公司.COMSOL Multiphysics有限元法多物理建模与分析M.北京:人民交通出版社,2007.3 乔文孝,车小花,张菲.边界条件对声波测井弯曲振子振动模态的影响J.中国科学(D辑),2008,38(z1):174-179.4 陈秋颖,丛健生,王秀明.流体负载对偶极声波测井换能器性能影响的数值分析
26、C/中国科学院声学研究所第四届青年学术交流会,北京,2012.5 张正鹏,刘玉凯,苏远大,等.考虑压电声源-井孔系统的随钻方位声波测井数值模拟J.测井技术,2020,44(1):1-7.6 唐晓明,郑传汉.定量测井声学M.赵晓敏,译.北京:石油工业出版社,2004.7 张海澜,王秀明,张碧星,等.井孔的声场和波M.北京:科学出版社,2004.8 栾桂冬,张金铎,王仁乾.压电换能器和换能器阵M.北京:北京大学出版社,2005.9 李淑君,王惠泉,赵文玉,等.基于COMSOL多物理场耦合仿真建模方法研究J.机械工程与自动化,2014(4):19-20,23.10 付勇,陈晔,张伟民.基于多物理场耦
27、合的夹心式压电换能器优化设计J.轻工机械,2018,36(5):1-8,13.11 陈雪莲,唐晓明.孔、裂隙并存地层中的声波测井理论及多极子声场特征J.地球物理学报,2012,55(6):2129-2140.12 黄飞,张炳军,王国平,等.数值模拟驱动的过钻具偶极声波测井仪设计J.测井技术,2022,46(2):158-163.13 张奔牛.传感和致动装置/结构中的多物理量耦合问题D.重庆:重庆大学,2003.(修改回稿日期:2022-09-11 编辑 林丽丽)(上接第 700 页)4 结 论(1)通过一系列的理论及数值模拟研究,确定相控雷达成像测井仪器工作频率及天线工作参数。(2)研究探头的
28、外径、馈电点间距、加载电阻的阻值及时域波形和S11参数的影响等,完成了全向发射探头和定向接收探头的设计。(3)设计研发出的相控雷达成像仪器经过现场测井试验,实现对井旁较远距离的小型地质构造探测,同时通过地层介电常数、电导率信息的定性分析,可对储层性质进行评价。该技术可以作为解决识别与评价裂缝性储层性质难题的手段之一,具有较好的应用前景。参考文献:1 MA C G,ZHAO Q,RAN L M,et al.Numerical study of borehole radar for cliff imaging J.Journal of Environmental and Engineering G
29、eophysics,2014,19(4):269-276.2 CHEW W C.Waves and fields in inhomogeneous media M.New York:Wiley-IEEE Press,1999.3 徐建华.层状媒质中的电磁场与电磁波M.北京:石油工业出版社,1997.4 粟毅,黄春琳,雷文太.探地雷达理论与应用M.北京:科学出版社,2006.5 马春光.瞬态脉冲雷达成像测井及实验研究D.成都:电子科技大学,2015.6 王小卫.海相碳酸盐岩非均质储层成像研究:以塔里木盆地碳酸盐岩资料处理为例D.武汉:中国地质大学(武汉),2011.7 吴优,息荣艳,潘笑天,等.基于相控电磁表面的W波段雷达系统研究J.雷达学报,2021,10(2):281-287.8 王诚至.基于震电测井的信号收发电路系统研究D.成都:电子科技大学,2022.(修改回稿日期:2022-10-08 编辑 林丽丽)