CN104481525B - 一种快测阵列感应测井系统及其测井方法 - Google Patents

一种快测阵列感应测井系统及其测井方法 Download PDF

Info

Publication number
CN104481525B
CN104481525B CN201410714371.6A CN201410714371A CN104481525B CN 104481525 B CN104481525 B CN 104481525B CN 201410714371 A CN201410714371 A CN 201410714371A CN 104481525 B CN104481525 B CN 104481525B
Authority
CN
China
Prior art keywords
circuit
unit
data
array
leather bag
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
Application number
CN201410714371.6A
Other languages
English (en)
Other versions
CN104481525A (zh
Inventor
陈涛
白彦
陈章龙
宋青山
范晓文
王水航
师光辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
China Petroleum Logging Co Ltd
Original Assignee
China National Petroleum Corp
China Petroleum Logging Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China National Petroleum Corp, China Petroleum Logging Co Ltd filed Critical China National Petroleum Corp
Priority to CN201410714371.6A priority Critical patent/CN104481525B/zh
Publication of CN104481525A publication Critical patent/CN104481525A/zh
Application granted granted Critical
Publication of CN104481525B publication Critical patent/CN104481525B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing 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)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种快测阵列感应测井系统及其测井方法,包括经由高性能集成电子线路、线圈系、压力平衡装置组成的井下测量装置,测井得到的快测阵列感应测量信号输送至地面采集模块,经过数据处理模块进行数据合成处理后通过记录仪显示单元显示,从而完成快测阵列感应仪器系统的信号测量。通过线圈系参数与结构优化、高性能集成电路系统、工作模式优化、厚膜集成与一体化设计等技术,在保证仪器纵向分辨率和径向探测深度不变的基础上,仪器长度由原来9.86m缩短到4.96m,测量精度由原来±1mS/m提高到±0.5mS/m。仪器现场应用表明,仪器稳定可靠,可操作性、易维护性大大提高。

Description

一种快测阵列感应测井系统及其测井方法
技术领域
本发明涉及一种用于快测平台探测裸眼井周围地层电阻率的快测阵列感应成像测井系统及其测井方法,属于快测平台电法领域的地球物理测井技术。
技术背景
随着油田开发方式的多样化,油藏状况越来越复杂,国外三大测井公司斯伦贝谢、贝克—阿特拉斯、哈里伯顿均先后推出了自己的快测平台Express Platform、LOG-IQ、FOCUS测井系统。随着国内钻井提速,井况越来越恶劣,研制长度短、时效快的快速测井平台成为油田勘探开发提速提效的迫切需求。国产EILog快测平台正是在此背景下开发,通过探测器优化、电路共用与集成,使用新的机电工艺设计,缩短了仪器长度,增强了测量的稳定性。目前,先进的快测阵列感应测井仪器应用广泛,它以电磁感应原理为理论依据,主要由电子线路、线圈系,电路共用与集成、计算机技术及数据处理等先进方法,通过多路遥测短节,把采集的大量数据送到地面,再经过计算机进行处理,得出具有5种不同探测深度和3种不同纵向分辨率曲线。快测阵列感应测井仪器的主要组成部分是线圈系结构,线圈系采用一个发射线圈、8个接收线圈阵列的三线圈系源距和工作频率与公开专利CN102042009 A线圈系参数一致,以保证仪器原三种纵向分辨率和五种径向探测深度不变。但是,常规的阵列感应测井仪线圈系芯轴,选用无磁的钛合金材料。
发明内容
为了高效地满足快测平台现场应用要求,本发明实施例提供了一种石油勘探设备快测阵列感应测井系统及其测井方法,仪器通过线圈系参数与结构优化、高性能集成电路系统、工作模式优化、厚膜集成与一体化设计等技术,在保证仪器纵向分辨率和径向探测深度不变的基础上,仪器长度由原来9.86m缩短到4.96m,测量精度由原来±1mS/m提高到±0.5mS/m。仪器现场应用表明,仪器稳定可靠,可操作性、易维护性大大提高。
本发明的目的是通过下述技术方案来实现的:
一种快测阵列感应测井系统,该系统包括:
高性能集成电子线路装置,包括厚膜集成电路单元模块和高性能采集处理单元模块,用于产生三电平控制信号提供给发射电路,启动发射控制电路;
线圈系,包括线圈系芯轴和线圈骨架与线圈屏蔽层,用于支撑阵列感应仪器线圈系总成;
压力平衡装置,采取皮囊压力平衡的方式与复合芯轴的结构相匹配设计,采用皮囊注油的方式,在井下高温高压环境下保持线圈系玻璃钢内外压力的平衡;
地面采集模块,用于实时采集、监测快测阵列感应仪器测井数据,并完成控制仪器刻度文件操作、车间刻度操作个温度传感器标定操作;
数据处理模块,用于对地面采集模块的原始测井数据进行处理,得到最终用户需求解释的地层电阻率测井曲线;
记录仪显示单元,用于记录显示地面采集模块和数据处理模块的地层信息曲线和图像数据,生成测井资料成果打印和导出、生成图片的功能;
经由高性能集成电子线路装置、线圈系、压力平衡装置组成的井下测量装置,将测井得到快测阵列感应测量信号输送至地面采集模块,经过数据处理模块进行数据合成处理后通过记录仪显示单元显示,从而完成阵列感应仪器系统的信号测量。
进一步地,所述线圈系与高性能集成电子线路装置相连,高性能集成电子线路装置包括发射驱动电路、前置放大和带通滤波电路、数据采集电路和电源管理单元;
线圈系接收由地层产生的测量信号,分别通过发射驱动电路、前置放大和带通滤波电路将处理后的信号传输至数据采集电路,数据采集处理电路将该测量信号进行信号选择、自动增益、A/D采集及数据处理再进行组织编码,并通过CAN接口电路将数据传送至井下遥测。
进一步地,所述发射驱动电路包括二级刻度电路和发射解码与驱动电路,二级刻度电路分别与线圈系和前置放大和带通滤波电路连接,发射解码与驱动电路分别连接线圈系和数据采集电路相连。
进一步地,所述前置放大和带通滤波电路包括若干个前置放大电路和带通滤波电路,各前置放大电路和带通滤波电路分别连接线圈系的8通道和数据采集电路,实现同步采集与处理,以及系统通讯。
进一步地,所述数据采集电路包括与发射解码与驱动电路连接的发射控制波形电路,与若干个前置放大电路和带通滤波电路相连的电平转换电路,各电平转换电路分别依次对应连接8个PGA可编程增益放大器、高速A/D转换电路;所述发射控制波形电路和高速A/D转换电路分别连接FPGA主控单元,FPGA主控单元通过CAN接口连接DSP微处理器;
DSP微处理器分别连接线圈系温度检测电路、自然电位SP、工作电压检测电路和电子仪温度检测电路。
进一步地,所述压力平衡装置包括皮囊支撑轴,及包裹在皮囊支撑轴外侧的皮囊,所述皮囊支撑轴的一端皮囊内镶装有皮囊左堵头,皮囊左堵头与皮囊支撑轴之间设有皮囊芯体(304);在皮囊支撑轴的另一端镶装有套有O型密封圈的皮囊右堵头,皮囊右堵头通过金属扎带连接球形堵头,皮囊左堵头通过金属扎带连接线圈系下管轴。
相应地,本发明还给出了一种快测阵列感应测井系统的测井方法,包括下述步骤:
步骤401,通过地面采集模块,完成实时采集快测阵列感应仪器测井数据,并同时将测井数据送入预处理单元;
步骤501,仪器在井下测量实时采集数据,经过地面采集模块,得到八组线圈系的24条实电压信号,送入预处理单元,调用刻度系数将电压信号转化为转化成实时测量的地层电导率信号;
步骤502,将预处理单元输出的地层电导率信号输入趋肤效应校正单元(502),通过三频趋肤效应校正处理算法,对24条实地层电导率信号进行校正得到8条测井数据曲线SEC1~SEC8;
步骤503,将趋肤效应校正单元得到8条数据曲线,输入到下一级井眼校正处理单元,井眼校正处理,依据快侧阵列感应仪器结构,利用三维正演模拟程序计算建立了大型居中和偏心井眼校正库;通过采用查表的形式调用井眼校正库,对趋肤效应校正单元502的8条曲线做井眼校正处理,完成对实际测井资料的仪器居中和偏心校正处理功能,最终输出8条去除井眼影响的数据曲线BHC1~BHC8;
步骤504,将井眼校正处理单元得到8条去除井眼影响的数据曲线,输入到下一级真分辨率聚焦单元,利用聚焦处理算法单元对8条信息数据进行井眼校正,通过调用事先算好的聚焦滤波器库进行真分辨率聚焦处理,得到5条不同探测深度的电阻率曲线T10、T20、T30、T60、T90;
步骤505,根据真分辨率聚焦单元的电阻率曲线T10、T20、T30、T60、T90数据,输入到下一级垂向分辨率匹配处理单元(505),利用8个子阵列3种频率的24条几何因子响应曲线经过三频趋肤效应校正,得到校正后的8条几何因子曲线,调用事先算好的聚焦滤波器库在相邻两条曲线中高分辨率曲线变换为低分辨率曲线的滤波器,调用分辨率匹配滤波器进行处理得到3组分辨率为0.3m、0.6m和1.2m的曲线,每一组曲线有5条探测深度为0.25m、0.50m、0.75m、1.50m和2.25m的曲线,输出合成处理后的3种分辨率5种探测深度的15条曲线处理,分别命名为AO10~AO90、AT10~AT90、AF10~AF90;
步骤506,最后将分辨率匹配处理后的分辨率为0.6m的曲线,探测深度为0.25m、0.50m、0.75m、1.50m和2.25m的曲线通过径向反演单元(506)进行处理,该反演处理单元采用常规通用的基于阵列感应仪器的3参数反演算法,处理后输出得到真地层电阻率Rt和侵入带电阻率Rxo以及侵入带半径r,即完成快测阵列感应测井过程。
进一步地,所述步骤502中,利用三频趋肤效应校正处理算法,利用下述表达式实现三频趋肤效应校正:
其中wi为第i频率的权重系数,σai为第i频率的视电导率值,E0,E1和E2为选择函数。
进一步地,所述步骤504中,利用聚焦处理算法单元对8条信息数据进行井眼校正:
其中,σpk,z)为合成后径向探测深度ρk,纵向深度z处的电导率;k=1,2,…,K,K是合成探测深度数目,这里表示为T10、T20、T30、T60、T90;为第j个子阵列在z点的测量值;j=1,2,…,J,J是子阵列个数;是聚焦合成到探测深度ρk时,第j个子阵列的滤波器;σ是地层的电导率;z=zmin~zmax是聚焦合成的窗口范围,它是聚焦合成滤波器是地层电导率的函数。
进一步地,所述步骤505中,调用事先算好的聚焦滤波器库在相邻两条曲线中高分辨率曲线变换为低分辨率曲线的滤波器,聚焦滤波器库通过下述方式实现:
设计匹配滤波器hvgfi,使gpvgfi变换为gpvgf(i+1),为:
gpvgf(i+1)(z)=hvgfi(z)*gpvgfi(z)
式中,i=1,2,…,4,gpvgf(i+1)(z)为将要分辨率匹配的目的函数。
步骤503中,井眼校正算法,在原来阵列感应仪器居中井眼校正算法基础上增加了偏心校正处理,实现了自适应居中和偏心校正处理功能
本发明实施提供的技术方案带来的有益效果是:
本发明采用的新型管材为主体、外表面覆盖玻璃钢的复合芯轴,其中芯轴是无磁铍青铜金属管长轴,外壁上有玻璃钢套,将金属管长轴与玻璃钢套复合粘接为一体构成阵列线圈系骨架,并且采用一种组合线圈的方式,将各个阵列的接收线圈和屏蔽线圈与安装的阵列线圈系一一对应,通过特殊三线圈系结构设计,解决了相邻两个接收线圈阵列主接收线圈与屏蔽线圈之间位置干涉技术难题;通过绕线工艺改进与实验优化,解决了快测阵列结构线圈各道直接耦合信号过大技术难题。其次,是快测阵列感应集成采集系统,采用开关电源、发射接收一体化设计、电路厚膜集成技术,去掉保温瓶,减少接头,实现了仪器长度与重量均缩小一半;通过实时二级刻度技术、系统内刻信号幅度与相位同步校正技术,解决了电路的温度漂移技术难题;通过采用高性能DSP与FPGA架构,倍增采集通道数,优化采集时序与DPSD算法,利用更多采集时间和获取更大数据量来提高测量精度,解决了高密度条件下强弱信号的影响及噪声抑制难题,并实现小信号测量精度提高。压力平衡装置,本发明改进了活塞式压力平衡短节,采取皮囊压力平衡的方式与复合芯轴的结构相匹配设计,由有效缩短整体仪器0.5m。
快测阵列感应仪器数据处理软件,采用在井下完成8组线圈3种频率信号的实部和虚部信号提取,得到24条实部和24条虚部电导率曲线,在地面通过软件,经过三频趋肤效应校正后进行自适应井眼校正实现真分分辨率聚焦与分辨率匹配合成处理,得到3种纵向分辨率、5种径向探测深度的电阻率曲线。多个接收线圈的测量信息包含有地层径向和纵向不同范围的地层电阻率信息及井眼信息,经过井眼环境影响校正和数据处理,得到地层径向视电阻率、地层流体性质的二维成像图,用于分析薄层和层内非均质性,描述泥浆侵入特征,判断油水层。
通过采取上述手段开发的适于快测平台的阵列感应成像测井仪已成功应用并开始产业化。仪器系统设计方法与电路一体化设计技术对于其他仪器设计具有重要参考价值。
附图说明
下面结合附图对本发明方法以作进一步的说明:
图1为快测阵列感应成像测井系统;
图2为电子仪主控采集处理单元;
图3为主控采集单元DSP&FPGA交互工作时序;
图4为厚膜集成电路模块单元;
图5为新型优化线圈系结构;
图6为压力平衡结构单元;
图7为地面测井资料数据处理框图;
图8为测井资料解释成果图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
如图1所示,本发明实施例提供一种快测阵列感应测井系统,该系统主要由6个环节组成:
高性能集成电子线路1,包括高性能采集处理单元模块和厚膜集成电路单元模块。用于产生三电平控制信号提供给发射电路,启动发射控制电路。发射信号经驱动后由发射线圈向地层发射,每个接收子阵列线圈接收由地层产生的电磁感应信号,通过高性能采集处理单元模块采集接收线圈产生的地层感应电动势,并将采集数据结果上传至地面,可以实时获得8个通道的主线圈测量信号。
线圈系2,具体包括线圈系芯轴和线圈骨架与线圈屏蔽层,用于支撑阵列感应仪器线圈系总成,其上安装有一体化线圈系结构和特殊的绕线技术的线圈系总成的所有阵列线圈骨架,起到很好的绝缘、抗冲击、减震、承力作用。
压力平衡装置3,采取皮囊压力平衡的方式与复合芯轴的结构相匹配设计,用于皮囊注油的方式,在井下高温高压环境下保持线圈系玻璃钢内外压力的平衡。
地面采集模块4,用于实时采集、监测快测阵列感应仪器测井数据,包括控制仪器刻度和测井。刻度完成三个任务,一是刻度文件操作,二是车间刻度操作,三是温度传感器标定操作;测井用于完成仪器配置、加载刻度及定时测、上测或下测组成,实时监视数据,实时看快测阵列感应仪器测井曲线。
数据处理模块5,用于对地面采集模块的原始测井数据进行处理得到最终用户需求解释的地层电阻率测井曲线。数据处理模块为一直线流程,即上一处理模块的输出为下一处理模块的输入。其中,预处理模块用到刻度系数,实现电压信号到地层电导率信号的转化;井眼校正模块用到井眼校正参数、井径和泥浆电阻率值,用于对测量数据的井眼影响校正。真分辨率聚焦处理和分辨率匹配处理模块用到聚焦滤波器库和匹配滤波器库,用于对上一级井眼校正输出数据进行真分辨率聚焦处理和分辨率匹配处理,最终完成用户需求的测井解释曲线。
记录仪显示单元6,用于记录显示地面采集模块和数据处理模块的地层信息曲线和图像数据,生成测井资料成果打印和导出、生成图片的功能。
经由高性能集成电子线路1、线圈系2、压力平衡装置3组成的井下测量装置,测井得到的快测阵列感应测量信号输送至地面采集模块4,经过数据处理模块5进行数据合成处理后通过记录仪显示单元显示6,从而完成快测阵列感应仪器测量地层信息包含有地层径向和纵向不同范围的地层电阻率信息及井眼信息,经过井眼环境影响校正和数据处理,得到地层径向视电阻率。
如图2所示,本发明实施例提供一种高性能集成电子线路,具体包括与线圈系2相连的高性能集成电子线路1。高性能集成电子线路1包括发射驱动电路102、前置放大和带通滤波电路105、数据采集电路103和电源管理单元104。
线圈系2接收由地层产生的测量信号,分别通过发射驱动电路102、前置放大和带通滤波电路105将处理后的信号传输至数据采集电路103,数据采集处理电路103将该测量信号进行信号选择、自动增益、A/D采集及数据处理再进行组织编码,并通过CAN接口电路将数据传送至井下遥测。
其中,发射驱动电路102包括二级刻度电路和发射解码与驱动电路,二级刻度电路分别与线圈系2和前置放大和带通滤波电路105连接,发射解码与驱动电路分别连接线圈系2和数据采集电路103相连。
前置放大和带通滤波电路105包括若干个前置放大电路和带通滤波电路,各前置放大电路和带通滤波电路分别连接线圈系2的8通道和数据采集电路103,实现同步采集与处理,以及系统通讯。
数据采集电路103包括与发射解码与驱动电路连接的发射控制波形电路,与若干个前置放大电路和带通滤波电路相连的电平转换电路,各电平转换电路分别依次对应连接8个PGA可编程增益放大器、高速A/D转换电路;所述发射控制波形电路和高速A/D转换电路分别连接FPGA主控单元,FPGA主控单元通过CAN接口连接DSP微处理器。
DSP微处理器分别连接线圈系温度检测电路、自然电位SP、工作电压检测电路和电子仪温度检测电路。
本发明采取厚膜集成电路模块,其中,采集处理模块单元作为高性能集成电子线路1系统控制核心,完成发射波形控制、电平转换、自动增益控制与校准、辅助参数测量、8通道同步采集与处理、系统通讯等功能。基于DSP及FPGA嵌入式架构设计,如图2中数据采集框所示。采用具有32bit高速高性能DSP和FPGA作为主控元件,DSP支持1.8V内核电压、运行速度150MHz,带有0等待周期片上SRAM和128K×16bit Flash及ADC,SCI、SPI、McBSP、eCAN等丰富片上外设;采用百万门级FPGA,速度350MHz,内核电压1.5V,144Kb双端口SRAM/1Kb Flash和多达300个用户I/O。利用高速FPGA的实时性实现8通道同步采集和利用高性能DSP灵活性及处理能力实现信号处理,使得在保证仪器测量功能前提下缩小电路尺寸成为可能。
主控采集处理单元工作流程为:电平转换电路完成前端8通道调理信号的电平转换,将双极性信号调整到适于AD采集的单极性信号;可编程增益放大电路实现各通道信号的自动增益调整放大,保证仪器信号测量动态范围;FPGA控制8道AD完成信号采集,累加降噪处理,将结果保存在内部设置双口RAM,通过中断通知DSP读取数据,完成8通道3种频率信号的数字相敏检波DSPD处理,形成24组实部和44组虚部数据;一个CAN通信周期内,DSP还要通过内置AD完成电子仪及线圈系温度、自然电位SP、工作电压等辅助参数的测量,DSP数据成帧后通过内置eCAN经由遥传完成数据上传。
其中,主控采集单元DSP&FPGA交互工作时序,如图3所示,DSP收到遥传命令启动FPGA开始采集,FPGA控制AD依次完成接收线圈主测量信号预采样、PGA增益设置、主测量信号采集、二级刻度信号采集、PGA*1档和PGA*m档信号的采集,PGA每完成一种信号的采集与累加处理后均中断通知DSP取走,并由DSP完成数据处理与保存,所有8个通道由FPGA控制同步完成对应信号采集过程,由DSP分时完成处理,最后数据成帧,等待接收数据帧命令后上传数据。
厚膜集成电路模块,主要包括发射电路和前放带通电路,如图4所示,前放带通电路,主要由测量与刻度切换电路、前置放大电路、带通滤波电路以及求和输出电路共四部分构成。测量与刻度切换电路主要用于对测量信号与刻度信号进行切换。前置放大电路主要用于对微弱信号进行放大,为了有效抵制噪声前放带通板采取差分输入,同时为了消除漂移电压,前置放大采用两级差分输入放大结构。前级采用深度电压串联负反馈电路,后级通过减法器实现差分输入。该电路具有高输入阻抗,很强的共模抑制能力和较小的输出漂移电压。发射电路,阵列感应发射信号为三电平方波信号,发射电流大小分别为:26.256kHz为1A,52.512kHz为1/4A,105.024kHz为1/16A,发射信号接收数据采集短节送来的时钟逻辑控制。为了实时刻度,发射模块设计了一个电流取样电路,其取样信号送到接收短节作为二级刻度信号源。发射模块由发射波形产生器、开关控制产生器、高压滤波和电流取样电路等部分组成。
所述厚膜集成电路单元模块,包括发射电路单元;前放带通电路单元。其中发射电路单元,设计了三种发射频率,远阵列使用低频信号,近阵列使用高频信号,中阵列使用中频信号。远阵列接收线圈离发射线圈最远,在相同的发射电磁强度情况下接收的感应信号最小,选取发射电磁强度最大,使信号幅度得到提高。相反,近阵列接收线圈离发射线圈最近,在相同的发射电磁强度情况下接收的感应信号最大,此时选取发射电磁强度最小,使接收信号幅度基本一致。这种等时交替变化的发射信号的幅度与频率的平方成反比。前放带通电路单元,由测量与刻度切换电路、前置放大电路、带通滤波电路以及求和输出电路共四部分构成,主要完成短阵列和长阵列共8道测量信号与二级刻度信号的前置放大和滤波处理。
主控采集单元产生同步发射控制波形,发射单元通过发射电路解码,将混频信号驱动后送至发射线圈T进行发射,8组接收线圈阵列接收经过地层后的二次感应信号,由前置放大电路进行低噪声放大、带通电路进行选频滤波,送主控单元完成采集、检波计算、数据成帧后通过遥测上传至地面进行合成聚焦处理。
本发明实施例,提供一种线圈系结构模块单元,如图5所示,所述仪器线圈系结构单元模块由线圈系芯轴和线圈骨架与线圈屏蔽层组成,其中,线圈系芯轴采用新型管材为主体、外表面覆盖玻璃钢的复合芯轴,其中芯轴是无磁铍青铜金属管长轴,外壁上有玻璃钢套,将金属管长轴与玻璃钢套复合粘接为一体构成阵列线圈系骨架,并且采用一种组合线圈的方式,将各个阵列的接收线圈和屏蔽线圈与安装的阵列线圈系一一对应,通过特殊三线圈系结构设计,线圈系由一个发射线圈、8组接收线圈阵列组成,其中每组接收阵列由主接收线圈和屏蔽线圈组成,与发射线圈一起组成三线圈系结构。为缩短线圈系长度,而且保证径向探测深度与纵向分辨率,8组接收阵列主线圈R1至R8源距分别为6、9、12、15、21、27、39、72inch,变为单边结构8组接收线圈中间6组接收线圈R1至R6与屏蔽线圈B2至B7线圈在物理位置上存在硬件干涉,通过一体化线圈系结构设计和采用特殊的绕线工艺,将原来11组线圈(发射T、8组主接收R1-R8、8组屏蔽线圈B1-B8),变为11组组合线圈(发射T、B1、R1B2、R2B3、R3B4、R4B5、R5B6、R6B7、R7、B8、R8)。单、双边结构线圈系布局图见图5,线圈系净长由2.5m缩短为1.8m,同时通过绕线工艺改进与布线实验优化,解决了单边结构线圈各道直接耦合信号过大技术难题。
如图6所示,所述压力平衡装置3包括皮囊301、左堵头302、皮囊支撑轴303、皮囊芯体304、O型密封圈305、皮囊右堵头306、金属扎带307。
在皮囊支撑轴303外侧包裹有皮囊301,皮囊支撑轴303的一端皮囊301内镶装有皮囊左堵头302,皮囊左堵头302与皮囊支撑轴303之间设有皮囊芯体304;在皮囊支撑轴303的另一端镶装有套有O型密封圈305的皮囊右堵头306,皮囊右堵头306通过金属扎带307连接球形堵头,皮囊左堵头302通过金属扎带307连接线圈系2下管轴。
安装本压力平衡装置3时,首先将皮囊左堵头302右端压入皮囊301内,然后将301和302这两个零件整体从皮囊支撑轴下端套入到皮囊支撑轴303上,再将皮囊芯体304从皮囊支撑轴303下端装入皮囊301内,接着把装有O型密封圈305的皮囊右堵头306从皮囊支撑轴下端装入,让皮囊右堵头306左端压入皮囊301内,最后用金属扎带307将皮囊两端与左右堵头连接处可靠扎紧,完成整个压力平衡装置的安装。
在测井过程中,井下温度与压力不断变化,快测阵列感应仪器内部硅油的压力也随之变化,通过皮囊301膨胀,使得仪器内外压力在小压差下保持平衡。当温度升高时,测井仪内部的硅油体积膨胀,内部压力大于外部压力,皮囊膨胀;当温度降低时,测井仪内部的硅油体积缩小,内部压力小于外部压力,皮囊收缩。在测井过程中,井下压力与温度共同作用在仪器上,通过皮囊膨胀来适应仪器内部硅油体积变化,从而保持仪器内外压力平衡。
本发明由高性能集成电子线路装置1、线圈系2、压力平衡装置3组成的井下测量装置总长度为4.96m,其中:高性能集成电子线路装置1采取电路模块分区布置,大电流布线通过屏蔽管,微弱信号放大模块一体屏蔽等措施。并且将所有组件、电路部分的设计都采用模块化设计,组装和维修简便,长度为1.5米、线圈系2,采取发射线圈上移,八个接收线圈系依次向下排列,并且每个线圈系的补偿线圈和上个线圈系的接收线圈采用了同槽复线绕制设计。该设计在解决了线圈位置复用的同时解决了线圈间在高温高压情况下产生形变所导致的相互影响,长度为2.66米、压力平衡装置3采取皮囊式压力平衡短节,缩短了仪器长度的同时,延长了维护的周期,长度为0.8米。
本发明实施例提供一种快测阵列感应数据处理方法,如图7所示,主要包括下述步骤:原始测井数据401、测井预处理501、趋肤效应校正502、井眼校正503、真分辨率聚焦504、垂直分辨率505、和径向反演507。
步骤401,通过地面采集模块(4),完成实时采集快测阵列感应仪器测井数据,并同时将测井数据送入预处理单元(501);
步骤501,仪器在井下测量实时采集数据,经过地面采集模块4,得到八组线圈系的24条实电压信号,送入预处理单元501,调用刻度系数将电压信号转化为转化成实时测量的地层电导率信号。
步骤502,将预处理单元501输出的地层电导率信号输入趋肤效应校正单元502,利用三频趋肤效应校正处理算法,对24条实地层电导率信号进行校正得到8条测井数据曲线SEC1~SEC8。
上述三频趋肤效应校正处理算法利用下述表达式实现三频趋肤效应校正其中wi为第i频率的权重系数,σai为第i频率的视电导率值,E0,E1和E2为选择函数。在实际测井中,视电导率与频率的关系是3个离散的点,不是简单的函数,因而必须解决离散数据的求导问题,本发明采用多项式拟合法求解导数,σaR(f)=x0+x1f+x2f2其中xi是待求的多项式系数,i=0,1,3是多项式的次数。已知频率fi对应的视电导率为σaR(fi)。
步骤503,将趋肤效应校正单元(502)得到8条数据曲线,输入到下一级井眼校正处理单元(503),井眼校正利用σbhc=σc-Δσ表达式进行校正处理,其中σc为仪器测量响应值、Δσ为井眼影响数值,由表达式Δσ=σbhmbh可以得到。其中σbh为σbh(hd,σmt,decc)数值计算的井眼模型响应值,即井眼校正库是依据快侧阵列感应仪器结构,利用三维正演模拟程序计算建立仪器居中和偏心情况下的井眼校正库、σmbh为数值计算的无井眼模型响应值。本实施例中建立了8个子阵列在18个井径hd、22个地层电导率σm、22个泥浆电导率σm、10个偏心decc距的大型居中和偏心井眼校正库。通过采用查表的形式调用井眼校正库在辅助测量参数泥浆电阻率σm、井径hd、偏心decc给定的基础上,给定计算精度要求,计算仪器测量响应与给定参数模型预测响应的最小二乘方数值,得到井眼影响数值Δσ。利用σbhc=σc-Δσ表达式进行校正处理,完成对实际测井资料的仪器居中和偏心校正处理功能,最终输出8条去除井眼影响的数据曲线BHC1~BHC8;
步骤504,将井眼校正处理单元(503)得到8条去除井眼影响的数据曲线,输入到下一级真分辨率聚焦单元,利用聚焦处理算法单元:
其中,σpk,z)为合成后径向探测深度ρk,纵向深度z处的电导率;k=1,2,…,K,K是合成探测深度数目这里表示为T10、T20、T30、T60、T90;为第j个子阵列在z点的测量值,j=1,2,…,J,J是子阵列个数;是聚焦合成到探测深度ρk时,第j个子阵列的滤波器,σ是地层的电导率;z=zmin~zmax是聚焦合成的窗口范围,它是聚焦合成滤波器是地层电导率的函数。
在快侧阵列感应滤波器设计中,将8个子阵列3种频率的24条几何因子响应曲线经过三频趋肤效应校正,得到校正后的8条几何因子曲线利用最优化理论得到不同背景电导率和不同探测深度下的各个子阵列的滤波器本实施例中建立了8个子阵列在11个背景电导率σ(0.001、0.01、0.02、0.05、0.1、0.2、0.4、0.6、0.8、1.0和2.0S/m)的真分辨率聚焦滤波器库。通过调用事先算好的聚焦滤波器库进行真分辨率聚焦处理,得到5条不同探测深度的电阻率曲线T10、T20、T30、T60、T90;
步骤505,根据真分辨率聚焦单元(504)的电阻率曲线T10、T20、T30、T60、T90数据,输入到下一级垂向分辨率匹配处理单元(505),首先利用要设计8个子阵列3种频率的24条几何因子响应曲线经过三频趋肤效应校正,得到校正后的8条几何因子曲线,调用事先算好的聚焦滤波器库在相邻两条曲线中高分辨率曲线变换为低分辨率曲线的滤波器,即设计匹配滤波器hvgfi,使gpvgfi变换为gpvgf(i+1),用数学公式表示为:
gpvgf(i+1)(z)=hvgfi(z)*gpvgfi(z)
式中,i=1,2,…,4,gpvgf(i+1)(z)为将要分辨率匹配的目的函数。
本实施例中建立了8个子阵列在11个背景电导率σ(0.001、0.01、0.02、0.05、0.1、0.2、0.4、0.6、0.8、1.0和2.0S/m)的分辨率匹配滤波器库
调用分辨率匹配滤波器hvgfi进行处理得到3组分辨率为0.3m、0.6m和1.2m的曲线,每一组曲线有5条探测深度为0.25m、0.50m、0.75m、1.50m和2.25m的曲线,输出合成处理后的3种分辨率5种探测深度的15条曲线处理,分别命名为AO10~AO90、AT10~AT90、AF10~AF90;
步骤506,最后将分辨率匹配处理后的分辨率为0.6m的曲线,探测深度为0.25m、0.50m、0.75m、1.50m和2.25m的曲线通过径向反演单元(506)进行处理,该反演处理单元采用常规通用的基于阵列感应仪器的3参数反演算法,处理后输出得到真地层电阻率Rt和侵入带电阻率Rxo以及侵入带半径r,即完成快测阵列感应测井过程。
本发明实施例提供的快测阵列感应数据处理单元的全部或部分内容可以通过软件编程实现,其软件程序存储在可读取的存储介质中,存储介质例如:计算机中的硬盘、光盘或软盘。通过本发明的一种快测阵列感应测井系统,实现了从仪器到地面资料采集4到数据处理单元5进行处理,完成处理解释成果曲线输出到记录显示单元6,其中显示单元6具备生成测井资料成果打印和导出、生成图片的功能。
本发明实施例提供一种快测阵列感应测井仪实测资料处理结果,如图8所示,在X井综合解释对比图,油层段响应表现为电阻率的升高,水层段表现为低阻负差异,干层及非渗透层曲线关系重合,体现了快测阵列感应测井仪对储层渗透性及含油水的识别能力。
以上实施例提供的技术方案中的全部或部分内容可以通过工程化机械部分、硬件电路和软件编程实现,其硬件包括实现仪器线圈系、电子仪、压力平衡;软件包括采集处理单元、数据处理单元,程序存储在可读取的存储介质中,存储介质例如:采集处理DSP、FPGA内存和计算机中的硬盘、光盘。以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种快测阵列感应测井系统,其特征在于,该系统包括:
高性能集成电子线路装置,包括厚膜集成电路单元模块和高性能采集处理单元模块,用于产生三电平控制信号提供给发射驱动电路,启动发射驱动电路;
线圈系,包括线圈系芯轴和线圈骨架与线圈屏蔽层,用于支撑阵列感应仪器线圈系总成,线圈系芯轴采用新型管材为主体、外表面覆盖玻璃钢的复合芯轴;
压力平衡装置,采取皮囊压力平衡的方式与复合芯轴的结构相匹配设计,采用皮囊注油的方式,在井下高温高压环境下保持线圈系玻璃钢内外压力的平衡;
地面采集模块,用于实时采集、监测快测阵列感应仪器测井数据,并完成控制仪器刻度文件操作、车间刻度操作和温度传感器标定操作;
数据处理模块,用于对地面采集模块的原始测井数据进行处理,得到最终用户需求解释的地层电阻率测井曲线;
记录仪显示单元,用于记录显示地面采集模块和数据处理模块的地层信息曲线和图像数据,具有生成测井资料成果打印和导出、生成图片的功能;
经由高性能集成电子线路装置、线圈系、压力平衡装置组成的井下测量装置,将测井得到快测阵列感应测量信号输送至地面采集模块,经过数据处理模块进行数据合成处理后通过记录仪显示单元显示,从而完成阵列感应测井系统的信号测量;
所述线圈系与高性能集成电子线路装置相连,高性能集成电子线路装置的厚膜集成电路单元模块包括发射驱动电路、前置放大和带通滤波电路;高性能采集处理单元模块包括数据采集电路,高性能集成电子线路装置还包括电源管理单元;
线圈系接收由地层产生的测量信号,分别通过发射驱动电路、前置放大和带通滤波电路将处理后的信号传输至数据采集电路,数据采集电路将该测量信号进行信号选择、自动增益、A/D采集及数据处理再进行组织编码,并通过CAN接口电路将数据传送至井下遥测;
所述发射驱动电路包括二级刻度电路和发射解码与驱动电路,二级刻度电路分别与线圈系和前置放大和带通滤波电路连接,发射解码与驱动电路分别连接线圈系和数据采集电路。
2.根据权利要求1所述的系统,其特征在于,所述前置放大和带通滤波电路包括若干个前置放大电路和若干个带通滤波电路,各前置放大和带通滤波电路分别连接线圈系的8通道和数据采集电路,实现同步采集与处理,以及系统通讯。
3.根据权利要求1所述的系统,其特征在于,所述数据采集电路包括与发射解码与驱动电路连接的发射控制波形电路,与若干个前置放大和带通滤波电路相连的电平转换电路,各电平转换电路分别依次对应连接8个PGA可编程增益放大器、高速A/D转换电路;所述发射控制波形电路和高速A/D转换电路分别连接FPGA主控单元,FPGA主控单元通过CAN接口连接DSP微处理器;
DSP微处理器分别连接线圈系温度检测电路、自然电位SP、工作电压检测电路和电子仪温度检测电路。
4.根据权利要求1所述的系统,其特征在于,所述压力平衡装置包括皮囊支撑轴,及包裹在皮囊支撑轴外侧的皮囊,所述皮囊支撑轴的一端皮囊内镶装有皮囊左堵头,皮囊左堵头与皮囊支撑轴之间设有皮囊芯体;在皮囊支撑轴的另一端镶装有套有O型密封圈的皮囊右堵头,皮囊右堵头通过金属扎带连接球形堵头,皮囊左堵头通过金属扎带连接线圈系下管轴。
5.一种基于权利要求1所述的快测阵列感应测井系统的测井方法,其特征在于,该方法包括下述步骤:
步骤401,通过地面采集模块,完成实时采集快测阵列感应仪器测井数据,并同时将测井数据送入预处理单元;
步骤501,仪器在井下测量实时采集数据,经过地面采集模块,得到八组线圈系的24条实电压信号,送入预处理单元,调用刻度系数将电压信号转化为转化成实时测量的地层电导率信号;
步骤502,将预处理单元输出的地层电导率信号输入趋肤效应校正单元,通过三频趋肤效应校正处理算法,对24条实地层电导率信号进行校正得到8条测井数据曲线SEC1~SEC8;
步骤503,将趋肤效应校正单元得到8条数据曲线,输入到下一级井眼校正处理单元,井眼校正处理,依据快测阵列感应仪器结构,利用三维正演模拟程序计算建立了大型居中和偏心井眼校正库;通过采用查表的形式调用井眼校正库,对趋肤效应校正单元的8条曲线做井眼校正处理,完成对实际测井资料的仪器居中和偏心校正处理功能,最终输出8条去除井眼影响的数据曲线BHC1~BHC8;
步骤504,将井眼校正处理单元得到8条去除井眼影响的数据曲线,输入到下一级真分辨率聚焦单元,利用聚焦处理算法单元对8条信息数据进行井眼校正,通过调用事先算好的聚焦滤波器库进行真分辨率聚焦处理,得到5条不同探测深度的电阻率曲线T10、T20、T30、T60、T90;
步骤505,根据真分辨率聚焦单元的电阻率曲线T10、T20、T30、T60、T90数据,输入到下一级垂向分辨率匹配处理单元,利用8个子阵列3种频率的24条几何因子响应曲线经过三频趋肤效应校正,得到校正后的8条几何因子曲线,调用事先算好的聚焦滤波器库将相邻两条曲线中高分辨率曲线变换为低分辨率曲线,即设计分辨率匹配滤波器,调用分辨率匹配滤波器进行处理得到3组分辨率为0.3m、0.6m和1.2m的曲线,每一组曲线有5条探测深度为0.25m、0.50m、0.75m、1.50m和2.25m的曲线,输出合成处理后的3种分辨率5种探测深度的15条曲线处理,分别命名为AO10~AO90、AT10~AT90、AF10~AF90;
步骤506,最后将分辨率匹配处理后的分辨率为0.6m的曲线,探测深度为0.25m、0.50m、0.75m、1.50m和2.25m的曲线通过径向反演单元进行处理,该反演处理单元采用常规通用的基于阵列感应仪器的3参数反演算法,处理后输出得到真地层电阻率Rt和侵入带电阻率Rxo以及侵入带半径r,即完成快测阵列感应测井过程。
6.根据权利要求5所述的快测阵列感应测井系统的测井方法,其特征在于,所述步骤502中,利用三频趋肤效应校正处理算法,利用下述表达式实现三频趋肤效应校正:
σ sec = E 0 Σ i = 1 3 w i σ a i + E 1 Σ i = 1 3 w i ( σ a i - 2 f i ∂ σ a i ∂ f i ) + E 2 Σ i = 1 3 w i ( σ a i - 4 f i 3 ∂ σ a i ∂ f i + 4 f i 2 3 ∂ 2 σ a i ∂ f i 2 ) ,
其中wi为第i频率的权重系数,σai为第i频率的视电导率值,E0,E1和E2为选择函数,fi为频率。
7.根据权利要求5所述的快测阵列感应测井系统的测井方法,其特征在于,所述步骤504中,利用聚焦处理算法单元对8条信息数据进行井眼校正:
σ p ( ρ k , z ) = Σ j = 1 J Σ z ′ = z min z max w j k ( z ′ , σ ) σ a j ( z - z ′ ) ,
其中,σpk,z)为合成后径向探测深度ρk,纵向深度z处的电导率;k=1,2,…,K,K是合成探测深度数目,这里表示为T10、T20、T30、T60、T90;为第j个子阵列在z点的测量值;j=1,2,…,J,J是子阵列个数;是聚焦合成到探测深度ρk时,第j个子阵列的滤波器;σ是地层的电导率;z=zmin~zmax是聚焦合成的窗口范围,它是聚焦合成滤波器的地层电导率的函数。
8.根据权利要求5所述的快测阵列感应测井系统的测井方法,其特征在于,所述步骤505中,调用事先算好的聚焦滤波器库聚焦滤波器库通过下述方式实现:
设计匹配滤波器hvgfi,使gpvgfi变换为gpvgf(i+1),为:
gpvgf(i+1)(z)=hvgfi(z)*gpvgfi(z)
式中,i=1,2,…,4,gpvgf(i+1)(z)为将要分辨率匹配的目的函数。
CN201410714371.6A 2014-11-28 2014-11-28 一种快测阵列感应测井系统及其测井方法 Active CN104481525B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410714371.6A CN104481525B (zh) 2014-11-28 2014-11-28 一种快测阵列感应测井系统及其测井方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410714371.6A CN104481525B (zh) 2014-11-28 2014-11-28 一种快测阵列感应测井系统及其测井方法

Publications (2)

Publication Number Publication Date
CN104481525A CN104481525A (zh) 2015-04-01
CN104481525B true CN104481525B (zh) 2017-06-06

Family

ID=52756117

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410714371.6A Active CN104481525B (zh) 2014-11-28 2014-11-28 一种快测阵列感应测井系统及其测井方法

Country Status (1)

Country Link
CN (1) CN104481525B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105572761B (zh) * 2015-12-14 2017-11-07 中国石油天然气集团公司 一种感应成像测井技术校正与标定系统
CN106525318B (zh) * 2016-12-16 2022-07-01 山东金洲科瑞节能科技有限公司 一种总线型压力测量装置
CN106837299B (zh) * 2017-01-19 2019-11-19 中国科学院声学研究所 一种井眼校正的系统及方法
CN107463722B (zh) * 2017-05-05 2020-05-05 西南石油大学 一种薄层电阻率测井响应校正的新方法
CN108049865B (zh) * 2017-09-12 2021-08-03 中国石油天然气集团公司 随钻感应测井电场刻度方法
CN110965983A (zh) * 2018-09-27 2020-04-07 中国石油化工股份有限公司 一种井下电磁探伤仪
CN109695446B (zh) * 2019-01-04 2021-06-08 电子科技大学 一种感应测井仪发射功率自适应调整装置
CN109695447B (zh) * 2019-01-04 2021-06-08 电子科技大学 一种感应测井仪发射功率自适应调整方法
CN109711109B (zh) * 2019-02-22 2023-12-22 中国石油天然气集团有限公司 一种电阻率测井仪器电极系结构参数智能优化的方法及装置
CN110513105B (zh) * 2019-07-30 2023-04-25 中国石油天然气集团有限公司 一种三维阵列感应测井仪器数据处理系统及方法
CN111255436B (zh) * 2020-01-17 2023-07-28 成都理工大学 一种改善hdil阵列感应测井仪测量曲线质量的方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5698982A (en) * 1996-03-18 1997-12-16 Computalog Research, Inc. Method and system for skin effect correction in a multiple transmit frequency induction logging system
CN2723677Y (zh) * 2004-08-20 2005-09-07 中国石油天然气集团公司 阵列感应成像测井仪
CN201288549Y (zh) * 2008-11-17 2009-08-12 陕西巨丰思源科技有限公司 数字聚焦阵列感应测井仪
CN102042009A (zh) * 2009-10-13 2011-05-04 中国石油天然气集团公司 一种测量地层电阻率的阵列感应测井线圈系
CN202510111U (zh) * 2012-01-16 2012-10-31 中国海洋石油总公司 应用于阵列感应测井仪的多通道数据采集电路
CN103643950A (zh) * 2013-12-23 2014-03-19 北京环鼎科技有限责任公司 高温高压小井眼阵列感应测井仪

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6476609B1 (en) * 1999-01-28 2002-11-05 Dresser Industries, Inc. Electromagnetic wave resistivity tool having a tilted antenna for geosteering within a desired payzone

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5698982A (en) * 1996-03-18 1997-12-16 Computalog Research, Inc. Method and system for skin effect correction in a multiple transmit frequency induction logging system
CN2723677Y (zh) * 2004-08-20 2005-09-07 中国石油天然气集团公司 阵列感应成像测井仪
CN201288549Y (zh) * 2008-11-17 2009-08-12 陕西巨丰思源科技有限公司 数字聚焦阵列感应测井仪
CN102042009A (zh) * 2009-10-13 2011-05-04 中国石油天然气集团公司 一种测量地层电阻率的阵列感应测井线圈系
CN202510111U (zh) * 2012-01-16 2012-10-31 中国海洋石油总公司 应用于阵列感应测井仪的多通道数据采集电路
CN103643950A (zh) * 2013-12-23 2014-03-19 北京环鼎科技有限责任公司 高温高压小井眼阵列感应测井仪

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈涛.阵列感应成像测井仪研究与设计.《西安石油大学硕士学位论文》.2008, *

Also Published As

Publication number Publication date
CN104481525A (zh) 2015-04-01

Similar Documents

Publication Publication Date Title
CN104481525B (zh) 一种快测阵列感应测井系统及其测井方法
CN102042009B (zh) 一种测量地层电阻率的阵列感应测井线圈系
CN109632575B (zh) 一种监测河岸带潜流交换速率的装置及其监测方法
CN106291719B (zh) 一种阵列人工源磁场频率测深方法
CN107630697A (zh) 基于随钻电磁波电阻率测井的地层电阻率联合反演方法
CN104903749A (zh) 具有多极传感器的深方位系统
CN103590823A (zh) 测量钻头处地层电阻率的装置
CN104865614A (zh) 基于可变骨架参数的复杂储层流体识别方法
CN116665067B (zh) 基于图神经网络的找矿靶区优选系统及其方法
CN106646661A (zh) 一种矿床水文地质综合勘查系统
CN108810838A (zh) 基于智能手机室内背景声感知的房间级定位方法
Raza et al. A survey on subsurface signal propagation
CN110400371A (zh) 一种水平构造地貌实体的三维模型构建方法
CN111927443B (zh) 一种基于电场测量地层真电阻率的测井系统与方法
CN204371322U (zh) 一种快测阵列感应测井系统
CN107632322B (zh) 一种适用于水域电法勘探的电缆系统的勘探方法
CN105891895A (zh) 一种确定天波传播特性的系统和方法
CN208999255U (zh) 一种监测河岸带潜流交换速率的装置
CN109154189A (zh) 具有减少的运动伪影的t2反演
CN203759264U (zh) 基于工频整周期触发的核磁共振双极性叠加消噪装置
CN103643950A (zh) 高温高压小井眼阵列感应测井仪
CN108363739A (zh) 一种基于稀疏采集的地震资料高低频拓展方法
CN107436448A (zh) 一种工程地震探测系统及数据读取系统
CN113359199A (zh) 一种基于聚焦磁场的井周电阻率测量方法
Rodriguez et al. Groundwater resource exploration and mapping methods: a review

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant