CN105239999B - 一种快速确定多频核磁共振测井仪井下发射功率的方法 - Google Patents
一种快速确定多频核磁共振测井仪井下发射功率的方法 Download PDFInfo
- Publication number
- CN105239999B CN105239999B CN201510563752.3A CN201510563752A CN105239999B CN 105239999 B CN105239999 B CN 105239999B CN 201510563752 A CN201510563752 A CN 201510563752A CN 105239999 B CN105239999 B CN 105239999B
- Authority
- CN
- China
- Prior art keywords
- magnetic field
- transmission power
- field intensity
- pulsed magnetic
- cal
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种快速确定多频核磁共振测井仪井下发射功率的方法,包括(1)确定每个工作频率对应刻度发射功率下的脉冲磁场强度最大值B1max;(2)按照各个频率刻度发射功率值GAcal第一次测试;采集各频率发射的脉冲磁场强度实测值B1和地层温度T;(3)计算所脉冲磁场强度实测值B1mod;(4)判断实时采集的脉冲磁场强度实测值B1mod和车间刻度脉冲磁场强度最大值B1max之间误差是否满足要求;(5)根据地层实时采集的脉冲磁场强度实测值B1mod,计算各个频率的功率值GAnext;直到满足要求。本发明解决了现有测井仪器由手工确定井下发射功率所导致时间长,不准确,测井质量低等问题。地层信息实时反馈,实现快速确定发射功率。该方法利于实现,自动进行准确寻找。
Description
技术领域
本发明涉及石油勘探测井技术领域,涉及一种快速确定多频核磁共振测井仪井下发射功率的方法。
背景技术
核磁共振测井是利用氢核在已知磁场中的核磁共振现象,探测地层孔隙及所含流体特征的一种方法。核磁共振测井测量的是氢核发生核磁共振后自由进动过程的衰减时间和振幅。振幅信息与探测区氢核的数量成正比,通过刻度可以获得地层孔隙度信息,它不受放射源和岩性的影响。
多频核磁共振测井仪是我国研制的具有自主知识产权居中测量型核磁共振测井仪。它主要由储能线路、电子线路、扶正器和探头四部分组成,其中探头由磁体和天线组成。该仪器采用5个频带共9个不同的频率进行测量,由永久磁铁产生均匀的静磁场,使地层中的氢原子核产生极化。仪器在井内居中测量,在井眼周围地层中形成以井轴为中心,直径为14-16″,厚度1mm,高24″,彼此之间相距1mm的9个圆柱壳。
仪器采用CPMG脉冲进行横向弛豫时间的测量,由仪器天线用一定功率发射对应的射频脉冲,使地层中氢核旋转所形成的磁化矢量分别扳倒90°和180°,形成CPMG脉冲序列,同时由天线接收回波弛豫信号。核磁共振测井仪器发射功率由于受到井眼流体的影响,以及在不同深度地层和地层水的电特性影响和温度的影响,仪器的发射功率会随着变化。这些变化会导致发射CPMG脉冲时不能够对地层中氢原子核磁化矢量实现扳转90°和180°,导致采集的回波弛豫信号信噪比低,不能在复杂油气藏中有效识别油气。目前核磁共振仪器一般采用手工调整的办法在确定发射功率,由于频率多,每个频率查找一次,导致寻找时间长,发射功率也不是最佳功率。因此有必要采用一种能够根据地层反馈快速确定核磁共振测井仪井下发射功率的方法来提高测井信号质量。
发明内容
本发明的目的在于针对现有多频核磁共振测井仪器由手工确定井下发射功率所导致测井时间长,查找功率不准确,间接增加了测井风险,降低了测井质量等问题。为此,本发明提出了一种快速确定核磁共振测井仪井下发射功率的办法,利于实现,可以在测井过程中自动进行准确的寻找。
本发明的目的是通过下述技术方案来实现的。
提供一种基于实时地层信息反馈,根据线性逼近原则计算功率,通过不断逼近,实现快速确定发射功率。
一种快速确定多频核磁共振测井仪井下发射功率的方法,包括下述步骤:
(1)在多频核磁共振测井仪器下放井下过程中,根据测井仪器测量周围地层的增益,查找对应的水箱刻度文件,确定每个工作频率对应刻度发射功率下的脉冲磁场强度最大值B1max;
(2)测井仪器按照每个频率刻度发射功率值GAcal发射射频磁场到地层,按照各个频率刻度发射功率值GAcal进行第一次测试;同时测井仪器采集各个不同地层段频率发射的脉冲磁场强度实测值B1和地层温度T;
(3)根据不同地层段频率发射的脉冲磁场强度实测值B1和地层温度T,计算所在地层段的温度情况下的脉冲磁场强度实测值B1mod;
(4)判断该段地层实时采集的脉冲磁场强度实测值B1mod和车间刻度脉冲磁场强度最大值B1max之间误差是否满足要求;如果满足要求,表明该段地层情况下使用该频率使用该功率值GA满足要求,否则,进行下一步;
(5)根据地层实时采集的脉冲磁场强度实测值B1mod,计算各个频率的功率值GAnext;然后重复步骤(2)-(4),直到满足要求。
进一步地,所述步骤(1)确定脉冲磁场强度最大值B1max,步骤如下:
(1a)通过车间刻度时测量的各个频率的脉冲磁场强度车间刻度值B’1和回波幅度Amp之间的关系,按照二次曲线拟合出关系式:
式中,a为二次项拟合系数,b为一次项拟合系数,c为常量;
(1b)确定脉冲磁场强度最大值B1max
B1max=(-b)/2a ②。
进一步地,所述步骤(2),按照刻度时的发射功率进行第一次测试过程如下:
2a)根据车间刻度射频磁场强度B1cal和频率刻度发射功率GAcal关系,拟合出关系式:
GAcal=Kcal×B1cal+C ③
式中,Kcal为线性拟合系数;C为线性拟合的常量;
2b)根据车间刻度脉冲磁场强度最大值B1max,得到各个频率的第一次测试发射功率GA’cal
GA’cal=Kcal×B1max+C; ④;
2c)以各个频率的第一次测试发射功率GA’cal发射射频磁场,同时采集对应地层条件下脉冲磁场强度实测值B1和地层温度T。
进一步地,所述步骤(3)确定所在地层段的温度情况下的脉冲磁场强度实测值B1mod,步骤如下:
(3a)通过车间温度实验确定温度校正系数,利用最小二乘法对地层温度T和脉冲磁场强度实测值B1曲线进行线性拟合,计算出对应的校正系数Kt;
(3b)在地层温度T的情况下,设Tc是刻度时的室内温度,则
B1mod=(1+Kt×(T-Tc)×B1 ⑤。
进一步地,所述步骤(5),计算各个频率的功率值GAnext,步骤如下:
(5a)根据步骤(1)(2)(3),按照车间刻度时发射功率GAcal发射射频磁场后,采集所在地层段的温度情况下的磁场强度B1mod;
(5b)再根据⑤式得到所在地层段的温度情况下的磁场强度B1mod,计算对应刻度条件下的功率值
GA=B1mod×Kcal+C ⑥;
(5c)计算出下一次的发射功率GAnext:
GAnext=GAcal+(GAcal-GA) ⑦;
经过前两次尝试后,再次测试则采用步骤如下:
(5d)设定前两次尝试所用的发射功率分别为GAn和GAn-1对应地层条件下的磁场强度分别为B1mod(n)和B1mod(n-1);
其中,n表示测试次数,n须≥2;
(5e)计算出发射功率值GAnext
(5f)按照计算出发射功率值GAnext发射射频磁场,同时采集在地层条件下的磁场强度,如果磁场强度与车间刻度脉冲磁场强度最大值B1max之间的误差是否满足5%的要求就设定当前测试的GAnext为最佳功率,否则继续步骤(5d)-(5f)。
本发明与现有技术相比,具有如下优点:
1.通过对刻度时的发射功率GA和脉冲磁场强度实测值B1进行拟合,加快了整个查找过程,经过大量测试一般情况下仅需查找3次;
2.该方法可以用软件实现自动查找,同时实现9个频率并行工作,查找的发射功率GA准确;
3.通过快速确定发射功率方法,能够减少人工查找的出错率,解决了查找时间长的问题,提高了测井施工的安全性和时效性,同时准确的发射功率也保证了测量回波数据质量。
附图说明
图1是多频核磁共振测井仪测井过程中应用本发明快速确定工作功率的流程图;
图2是本发明井下发射功率计算流程图;
图3是本发明引用的多频核磁共振测井仪车间刻度图;
图4是本发明引用的多频核磁共振测井仪温度和射频磁场强度补偿关系图;
图5是本发明引用的多频核磁共振测井仪射频磁场强度和发射功率关系图。
具体实施方式
以下参照附图对本发明作进一步详细描述。
参照图1所示,本发明的具体步骤如下:
步骤一,在多频核磁共振测井仪器下放井下过程中,根据测井仪器测量周围地层的增益,查找对应的水箱车间刻度数据,确定各工作频率对应刻度发射功率下的90度扳倒角所需脉冲磁场强度最大值B1max。
加载车间刻度文件,读取各个频率刻度功率值GA和车间刻度时测量的各个频率的脉冲磁场强度B’1等参数,通过车间刻度时测量的各个频率的B1和回波幅度Amp之间的关系,按照二次曲线拟合出关系式:
由公式B1max=(-b)/2a ②
确定对应的射频磁场强度最大值B1max。如图3所示。式中,a为二次项拟合系数,b为一次项拟合系数,c为常量。
步骤二,测井仪器按照刻度时的发射功率进行第一次测试。
首先,根据车间刻度射频磁场强度和发射功率关系,拟合出关系式
GAcal=Kcal×B1cal+C ③
式中,Kcal为线性拟合系数;C为线性拟合的常量;
如图5所示,以及步骤一中确定的B1max,代入计算可以得到第一次测试的发射功率GA’cal
GA’cal=Kcal×B1max+C; ④;
其次,以各个频率的第一次测试发射功率GA’cal发射射频磁场,同时采集对应地层条件下脉冲磁场强度实测值B1和地层温度T。
步骤三,计算对应地层条件下的B1mod。
根据车间刻度射频磁场强度和温度补偿关系,见图4,利用最小二乘法对地层温度T和脉冲磁场强度实测值B1曲线进行拟合,计算出对应的校正系数Kt。
在地层温度T的情况下,Tc是刻度时的室内温度,可以确定
B1mod=(1+Kt×(T-Tc)×B1 ⑤。
步骤四,比较脉冲磁场强度实测值B1mod和车间刻度脉冲磁场强度最大值刻度B1max之间的误差是否满足5%的要求,如果满足要求,表明该段地层情况下使用该频率刻度功率值GA满足要求,否则继续步骤五。
步骤五,根据地层实时采集的脉冲磁场强度实测值B1mod,计算各个频率的功率值GAnext;然后重复步骤(2)-(4),直到满足要求。见图2。
(5a)根据步骤(1)(2)(3),设定按照车间刻度时发射功率GAcal发射射频磁场后,采集对应地层条件下所在地层段的温度情况下的的磁场强度B1mod。
(5b)再根据步骤二关系式GAcal=Kcal×B1cal+C
可以得到地层条件下B′1mod所对应刻度条件下的功率值
GA=B1mod×Kcal+C ⑥;
(5c)计算出下一次的发射功率GAnext:
GAnext=GAcal+(GAcal-GA) ⑦;
经过前两次尝试后,再次测试则采用步骤如下:
(5d)设定前两次尝试所用的发射功率分别为GAn和GAn-1,对应地层条件下的磁场强度分别为B1mod(n)和B1mod(n-1);
其中,n表示测试次数,n须≥2;
(5e)计算出发射功率值GAnext:
(5f)按照计算出发射功率值GAnext发射射频磁场,同时采集在地层条件下的磁场强度,如果磁场强度与车间刻度脉冲磁场强度最大值B1max之间的误差是否满足5%的要求就设定当前测试的GAnext为最佳功率,否则继续步骤(5d)到(5f)。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种快速确定多频核磁共振测井仪井下发射功率的方法,其特征在于,该方法包括下述步骤:
(1)在多频核磁共振测井仪器下放井下过程中,根据测井仪器测量周围地层的增益,查找对应的水箱刻度文件,确定每个工作频率对应刻度发射功率下的车间刻度脉冲磁场强度最大值B1max;
(2)测井仪器按照每个频率刻度发射功率值GAcal发射射频磁场到对应的地层段,按照各个频率刻度发射功率值GAcal进行第一次测试;同时测井仪器采集各个不同地层段对应的频率发射的脉冲磁场强度实测值B1和地层温度T;
(3)根据不同地层段对应的频率发射的脉冲磁场强度实测值B1和地层温度T,计算所在地层段的温度情况下的脉冲磁场强度实测值B1mod;
(4)判断所在地层段的温度情况下的脉冲磁场强度实测值B1mod和车间刻度脉冲磁场强度最大值B1max之间误差是否满足要求;如果满足要求,表明该地层段的温度情况下使用该频率使用该功率值GA满足要求,否则,进行下一步;
(5)根据所在地层段的温度情况下的脉冲磁场强度实测值B1mod,计算各个频率的发射功率值GAnext;
具体步骤如下:
判断是否第一次测试,如果是则:
(5a)根据步骤(1)(2)(3),按照频率刻度发射功率值GAcal发射射频磁场后,计算所在地层段的温度情况下的脉冲磁场强度实测值B1mod;
(5b)再根据车间刻度射频磁场强度B1cal和频率刻度发射功率值GAcal关系得到的线性拟合系数,再结合所在地层段的温度情况下的脉冲磁场强度实测值B1mod,计算对应刻度条件下的功率值
GA=B1mod×Kcal+C ⑥;
Kcal为线性拟合系数;C为线性拟合的常量;
(5c)计算出下一次的发射功率值GAnext:
GAnext=GAcal+(GAcal‐GA) ⑦;
如果不是则:经过前两次尝试后,再次测试则采用步骤如下:
(5d)设定前两次尝试所用的发射功率分别为GAn和GAn-1对应地层条件下的磁场强度分别为B1mod(n)和B1mod(n-1);
其中,n表示测试次数,n须≥2;
(5e)计算出发射功率值GAnext
(5f)按照计算出发射功率值GAnext发射射频磁场,同时采集在地层条件下的磁场强度,如果磁场强度与车间刻度脉冲磁场强度最大值B1max之间的误差满足5%的要求,就设定当前测试的发射功率值GAnext为最佳功率,否则继续步骤(5d)-(5f)。
2.根据权利要求1所述的快速确定多频核磁共振测井仪井下发射功率的方法,其特征在于,所述步骤(1)确定每个工作频率对应刻度发射功率下的车间刻度脉冲磁场强度最大值B1max,步骤如下:
(1a)通过车间刻度时测量的各个频率的脉冲磁场强度车间刻度值B’1和回波幅度Amp之间的关系,按照二次曲线拟合出关系式:
式中,a为二次项拟合系数,b为一次项拟合系数,c为常量;
(1b)确定车间刻度脉冲磁场强度最大值B1max
B1max=(-b)/2a ②。
3.根据权利要求1所述的快速确定多频核磁共振测井仪井下发射功率的方法,其特征在于,所述步骤(2),按照刻度时的发射功率进行第一次测试过程如下:
2a)根据车间刻度射频磁场强度B1cal和频率刻度发射功率值GAcal关系,拟合出关系式:
GAcal=Kcal×B1cal+C ③
2b)根据车间刻度脉冲磁场强度最大值B1max,得到各个频率的第一次测试发射功率GA’cal
GA’cal=Kcal×B1max+C; ④;
2c)以各个频率的第一次测试发射功率GA’cal发射射频磁场,同时采集各个不同地层段对应的频率发射的脉冲磁场强度实测值B1和地层温度T。
4.根据权利要求1所述的快速确定多频核磁共振测井仪井下发射功率的方法,其特征在于,所述步骤(3)计算所在地层段的温度情况下的脉冲磁场强度实测值B1mod,步骤如下:
(3a)通过车间温度实验确定温度校正系数,利用最小二乘法对地层温度T和各个不同地层段对应的频率发射的脉冲磁场强度实测值B1曲线进行线性拟合,计算出对应的校正系数Kt;
(3b)在地层温度T的情况下,设Tc是刻度时的室内温度,则
B1mod=(1+Kt×(T-Tc)×B1 ⑤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510563752.3A CN105239999B (zh) | 2015-09-07 | 2015-09-07 | 一种快速确定多频核磁共振测井仪井下发射功率的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510563752.3A CN105239999B (zh) | 2015-09-07 | 2015-09-07 | 一种快速确定多频核磁共振测井仪井下发射功率的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105239999A CN105239999A (zh) | 2016-01-13 |
CN105239999B true CN105239999B (zh) | 2018-09-04 |
Family
ID=55037792
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510563752.3A Active CN105239999B (zh) | 2015-09-07 | 2015-09-07 | 一种快速确定多频核磁共振测井仪井下发射功率的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105239999B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111538099B (zh) * | 2020-05-06 | 2023-05-26 | 中国石油天然气集团有限公司 | 一种核磁共振测井速度确定方法 |
CN114412456B (zh) * | 2022-01-14 | 2022-09-02 | 中国科学院地质与地球物理研究所 | 一种精细化核磁共振测井脉冲序列刻度方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5055787A (en) * | 1986-08-27 | 1991-10-08 | Schlumberger Technology Corporation | Borehole measurement of NMR characteristics of earth formations |
US5451873A (en) * | 1993-10-05 | 1995-09-19 | Schlumberger Technology Corporation | Method and apparatus for determining the in situ larmor frequency of a wellbore NMR tool to compensate for accumulation of magnetic material on the magnet housing of the tool |
CN1763563A (zh) * | 2004-10-21 | 2006-04-26 | 北京大学 | 用岩石核磁共振弛豫信号测量地层岩石物性的设备及测量方法 |
CN201060934Y (zh) * | 2007-04-29 | 2008-05-14 | 中国石油天然气集团公司 | 居中型核磁共振测井仪探头天线 |
CN102418519A (zh) * | 2011-12-31 | 2012-04-18 | 北京环鼎科技有限责任公司 | 核磁共振测井仪刻度箱 |
CN103424420A (zh) * | 2013-01-31 | 2013-12-04 | 上海理工大学 | 一种基于拟合的核磁共振信号处理方法 |
-
2015
- 2015-09-07 CN CN201510563752.3A patent/CN105239999B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5055787A (en) * | 1986-08-27 | 1991-10-08 | Schlumberger Technology Corporation | Borehole measurement of NMR characteristics of earth formations |
US5451873A (en) * | 1993-10-05 | 1995-09-19 | Schlumberger Technology Corporation | Method and apparatus for determining the in situ larmor frequency of a wellbore NMR tool to compensate for accumulation of magnetic material on the magnet housing of the tool |
CN1763563A (zh) * | 2004-10-21 | 2006-04-26 | 北京大学 | 用岩石核磁共振弛豫信号测量地层岩石物性的设备及测量方法 |
CN201060934Y (zh) * | 2007-04-29 | 2008-05-14 | 中国石油天然气集团公司 | 居中型核磁共振测井仪探头天线 |
CN102418519A (zh) * | 2011-12-31 | 2012-04-18 | 北京环鼎科技有限责任公司 | 核磁共振测井仪刻度箱 |
CN103424420A (zh) * | 2013-01-31 | 2013-12-04 | 上海理工大学 | 一种基于拟合的核磁共振信号处理方法 |
Non-Patent Citations (1)
Title |
---|
多频核磁共振测井仪实时数据采集的质量控制;陈江浩等;《测井技术》;20100630;第34卷(第3期);第279-281、283页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105239999A (zh) | 2016-01-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9625601B2 (en) | Nuclear magnetic resonance rock sample analysis method and instrument with constant gradient field | |
US6255818B1 (en) | Method and apparatus for performing magnetic resonance measurements | |
US8686724B2 (en) | System for emulating nuclear magnetic resonance well logging tool diffusion editing measurements on a bench-top nuclear magnetic resonance spectrometer for laboratory-scale rock core analysis | |
RU2401442C2 (ru) | Способ определения ошибки зонда для прибора на основе индукции или распространения с поперечными или трехосными массивами | |
US6559640B2 (en) | NMR apparatus and method utilizing pulsed static magnetic fields | |
CN102648424B (zh) | 可用于采矿和矿物勘探的方法和设备 | |
US10107930B2 (en) | Hybrid saturation recovery-inversion recovery pulse sequence for improved NMR logging of boreholes | |
US20180321169A1 (en) | Selective Saturation Pulse For Downhole Measurements | |
NO342087B1 (no) | Fremgangsmåte og apparat for NMR logging med godt signal/støyforhold | |
WO2005119300A1 (en) | Method of detecting, quantifying and correcting borehole contaminations from multi-frequency, multi-sensitive volume nmr logging data | |
CN107655922A (zh) | 一种基于核磁共振的烃源岩含水量无损检测的方法 | |
NO317464B1 (no) | Estimering av permabilitet | |
AU2017204581A1 (en) | Nuclear magnetic resonance tool calibration | |
CN105239999B (zh) | 一种快速确定多频核磁共振测井仪井下发射功率的方法 | |
CN108873082B (zh) | 一种考虑弛豫组分区间的致密岩心核磁共振孔隙度校正方法 | |
US10386530B2 (en) | NMR probe and methods of use | |
US9223048B2 (en) | System and method to detect a fluid flow without a tipping pulse | |
CN110410058B (zh) | 一种校正岩心实验结果刻度二维核磁测井的方法 | |
RU90225U1 (ru) | Устройство ядерно-магнитного каротажа | |
Rastegarnia et al. | Application of TDA technique to estimate the hydrocarbon saturation using MRIL Data: A Case study for a Southern Iranian Oilfield | |
CN108507923A (zh) | 岩心核磁共振分析仪孔隙度测量精准度检验及校正方法 | |
RU2361247C1 (ru) | Способ ядерно-магнитного каротажа и устройство ядерно-магнитного каротажа | |
CN209088048U (zh) | 用于核磁共振流体分析仪的天线 | |
CN111058828A (zh) | 基于表面核磁共振技术的流体t1弛豫谱的快速测量***与方法 | |
US11795811B2 (en) | Method and system to improve open loop systems |
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 |