CN112855118B - 随钻核磁共振探头天线配置方法、装置、设备及存储介质 - Google Patents

随钻核磁共振探头天线配置方法、装置、设备及存储介质 Download PDF

Info

Publication number
CN112855118B
CN112855118B CN202110053137.3A CN202110053137A CN112855118B CN 112855118 B CN112855118 B CN 112855118B CN 202110053137 A CN202110053137 A CN 202110053137A CN 112855118 B CN112855118 B CN 112855118B
Authority
CN
China
Prior art keywords
magnetic field
inductance
drilling
antenna
nuclear magnetic
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
CN202110053137.3A
Other languages
English (en)
Other versions
CN112855118A (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 University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
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 University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN202110053137.3A priority Critical patent/CN112855118B/zh
Publication of CN112855118A publication Critical patent/CN112855118A/zh
Application granted granted Critical
Publication of CN112855118B publication Critical patent/CN112855118B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • E21B47/12Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Data Mining & Analysis (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Computational Mathematics (AREA)
  • Geophysics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Theoretical Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Pure & Applied Mathematics (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明实施例提供一种随钻核磁共振探头天线配置方法、装置、设备及存储介质,所述方法包括:采集随钻核磁共振探头的静磁场,根据静磁场确定静磁场对应的目标射频磁场;根据目标射频磁场确定电感值最小时,随钻核磁共振探头天线的天线参数;根据天线参数构建随钻核磁共振探头天线,以使随钻核磁共振探头天线产生最优射频磁场,使随钻核磁共振探头能够根据最优射频磁场激发共振信号,并接收共振信号对应的反馈信号,基于反馈信号获取原状地层信息。本发明实施例能够使设计得到的随钻核磁共振探头天线在实际应用中产生与理论中的目标射频磁场相似的最优射频磁场,进而使随钻核磁共振测井仪可以准确地探测原状地层信息。

Description

随钻核磁共振探头天线配置方法、装置、设备及存储介质
技术领域
本发明涉及核磁共振测井技术领域,尤其涉及一种随钻核磁共振探头天线配置方法、装置、设备及存储介质。
背景技术
随着石油开发的不断进行,石油勘探技术也越来越重要。石油勘探领域中,随钻核磁共振测井作为核磁共振测井的最新发展手段,能够提供更准确的原状地层信息以及更可靠的油藏评价参数。
现有技术中,随钻核磁共振测井仪的探头天线的主流设计方法是正向设计方法,即应用电磁仿真软件对天线的结构参数进行优化仿真。
然而,上述方法依赖于设计人员的主观经验判断,导致设计得到的天线在实际应用中产生的射频磁场与理论中的目标射频磁场存在较大差异,进而造成随钻核磁共振测井仪无法准确地探测原状地层信息等。
发明内容
本发明实施例提供一种随钻核磁共振探头天线配置方法、装置、设备及存储介质,以解决正向设计方法中设计得到的天线产生的射频磁场与理论中的目标射频磁场存在较大差异,进而造成随钻核磁共振测井仪无法准确地探测原状地层信息的技术问题。
第一方面,本发明实施例提供一种随钻核磁共振探头天线配置方法,所述随钻核磁共振探头天线配置方法应用于随钻核磁共振探头天线配置装置,所述配置装置用于控制安装在随钻核磁共振探头上的随钻核磁共振探头天线工作,以使所述随钻核磁共振探头根据所述随钻核磁共振探头天线产生的射频磁场激发共振信号,所述方法包括:
采集随钻核磁共振探头的静磁场,根据所述静磁场确定所述静磁场对应的目标射频磁场;
根据所述目标射频磁场确定电感值最小时,随钻核磁共振探头天线的天线参数;
根据所述天线参数构建随钻核磁共振探头天线,以使所述随钻核磁共振探头天线产生最优射频磁场,使随钻核磁共振探头能够根据所述最优射频磁场激发共振信号,并接收所述共振信号对应的反馈信号,基于所述反馈信号获取原状地层信息。
在一种可能的实施方式中,所述根据所述目标射频磁场确定电感值最小时,随钻核磁共振探头天线的天线参数,包括:
根据目标射频磁场建立电感函数,并确定所述电感函数对电感的约束条件;
利用所述约束条件对所述电感函数进行约束,确定所述电感函数的电感值最小时随钻核磁共振探头天线的天线参数。
在一种可能的实施方式中,所述根据目标射频磁场建立电感函数,并确定所述电感函数对电感的约束条件,包括:
根据目标射频磁场构建电流密度关于电感的残差约束项,根据所述残差约束项建立电感函数;
确定所述电感函数对电感的约束条件。
在一种可能的实施方式中,所述利用所述约束条件对所述电感函数进行约束,确定所述电感函数的电感值最小时随钻核磁共振探头天线的天线参数,包括:
根据所述约束条件,对电感函数进行偏导计算,确定所述电感函数的最优正则化参数;
根据所述最优正则化参数与所述电感函数,确定所述随钻核磁共振探头天线的天线参数。
在一种可能的实施方式中,所述根据所述约束条件,对电感函数进行偏导计算,确定所述电感函数的最优正则化参数,包括:
根据预设的初始正则化参数,判断电感函数的偏导值是否等于预设阈值;
若是,则此时电感函数对应的正则化参数为最优正则化参数;
若否,则更新所述正则化参数,直至所述电感函数的偏导值等于预设阈值。
在一种可能的实施方式中,若电感函数的偏导值不等于预设阈值,则更新所述正则化参数,直至所述电感函数的偏导值等于预设阈值,包括:
若电感函数的偏导值小于预设阈值,则增大所述正则化参数,直至所述电感函数的偏导值等于预设阈值;
若电感函数的偏导值大于预设阈值,则减小所述正则化参数,直至所述电感函数的偏导值等于预设阈值。
在一种可能的实施方式中,所述根据所述最优正则化参数与所述电感函数,确定所述随钻核磁共振探头天线的天线参数,包括:
根据所述最优正则化参数对所述电感函数进行求解,确定电流密度;
根据所述电流密度对流函数进行求解,确定所述随钻核磁共振探头天线的天线参数。
第二方面,本发明实施例提供一种随钻核磁共振探头天线配置装置,包括:
采集模块,用于采集随钻核磁共振探头的静磁场,根据所述静磁场确定所述静磁场对应的目标射频磁场;
确定模块,用于根据所述目标射频磁场确定电感值最小时,随钻核磁共振探头天线的天线参数;
执行模块,用于根据所述天线参数构建随钻核磁共振探头天线,以使所述随钻核磁共振探头天线产生最优射频磁场,使随钻核磁共振探头能够根据所述最优射频磁场激发共振信号,并接收所述共振信号对应的反馈信号,基于所述反馈信号获取原状地层信息。
第三方面,本发明实施例提供一种随钻核磁共振探头天线配置设备,包括:存储器和至少一个处理器;
所述存储器存储计算机执行指令;
所述至少一个处理器执行所述存储器存储的计算机执行指令,使得所述至少一个处理器执行如第一方面任一项所述的随钻核磁共振探头天线配置方法。
第四方面,本发明实施例提供一种计算机可读存储介质,所述计算机可读存储介质中存储有计算机执行指令,当处理器执行所述计算机执行指令时,实现如第一方面任一项所述的方法。
本发明实施例提供的随钻核磁共振探头天线配置方法、装置、设备及存储介质,通过采集随钻核磁共振探头的静磁场,并确定静磁场对应的目标射频磁场,根据目标射频磁场确定电感值最小时,随钻核磁共振探头天线的天线参数,根据天线参数构建随钻核磁共振探头天线,能够使设计得到的随钻核磁共振探头天线在实际应用中产生与理论中的目标射频磁场相似的最优射频磁场,使随钻核磁共振测井仪可以准确地探测原状地层信息。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种应用场景示意图;
图2为本发明实施例提供的一种随钻核磁共振探头天线配置方法的流程示意图;
图3为本发明实施例提供的一种静磁场数据的示意图;
图4为本发明实施例提供的另一种随钻核磁共振探头天线配置方法的流程示意图;
图5为本发明实施例提供的一种随钻核磁共振探头天线的天线等效模型的结构示意图;
图6为本发明实施例提供的一种随钻核磁共振探头天线配置装置的结构示意图;
图7为本发明实施例提供的一种随钻核磁共振探头天线配置设备的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
石油勘探领域中,相比于电缆核磁共振测井,随钻核磁共振测井作为核磁共振测井的最新发展手段,能够提供更为准确的原状地层信息例如孔隙度、饱和度、渗透率等,以及更加可靠的油藏评价参数例如可动烃含量等。良好的探测特性以及信噪比对随钻核磁共振测井有重要意义。
目前,随钻核磁共振探头天线的主流设计方法是正向设计方法,即应用电磁仿真软件对天线的结构参数进行优化仿真,但该类方法依赖于设计人员的主观经验判断,导致设计得到的天线在实际应用中产生的射频磁场与理论中的目标射频磁场存在较大的差异,进而造成随钻核磁共振测井仪无法准确地探测原状地层信息等。
为了解决上述问题,本实施例采用逆向设计方法对随钻核磁共振探头天线进行设计与配置,即根据目标射频磁场确定随钻核磁共振探头天线的天线参数,进而构建随钻核磁共振探头天线。通过确定与采集到的随钻核磁共振探头的静磁场对应的目标射频磁场,并根据目标射频磁场确定在电感值最小时,随钻核磁共振探头天线的天线参数,根据天线参数构建随钻核磁共振探头天线,能够使得到的随钻核磁共振探头天线在实际应用中产生与理论中的目标射频磁场相似的最优射频磁场,从而可以最大化随钻核磁共振探头所激发的共振信号,使随钻核磁共振测井仪可以准确地探测原状地层信息。
图1为本发明实施例提供的一种应用场景示意图。如图1所示,随钻核磁共振探头天线配置设备采集随钻核磁共振测井仪上的随钻核磁共振探头的相关数据,并根据上述相关数据进行相关运算处理,控制安装在随钻核磁共振探头上的随钻核磁共振探头天线工作。
图2为本发明实施例提供的一种随钻核磁共振探头天线配置方法的流程示意图。本发明实施例中方法的执行主体可以为随钻核磁共振探头天线配置设备。如图2所示,本实施例中的方法,可以包括:
步骤201、采集随钻核磁共振探头的静磁场,根据所述静磁场确定所述静磁场对应的目标射频磁场。
本实施例中,随钻核磁共振探头的静磁场为该探头的永磁体产生的磁场,目标射频磁场为理论中,随钻核磁共振探头天线所要产生的理想射频磁场,该目标射频磁场理论上与探头的永磁体产生的静磁场匹配后产生的共振信号最大。随钻核磁共振探头的静磁场可以通过高斯计进行采集。
其中可选的,执行根据静磁场确定静磁场对应的目标射频磁场的步骤,具体可以为选择过随钻核磁共振测井仪中心轴的某一平面上的静磁场,根据核磁共振条件等,基于目标场法确定与该静磁场对应的目标射频磁场。
图3为本发明实施例提供的一种静磁场数据的示意图。如图3所示,横坐标为目标点距随钻核磁共振测井仪中心轴的径向距离,纵坐标为静磁场强度,表示静磁场在随钻核磁共振测井仪中心轴的径向上的强度分布。目标点为在过随钻核磁共振测井仪中心轴的某一平面上的静磁场中,沿随钻核磁共振测井仪中心轴辐射向外的某一射线上,依据不同距离离散选取的点。
具体的,根据核磁共振条件,基于目标场法,确定与上述静磁场数据对应的目标射频磁场数据,进而确定目标射频磁场。目标射频磁场可以表示为一个向量,该向量与静磁场对应,向量的维度即为目标射频磁场数据的个数。上述核磁共振条件包括:目标射频磁场频率为与其对应的静磁场强度对应的拉莫尔频率;目标射频磁场方向与静磁场方向正交。上述目标场法为一种解析算法,通过目标场法可以推导出磁场与线圈之间的关系。目标射频磁场数据与静磁场数据相对应,选取的静磁场数据的个数越多,即选取的目标点越多,则对应的目标射频磁场数据的个数越多,目标射频磁场的精确度越高,例如选取的静磁场数据的个数可以为600,如果要提高目标射频磁场的精确度,则需要增大选取的静磁场数据的个数。
步骤202、根据所述目标射频磁场确定电感值最小时,随钻核磁共振探头天线的天线参数。
本实施例中,以电感作为约束项参数进行说明,上述约束项参数与天线的线圈储能相关,即约束项参数为天线优化设计实际要考虑的参数,天线的线圈储能与线圈的电感、阻抗等成正相关关系,由此可知,约束项参数可以为电感、阻抗等。当电感值最小时,求取对应的天线参数,根据上述天线参数构建的随钻核磁共振探头天线可以产生与理论中的目标射频磁场相似的最优射频磁场。
进一步的,预设的天线模型结构为螺线管型结构,相应的,上述天线参数为可变参数,可以包括天线的线圈线宽,线圈间距,螺线管直径以及线圈匝数等。根据天线等效模型可以计算确定天线参数。
步骤203、根据所述天线参数构建随钻核磁共振探头天线,以使所述随钻核磁共振探头天线产生最优射频磁场,使随钻核磁共振探头能够根据所述最优射频磁场激发共振信号,并接收所述共振信号对应的反馈信号,基于所述反馈信号获取原状地层信息。
具体的,根据得到的天线参数配置的随钻核磁共振探头天线,能够产生最优射频磁场,该最优射频磁场能与随钻核磁共振探头产生的静磁场匹配激发共振信号。上述最优射频磁场为随钻核磁共振探头天线在实际应用中产生的与理论中的目标射频磁场相似的射频磁场。
本实施例提供的随钻核磁共振探头天线配置方法,通过采集随钻核磁共振探头的静磁场,并确定该静磁场对应的目标射频磁场,根据目标射频磁场确定电感值最小时,随钻核磁共振探头天线的天线参数,根据天线参数构建随钻核磁共振探头天线,能够使设计得到的随钻核磁共振探头天线在实际应用中产生与理论中的目标射频磁场相似的最优射频磁场,使随钻核磁共振测井仪可以准确地探测原状地层信息。
为更准确地在电感值最小时确定随钻核磁共振探头天线的天线参数,本发明实施例还将通过引入电感函数并对其进行偏导计算,确定最优正则化参数来确定电感值最小时,随钻核磁共振探头天线的天线参数。
图4为本发明实施例提供的另一种随钻核磁共振探头天线配置方法的流程示意图。如图4所示,本实施例中的方法,可以包括:
步骤401、采集随钻核磁共振探头的静磁场,根据所述静磁场确定所述静磁场对应的目标射频磁场。
本实施例中步骤401的具体实现过程和原理可以参见前述实施例,此处不再赘述。
步骤402、根据目标射频磁场构建电流密度关于电感的残差约束项,根据所述残差约束项建立电感函数。
具体的,基于目标场法,构建电流密度关于电感的残差约束项,即用电流密度表达电感,电感函数可以包括残差约束项与正则化项。残差约束项为根据目标场法推导出的磁场线性度偏差。
电感函数计算公式为:
Figure BDA0002899724290000071
式中,Mc为电流密度计算得到的射频磁场值,Bz为设定的目标射频磁场,Mc-Bz为残差,λ为正则化参数,W为与电流密度相关的正定矩阵,c为电流密度展开式的傅里叶系数。Φ的第一项为残差约束项,即
Figure BDA0002899724290000081
为残差约束项,第二项为正则化项,即
Figure BDA0002899724290000082
为正则化项。上述残差约束项采用l2范数,正则化项采用Tikhonov正则化项。
步骤403、确定所述电感函数对电感的约束条件。
本实施例中,电感函数对电感的约束条件为在电感值最小时,确定随钻核磁共振探头天线的天线参数。
步骤404、根据所述约束条件,对电感函数进行偏导计算,确定所述电感函数的最优正则化参数。
其中可选的,执行根据所述约束条件,对电感函数进行偏导计算,确定所述电感函数的最优正则化参数的步骤,具体可包括:根据预设的初始正则化参数,判断电感函数的偏导值是否等于预设阈值;若是,则此时电感函数对应的正则化参数为最优正则化参数;若否,则更新所述正则化参数,直至所述电感函数的偏导值等于预设阈值。
本实施例中,电感值最小时对应电感函数取极小值时,即可对电感函数进行偏导计算,在电感函数取极小值时确定电感函数的最优正则化参数。
具体的,更改正则化参数的值,当电感函数取极小值或接近极小值时,电感函数所对应的正则化参数即为最优正则化参数,即电感函数的偏导值为0或接近0时,电感函数所对应的正则化参数即为最优正则化参数。
初始正则化参数可以根据经验进行预设,以减少迭代的次数。电感函数随正则化参数变化的斜率公式,即电感函数对电流密度展开式的傅里叶系数c的偏导公式为:
Figure BDA0002899724290000083
进一步的,
Figure BDA0002899724290000084
的值即为电感函数的偏导值,将
Figure BDA0002899724290000085
的值与预设阈值进行比较,预设阈值可以为大于0且小于1的数值,例如,预设阈值可以为0.1。
若本次
Figure BDA0002899724290000086
的值等于预设阈值,则取此时电感函数对应的正则化参数为最优正则化参数。
若本次
Figure BDA0002899724290000087
的值小于预设阈值,则增大本次正则化参数,例如将本次正则化参数乘以某一大于1的倍数,该倍数可以为预设阈值除以本次
Figure BDA0002899724290000088
的值得到的数值,以进行下一次迭代,直至
Figure BDA0002899724290000091
的值等于预设阈值。
若本次
Figure BDA0002899724290000092
的值大于预设阈值,则减小本次正则化参数,例如将本次正则化参数乘以某一小于1的倍数,该倍数可以为预设阈值除以本次
Figure BDA0002899724290000093
的值得到的数值,以进行下一次迭代,直至
Figure BDA0002899724290000094
的值等于预设阈值。
步骤405、根据所述最优正则化参数对所述电感函数进行求解,确定电流密度。
具体的,将最优正则化参数带入电感函数中,对电感函数进行求解,确定电流密度展开式的傅里叶系数c,进而根据电流密度展开式的傅里叶系数c确定对应的电流密度j。
电流密度j的计算公式为:
Figure BDA0002899724290000095
Figure BDA0002899724290000096
其中,jθ(θ′,z′)表示电流密度j在角向方向θ的分量,jz(θ′,z′)表示电流密度j在轴向方向z的分量,m、n表示傅里叶级数项数,L表示线圈的高度,a表示电流密度j分布的圆柱面半径。
步骤406、根据所述电流密度对流函数进行求解,确定所述随钻核磁共振探头天线的天线参数。
对电流密度j直接求解较为复杂,为简化计算流程,引入流函数,根据电流密度j对流函数进行求解,上述电流密度j为矢量,流函数为标量。
电流密度j与流函数Ψ(θ′,z′)的关系可表示为:
Figure BDA0002899724290000097
其中,er表示矢量势的方向,r′为原点到对应的电流密度位置的位置矢量。由公式(5)可得到电流密度j在角向方向θ的分量和轴向方向z的分量:
Figure BDA0002899724290000098
Figure BDA0002899724290000099
根据公式(3)-(7),得到流函数Ψ(θ′,z′)公式:
Figure BDA00028997242900000910
对流函数Ψ(θ′,z′)进行求解得到流函数的傅里叶级数表达式,进而可以得到等值线路径,等值线路径即为天线的绕组形式。根据天线的绕组形式进一步得到天线等效模型,即天线等效电流环模型,进而根据上述天线等效电流环模型计算得到天线参数。
上述具体实现方式为本领域现有技术,本申请对此不进行具体限制。
步骤407、根据所述天线参数构建随钻核磁共振探头天线,以使所述随钻核磁共振探头天线产生最优射频磁场,使随钻核磁共振探头能够根据所述最优射频磁场激发共振信号,并接收所述共振信号对应的反馈信号,基于所述反馈信号获取原状地层信息。
本实施例中步骤407的具体实现过程和原理可以参见前述实施例,此处不再赘述。
图5为本发明实施例提供的一种随钻核磁共振探头天线的天线等效模型的结构示意图。随钻核磁共振探头天线为螺线管型结构,如图5所示,为方便观察以及计算,将随钻核磁共振探头天线的螺线线圈形式等效为独立线圈形式。随钻核磁共振探头天线的天线参数为可变参数,可以包括线圈线宽,线圈间距,螺线管直径与线圈匝数等,根据天线等效模型可以计算得出天线参数,根据天线参数可以确定天线线圈在圆柱面上的分布位置,即确定各匝天线的绕线位置。
本实施例提供的随钻核磁共振探头天线配置方法,通过采集随钻核磁共振探头的静磁场,确定该静磁场对应的目标射频磁场,然后根据电流密度与电感的关系构造残差约束项,进而建立电感函数,确定电感函数对电感的约束条件,根据约束条件对电感函数进行偏导计算,并基于迭代法确定电感函数的最优正则化参数,然后根据最优正则化参数对电感函数进行求解,确定电流密度,并根据电流密度对流函数进行求解,确定随钻核磁共振探头天线的天线参数,根据该天线参数构建随钻核磁共振探头天线,能够通过逆向设计方法对随钻核磁共振探头天线进行设计,使设计得到的随钻核磁共振探头天线在实际应用中产生与理论中的目标射频磁场相似的最优射频磁场,提高随钻核磁共振探头的信噪比,进而使随钻核磁共振测井仪可以准确地探测原状地层信息。
图6为本发明实施例提供的一种随钻核磁共振探头天线配置装置的结构示意图。如图6所示,本实施例提供的随钻核磁共振探头天线配置装置,可以包括:采集模块61、确定模块62以及执行模块63。
采集模块61,用于采集随钻核磁共振探头的静磁场,根据所述静磁场确定所述静磁场对应的目标射频磁场;
确定模块62,用于根据所述目标射频磁场确定电感值最小时,随钻核磁共振探头天线的天线参数;
执行模块63,用于根据所述天线参数构建随钻核磁共振探头天线,以使所述随钻核磁共振探头天线产生最优射频磁场,使随钻核磁共振探头能够根据所述最优射频磁场激发共振信号,并接收所述共振信号对应的反馈信号,基于所述反馈信号获取原状地层信息。
在一种可选的实现方式中,所述确定模块62具体用于:
根据目标射频磁场建立电感函数,并确定所述电感函数对电感的约束条件;
利用所述约束条件对所述电感函数进行约束,确定所述电感函数的电感值最小时随钻核磁共振探头天线的天线参数。
在一种可选的实现方式中,所述确定模块62在根据目标射频磁场建立电感函数,并确定所述电感函数对电感的约束条件时,具体用于:
根据目标射频磁场构建电流密度关于电感的残差约束项,根据所述残差约束项建立电感函数;
确定所述电感函数对电感的约束条件。
在一种可选的实现方式中,所述确定模块62在利用所述约束条件对所述电感函数进行约束,确定所述电感函数的电感值最小时随钻核磁共振探头天线的天线参数时,具体用于:
根据所述约束条件,对电感函数进行偏导计算,确定所述电感函数的最优正则化参数;
根据所述最优正则化参数与所述电感函数,确定所述随钻核磁共振探头天线的天线参数。
在一种可选的实现方式中,所述确定模块62在根据所述约束条件,对电感函数进行偏导计算,确定所述电感函数的最优正则化参数时,具体用于:
根据预设的初始正则化参数,判断电感函数的偏导值是否等于预设阈值;
若是,则此时电感函数对应的正则化参数为最优正则化参数;
若否,则更新所述正则化参数,直至所述电感函数的偏导值等于预设阈值。
在一种可选的实现方式中,所述确定模块62在电感函数的偏导值不等于预设阈值,更新所述正则化参数,直至所述电感函数的偏导值等于预设阈值时,具体用于:
若电感函数的偏导值小于预设阈值,则增大所述正则化参数,直至所述电感函数的偏导值等于预设阈值;
若电感函数的偏导值大于预设阈值,则减小所述正则化参数,直至所述电感函数的偏导值等于预设阈值。
在一种可选的实现方式中,所述确定模块62在根据所述最优正则化参数与所述电感函数,确定所述随钻核磁共振探头天线的天线参数时,具体用于:
根据所述最优正则化参数对所述电感函数进行求解,确定电流密度;
根据所述电流密度对流函数进行求解,确定所述随钻核磁共振探头天线的天线参数。
本实施例提供的随钻核磁共振探头天线配置装置,可以执行上述方法实施例的技术方案,其实现原理和技术效果类似,此处不再赘述。
图7为本发明实施例提供的一种随钻核磁共振探头天线配置设备的结构示意图。如图7所示,本实施例提供的随钻核磁共振探头天线配置设备,包括:存储器71和至少一个处理器72;
所述存储器71存储计算机执行指令;
所述至少一个处理器72执行所述存储器71存储的计算机执行指令,使得所述至少一个处理器执行如上述任一实施例所述的随钻核磁共振探头天线配置方法。
其中,存储器71和处理器72可以通过总线73连接。
本实施例提供的随钻核磁共振探头天线配置设备的具体实现原理和效果可以参见图1-图5所示实施例对应的相关描述和效果,此处不做过多赘述。
本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质中存储有计算机执行指令,当处理器执行所述计算机执行指令时,实现如上述任一实施例所述的方法。
其中,计算机可读存储介质可以是ROM、随机存取存储器(RAM)、CD-ROM、磁带、软盘和光数据存储设备等。
在本发明所提供的几个实施例中,应该理解到,所揭露的装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,模块的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个模块或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或模块的间接耦合或通信连接,可以是电性,机械或其它的形式。
上述存储介质可以是由任何类型的易失性或非易失性存储设备或者它们的组合实现,如静态随机存取存储器(SRAM),电可擦除可编程只读存储器(EEPROM),可擦除可编程只读存储器(EPROM),可编程只读存储器(PROM),只读存储器(ROM),磁存储器,快闪存储器,磁盘或光盘。存储介质可以是通用或专用计算机能够存取的任何可用介质。
一种示例性的存储介质耦合至处理器,从而使处理器能够从该存储介质读取信息,且可向该存储介质写入信息。当然,存储介质也可以是处理器的组成部分。处理器和存储介质可以位于专用集成电路(Application Specific Integrated Circuits,简称ASIC)中。当然,处理器和存储介质也可以作为分立组件存在于电子设备或主控设备中。
最后应说明的是:以上各实施例仅用以说明本申请的技术方案,而非对其限制;尽管参照前述各实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本申请各实施例技术方案的范围。

Claims (7)

1.一种随钻核磁共振探头天线配置方法,其特征在于,所述随钻核磁共振探头天线配置方法应用于随钻核磁共振探头天线配置装置,所述配置装置用于控制安装在随钻核磁共振探头上的随钻核磁共振探头天线工作,以使所述随钻核磁共振探头根据所述随钻核磁共振探头天线产生的射频磁场激发共振信号;
所述方法包括:
采集随钻核磁共振探头的静磁场,根据所述静磁场确定所述静磁场对应的目标射频磁场;
根据所述目标射频磁场确定电感值最小时,随钻核磁共振探头天线的天线参数,包括:
根据目标射频磁场建立电感函数,并确定所述电感函数对电感的约束条件;
所述根据目标射频磁场建立电感函数,并确定所述电感函数对电感的约束条件,包括:根据目标射频磁场构建电流密度关于电感的残差约束项,根据所述残差约束项建立电感函数;确定所述电感函数对电感的约束条件;
其中,电感函数计算公式为:
Figure FDA0004065728750000011
Figure FDA0004065728750000012
Mc为电流密度计算得到的射频磁场值,Bz为设定的目标射频磁场,Mc-Bz为残差,λ为正则化参数,W为与电流密度相关的正定矩阵,c为电流密度展开式的傅里叶系数,
Figure FDA0004065728750000013
为残差约束项,
Figure FDA0004065728750000014
为正则化项;
利用所述约束条件对所述电感函数进行约束,确定所述电感函数的电感值最小时随钻核磁共振探头天线的天线参数,包括:
根据所述约束条件,对电感函数进行偏导计算,确定所述电感函数的最优正则化参数;其中,所述电感函数的偏导公式为:
Figure FDA0004065728750000015
Figure FDA0004065728750000021
Figure FDA0004065728750000022
的值即为电感函数的偏导值,MT为磁化矢量的转置矩阵,WT为与电流密度相关的正定矩阵的转置矩阵;
根据所述最优正则化参数与所述电感函数,确定所述随钻核磁共振探头天线的天线参数;
根据所述天线参数构建随钻核磁共振探头天线,以使所述随钻核磁共振探头天线产生最优射频磁场,使随钻核磁共振探头能够根据所述最优射频磁场激发共振信号,并接收所述共振信号对应的反馈信号,基于所述反馈信号获取原状地层信息。
2.根据权利要求1所述的方法,其特征在于,所述根据所述约束条件,对电感函数进行偏导计算,确定所述电感函数的最优正则化参数,包括:
根据预设的初始正则化参数,判断电感函数的偏导值是否等于预设阈值;
若是,则此时电感函数对应的正则化参数为最优正则化参数;
若否,则更新所述正则化参数,直至所述电感函数的偏导值等于预设阈值。
3.根据权利要求2所述的方法,其特征在于,若电感函数的偏导值不等于预设阈值,则更新所述正则化参数,直至所述电感函数的偏导值等于预设阈值,包括:
若电感函数的偏导值小于预设阈值,则增大所述正则化参数,直至所述电感函数的偏导值等于预设阈值;
若电感函数的偏导值大于预设阈值,则减小所述正则化参数,直至所述电感函数的偏导值等于预设阈值。
4.根据权利要求1所述的方法,其特征在于,所述根据所述最优正则化参数与所述电感函数,确定所述随钻核磁共振探头天线的天线参数,包括:
根据所述最优正则化参数对所述电感函数进行求解,确定电流密度;
根据所述电流密度对流函数进行求解,确定所述随钻核磁共振探头天线的天线参数。
5.一种随钻核磁共振探头天线配置装置,其特征在于,包括:
采集模块,用于采集随钻核磁共振探头的静磁场,根据所述静磁场确定所述静磁场对应的目标射频磁场;
确定模块,用于根据所述目标射频磁场确定电感值最小时,随钻核磁共振探头天线的天线参数;
执行模块,用于根据所述天线参数构建随钻核磁共振探头天线,以使所述随钻核磁共振探头天线产生最优射频磁场,使随钻核磁共振探头能够根据所述最优射频磁场激发共振信号,并接收所述共振信号对应的反馈信号,基于所述反馈信号获取原状地层信息;
所述确定模块,具体用于根据目标射频磁场建立电感函数,并确定所述电感函数对电感的约束条件;所述根据目标射频磁场建立电感函数,并确定所述电感函数对电感的约束条件,包括:根据目标射频磁场构建电流密度关于电感的残差约束项,根据所述残差约束项建立电感函数;确定所述电感函数对电感的约束条件;
其中,电感函数计算公式为:
Figure FDA0004065728750000031
Figure FDA0004065728750000032
Mc为电流密度计算得到的射频磁场值,Bz为设定的目标射频磁场,Mc-Bz为残差,λ为正则化参数,W为与电流密度相关的正定矩阵,c为电流密度展开式的傅里叶系数,
Figure FDA0004065728750000033
为残差约束项,
Figure FDA0004065728750000034
为正则化项;
根据所述约束条件,对电感函数进行偏导计算,确定所述电感函数的最优正则化参数;其中,所述电感函数的偏导公式为:
Figure FDA0004065728750000035
Figure FDA0004065728750000036
Figure FDA0004065728750000037
的值即为电感函数的偏导值,MT为磁化矢量的转置矩阵,WT为与电流密度相关的正定矩阵的转置矩阵;
根据所述最优正则化参数与所述电感函数,确定所述随钻核磁共振探头天线的天线参数。
6.一种随钻核磁共振探头天线配置设备,其特征在于,包括:存储器和至少一个处理器;
所述存储器存储计算机执行指令;
所述至少一个处理器执行所述存储器存储的计算机执行指令,使得所述至少一个处理器执行如权利要求1-4任一项所述的随钻核磁共振探头天线配置方法。
7.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质中存储有计算机执行指令,当处理器执行所述计算机执行指令时,实现如权利要求1-4任一项所述的方法。
CN202110053137.3A 2021-01-15 2021-01-15 随钻核磁共振探头天线配置方法、装置、设备及存储介质 Active CN112855118B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110053137.3A CN112855118B (zh) 2021-01-15 2021-01-15 随钻核磁共振探头天线配置方法、装置、设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110053137.3A CN112855118B (zh) 2021-01-15 2021-01-15 随钻核磁共振探头天线配置方法、装置、设备及存储介质

Publications (2)

Publication Number Publication Date
CN112855118A CN112855118A (zh) 2021-05-28
CN112855118B true CN112855118B (zh) 2023-03-14

Family

ID=76006504

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110053137.3A Active CN112855118B (zh) 2021-01-15 2021-01-15 随钻核磁共振探头天线配置方法、装置、设备及存储介质

Country Status (1)

Country Link
CN (1) CN112855118B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NO976053D0 (no) * 1997-01-29 1997-12-22 Western Atlas Int Inc Kombinert kjernemagnetisk resonans- og elektromagnetisk induksjonresistivitets-brönnloggeinstrument og -fremgangsmåte
CN101996273A (zh) * 2010-11-29 2011-03-30 浙江大学 Mri系统梯度线圈的有限差分设计方法
CN102376410A (zh) * 2011-11-10 2012-03-14 中国科学院电工研究所 一种开放式核磁共振超导磁体的设计方法
CN105718677A (zh) * 2016-01-22 2016-06-29 中国科学院电工研究所 自屏蔽超导核磁共振成像系统梯度线圈设计方法
CN106291723A (zh) * 2016-07-25 2017-01-04 中国石油大学(北京) 基于双参数正则化的核磁共振回波数据反演方法及装置
CN108872896A (zh) * 2018-05-29 2018-11-23 河北惠仁医疗设备科技有限公司 一种双平面磁共振成像系统梯度线圈的设计方法
CN109031437A (zh) * 2018-08-10 2018-12-18 中国石油大学(北京) 随钻核磁共振自旋回波校正方法与装置

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2983448C (en) * 2015-05-12 2021-03-09 Hyperfine Research, Inc. Radio frequency coil methods and apparatus
US10261210B2 (en) * 2016-03-09 2019-04-16 Baker Hughes, A Ge Company, Llc Method and apparatus for active suppression of pipe signals in transient electromagnetic measurements
US10802176B2 (en) * 2017-02-14 2020-10-13 Baker Hughes, A Ge Company, Llc Methods and devices for magnetic resonance measurements using decoupled transmit antennas

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NO976053D0 (no) * 1997-01-29 1997-12-22 Western Atlas Int Inc Kombinert kjernemagnetisk resonans- og elektromagnetisk induksjonresistivitets-brönnloggeinstrument og -fremgangsmåte
CN101996273A (zh) * 2010-11-29 2011-03-30 浙江大学 Mri系统梯度线圈的有限差分设计方法
CN102376410A (zh) * 2011-11-10 2012-03-14 中国科学院电工研究所 一种开放式核磁共振超导磁体的设计方法
CN105718677A (zh) * 2016-01-22 2016-06-29 中国科学院电工研究所 自屏蔽超导核磁共振成像系统梯度线圈设计方法
CN106291723A (zh) * 2016-07-25 2017-01-04 中国石油大学(北京) 基于双参数正则化的核磁共振回波数据反演方法及装置
CN108872896A (zh) * 2018-05-29 2018-11-23 河北惠仁医疗设备科技有限公司 一种双平面磁共振成像系统梯度线圈的设计方法
CN109031437A (zh) * 2018-08-10 2018-12-18 中国石油大学(北京) 随钻核磁共振自旋回波校正方法与装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
实际电磁场环境中低场脉冲核磁共振信号的分析计算;骆逸锋;《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》;20041215;1-5,41-44 *
骆逸锋.实际电磁场环境中低场脉冲核磁共振信号的分析计算.《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》.2004, *

Also Published As

Publication number Publication date
CN112855118A (zh) 2021-05-28

Similar Documents

Publication Publication Date Title
CN103323801B (zh) 磁共振系统控制序列的确定
US6393363B1 (en) Method and apparatus for the measurement of the electrical resistivity of geologic formations employing modeling data
CA2488144C (en) Method for determining sonde error for an induction or propagation tool with transverse or triaxial arrays
EP2105765A1 (en) Simultaneous inversion of induction data for dielectric permittivity and electric conductivity
US20050156602A1 (en) System, apparatus, and method for conducting electromagnetic induction surveys
WO2016184779A1 (en) Time-domain mri
US10359539B2 (en) Calibration method for multi-component induction tools
CN105240001A (zh) 核磁共振测井孔隙度校正方法及装置
US10001578B2 (en) Nuclear magnetic resonance refocusing pulses for inhomogeneous magnetic fields
US10254433B2 (en) Nuclear magnetic resonance (NMR) logging tool calibration
US11422282B2 (en) Circular polarization correction in nuclear magnetic resonance (NMR) logging
CN112855118B (zh) 随钻核磁共振探头天线配置方法、装置、设备及存储介质
CN110029990B (zh) 一种核磁共振测井方法和装置
CN109828307B (zh) 一种瞬变电磁多频率融合的探测方法及应用
US11143781B2 (en) Accounting for tool based effects in nuclear magnetic resonance logging data
Colinas et al. Application of the Tikhonov regularization technique to the equivalent magnetic currents near-field technique
CN111965719A (zh) 一种弛豫时间测量方法和装置
US11655702B2 (en) Space mapping optimization to characterize multiple concentric pipes
CN112773396A (zh) 基于全波形反演的医学成像方法、计算机设备及存储介质
US8694259B2 (en) Simultaneous inversion of induction data for dielectric permittivity and electric conductivity
CA2383245C (en) Method and apparatus of electrical resistivity measurements in geological formations employing modeling data
US11567157B2 (en) Calibration of a magnetic resonance device and estimating a specific absorption rate
RU2431859C2 (ru) Способ и устройство для определения дефектов изготовления, сборки и установки магнитных систем
US10408964B2 (en) Multiple detector NMR tool system and inversion methods for NMR logging
Taybi et al. Geometrical and Analytical Methods for Probes Correction applied to Near Field Antennas Measurement Systems

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