CN117434614A - 一种地质体的总磁化方向的确定方法、装置、设备和介质 - Google Patents

一种地质体的总磁化方向的确定方法、装置、设备和介质 Download PDF

Info

Publication number
CN117434614A
CN117434614A CN202311426553.9A CN202311426553A CN117434614A CN 117434614 A CN117434614 A CN 117434614A CN 202311426553 A CN202311426553 A CN 202311426553A CN 117434614 A CN117434614 A CN 117434614A
Authority
CN
China
Prior art keywords
magnetic
magnetization direction
current
determining
component 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.)
Granted
Application number
CN202311426553.9A
Other languages
English (en)
Other versions
CN117434614B (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 Aero Geophysical Survey and Remote Sensing Center for Natural Resources
Original Assignee
China Aero Geophysical Survey and Remote Sensing Center for Natural Resources
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 Aero Geophysical Survey and Remote Sensing Center for Natural Resources filed Critical China Aero Geophysical Survey and Remote Sensing Center for Natural Resources
Priority to CN202311426553.9A priority Critical patent/CN117434614B/zh
Publication of CN117434614A publication Critical patent/CN117434614A/zh
Application granted granted Critical
Publication of CN117434614B publication Critical patent/CN117434614B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/40—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for measuring magnetic field characteristics of the earth
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Physics & Mathematics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种地质体的总磁化方向的确定方法、装置、设备和介质。该方法包括：获取目标地质体的实测磁三分量数据，确定当前次处理周期使用的当前磁化方向，基于第一参考磁化方向对实测磁三分量数据中第三方向磁分量数据进行化极处理获得化极磁异常数据，基于当前磁化方向和实测磁三分量数据确定当前磁模量，并基于当前磁模量与化极磁异常数据确定下一磁化方向，进一步判断第一参数是否满足预设条件，从而确定是否对磁化方向做进一步的更新以准确确定目标地质体的总磁化方向。本申请技术方案解决了因分量转化而导致磁异常模量计算不准确的问题，也避免了引入由观测点磁场方向不准确带来的系统误差，提高了地质体的总磁化方向的准确性。

Description

一种地质体的总磁化方向的确定方法、装置、设备和介质
技术领域
本发明涉及磁力勘探技术领域，尤其涉及一种地质体的总磁化方向的确定方法、装置、设备和介质。
背景技术
磁力勘探应用广泛，尤其在隐伏富铁矿勘查方面具有很好的应用效果。准确而快速地估计地下磁性体的总磁化方向，对于地球物理磁测资料的解释、处理和反演具有重要意义。
目前互相关法是计算磁化方向的方法中比较快捷和稳定的方法，这种方法一般会通过计算不同总磁化方向下得到的化极场与选定的磁异常模量之间的线性相关性来获取最佳的磁化方向，但是此种方式需要对磁异常模量进行分量转化，无法保证磁异常模量的准确性，且磁倾角对磁异常模量的中心对称性有较大的影响，会导致磁异常模量的极大值向磁化方向的反方向偏移，从而导致总磁化方向的确定不准确。
发明内容
本发明提供了一种地质体的总磁化方向的确定方法、装置、设备和介质，以解决因分量转化而导致磁异常模量计算不准确的问题，也避免了引入由观测点磁场方向不准确带来的系统误差，提高了地质体的总磁化方向的准确性。
根据本发明的一方面，提供了一种地质体的总磁化方向的确定方法，该方法包括：
获取目标地质体的实测磁三分量数据，所述实测磁三分量数据包括第一方向磁分量数据、第二方向磁分量数据和第三方向磁分量数据，所述第一方向为实测时正北方向，第二方向为实测时正东方向，第三方向为实测时垂直向下的方向；
确定当前次处理周期使用的当前磁化方向，所述当前磁化方向包括初始磁偏角和当前磁倾角，所述初始磁偏角为初始磁化方向包含的一个方位角度，所述初始磁化方向为初次根据所述实测磁三分量数据所确定的磁化方向；
基于第一参考磁化方向对所述实测磁三分量数据中第三方向磁分量数据进行化极处理获得化极磁异常数据，所述第一参考磁化方向包括所述初始磁偏角和预设磁倾角；
基于所述当前磁化方向和所述实测磁三分量数据确定当前磁模量，并基于所述当前磁模量与所述化极磁异常数据确定下一磁化方向，以供进入下一次处理周期加载使用所述下一磁化方向，所述下一磁化方向包括初始磁偏角和下一磁倾角；
确定第一参数，判断第一参数是否满足预设条件，所述第一参数为所述当前磁倾角与下一磁倾角的差值的绝对值，所述预设条件为所述第一参数小于或等于预设值；
若所述第一参数不满足所述预设条件，则将所述下一磁化方向作为当前磁化方向，并基于所述当前磁化方向和所述实测磁三分量数据确定当前磁模量，以及基于所述化极磁异常数据与所述当前磁模量确定下一磁化方向，直至所述第一参数满足预设条件；
若所述第一参数满足所述预设条件，则将下一磁化方向作为所述目标地质体的总磁化方向。
根据本发明的另一方面，提供了一种地质体的总磁化方向的确定装置，该装置包括：
第一数据确定模块，用于获取目标地质体的实测磁三分量数据，所述实测磁三分量数据包括第一方向磁分量数据、第二方向磁分量数据和第三方向磁分量数据，所述第一方向为实测时正北方向，第二方向为实测时正东方向，第三方向为实测时垂直向下的方向；
磁化方向确定模块，用于确定当前次处理周期使用的当前磁化方向，所述当前磁化方向包括初始磁偏角和当前磁倾角，所述初始磁偏角为初始磁化方向包含的一个方位角度，所述初始磁化方向为初次根据所述实测磁三分量数据所确定的磁化方向；
第二数据确定模块，用于基于第一参考磁化方向对所述实测磁三分量数据中第三方向磁分量数据进行化极处理获得化极磁异常数据，所述第一参考磁化方向包括所述初始磁偏角和预设磁倾角；
信息确定模块，用于基于所述当前磁化方向和所述实测磁三分量数据确定当前磁模量，并基于所述当前磁模量与所述化极磁异常数据确定下一磁化方向，以供进入下一次处理周期加载使用所述下一磁化方向，所述下一磁化方向包括初始磁偏角和下一磁倾角；
判断模块，用于确定第一参数，判断第一参数是否满足预设条件，所述第一参数为所述当前磁倾角与下一磁倾角的差值的绝对值，所述预设条件为所述第一参数小于或等于预设值；
第一条件确定模块，用于若所述第一参数不满足所述预设条件，则将所述下一磁化方向作为当前磁化方向，并基于所述当前磁化方向和所述实测磁三分量数据确定当前磁模量，以及基于所述化极磁异常数据与所述当前磁模量确定下一磁化方向，直至所述第一参数满足预设条件；
第二条件确定模块，用于若所述第一参数满足所述预设条件，则将下一磁化方向作为所述目标地质体的总磁化方向。
根据本发明的另一方面，提供了一种电子设备，所述电子设备包括：
至少一个处理器；以及
与所述至少一个处理器通信连接的存储器；其中，
所述存储器存储有可被所述至少一个处理器执行的计算机程序，所述计算机程序被所述至少一个处理器执行，以使所述至少一个处理器能够执行本发明任一实施例所述的一种地质体的总磁化方向的确定方法。
根据本发明的另一方面，提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于使处理器执行时实现本发明任一实施例所述的一种地质体的总磁化方向的确定方法。
本发明实施例的技术方案，确定目标地质体的实测磁三分量数据，然后确定当前次处理周期使用的当前磁化方向，基于第一参考磁化方向对实测磁三分量数据中第三方向磁分量数据进行化极处理获得化极磁异常数据，基于当前磁化方向和实测磁三分量数据确定当前磁模量，并基于当前磁模量与化极磁异常数据确定下一磁化方向，进一步判断第一参数是否满足预设条件，第一参数为当前磁倾角与下一磁倾角的差值的绝对值，预设条件为第一参数小于或等于预设值，从而确定是否对磁化方向做进一步的更新以准确确定目标地质体的总磁化方向，解决了因分量转化而导致磁异常模量计算不准确的问题，也避免了引入由观测点磁场方向不准确带来的系统误差，提高了地质体的总磁化方向的准确性。
应当理解，本部分所描述的内容并非旨在标识本发明的实施例的关键或重要特征，也不用于限制本发明的范围。本发明的其它特征将通过以下的说明书而变得容易理解。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是根据本发明实施例提供的一种地质体的总磁化方向的确定方法的流程图；
图2是根据本发明实施例提供的一种地质体的总磁化方向的确定方法的流程图；
图3是根据本发明实施例所适用的实测磁三分量数据的示例图；
图4是根据本发明实施例所适用的磁异常模量数据的示例图；
图5是根据本发明实施例所适用的参考磁化方向以网格形式表示的示例图；
图6是根据本发明实施例所适用的第二相关系数分布图；
图7是根据本发明实施例所适用的第一相关系数分布图；
图8是根据本发明实施例提供的一种地质体的总磁化方向的确定装置的结构示意图；
图9是根据本发明实施例提供的实现本发明实施例的一种地质体的总磁化方向的确定方法的电子设备的结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
实施例一
图1为本发明实施例提供的一种地质体的总磁化方向的确定方法的流程图，本实施例可适用于对地质体的总磁化方向进行确定的情况，该方法可以由地质体的总磁化方向的确定装置来执行，该地质体的总磁化方向的确定装置可以采用硬件和/或软件的形式实现，该地质体的总磁化方向的确定装置可配置于任何具有网络通信功能的电子设备中。如图1所示，该方法包括：
S110、获取目标地质体的实测磁三分量数据。
其中，实测磁三分量数据包括第一方向磁分量数据、第二方向磁分量数据和第三方向磁分量数据，第一方向为实测时正北方向，第二方向为实测时正东方向，第三方向为实测时垂直向下的方向。第一方向、第二方向和第三方向互相垂直。
具体的，实测磁三分量数据包括的第一方向磁分量数据、第二方向磁分量数据和第三方向磁分量数据可以分别用ΔBx，ΔBy和ΔBz表示。
S120、确定当前次处理周期使用的当前磁化方向。
其中，当前磁化方向包括初始磁偏角和当前磁倾角，初始磁偏角为初始磁化方向包含的一个方位角度，初始磁化方向为初次根据实测磁三分量数据所确定的磁化方向。
具体的，当前磁化方向决定如何对实测磁三分量数据进行垂直磁化处理，因此如何确定当前次处理周期使用的当前磁化方向十分重要，具体包括如下步骤A1-A2：
步骤A1、若当前次处理周期为第一个处理周期，则当前磁化方向为初始磁化方向，初始磁化方向包括初始磁偏角和初始磁倾角。
具体的，初始磁化方向的确定过程，包括步骤B1-B2：
步骤B1、基于实测磁三分量数据确定磁异常模量，并基于第二参考磁化方向对第三方向磁分量数据进行化极处理获得初始化极磁异常数据，第二参考磁化方向包括预设磁偏角和预设磁倾角。
具体的，获取ΔBx，ΔBy和ΔBz，将实测磁三分量数据按照如下公式求得磁异常模量Ta：
基于第二参考磁化方向确定初始波数域转换因子基于初始波数域转换因子，对第三方向磁分量数据ΔBz进行化极处理获得初始化极磁异常数据ΔB0 z⊥。其中，初始波数域转换因子可如下表示：
其中，ωx和ωy分别为x方向(正北方向)和y方向(正东方向)的波数， 为总磁化方向的单位矢量，即/>其中d、i为假定的磁化方向的预设磁偏角、预设磁倾角，即第二参考磁化方向，预设磁偏角在0°到360°之间，预设磁倾角在-90°到90°之间。与常规的总磁场异常ΔT的转换因子不同，该转换因子不需要输入环境地磁场的方向，即不受总磁场方向的影响。
步骤B2、基于磁异常模量和初始化极磁异常数据确定初始磁化方向，初始磁化方向包括初始磁偏角d0和初始磁倾角i0。
具体的，确定第二关联关系，基于磁异常模量、初始化极磁异常数据和第二关联关系确定初始相关系数，并将初始相关系数对应的磁化方向作为初始磁化方向(d0,i0)，初始相关系数为第二相关系数中的最大值，第二关联关系为第二相关系数磁异常模量Ta和初始化极磁异常数据ΔB0 z⊥之间的相关关系。其中，第二关联关系可如下表示：
且初始磁化方向确定后，磁偏角是准确的，而在三维情况下，磁倾角对磁异常模量Ta的中心对称性的影响是不能忽略的，因此需要对磁倾角进行进一步的修正，即此后周期的磁偏角一直为初始磁偏角，而磁倾角需进一步的修正。也即如此第一参考磁化方向是由初始磁偏角和预设磁倾角假定，预设磁倾角在-90°到90°之间。
步骤A2、若当前次处理周期不为第一个处理周期，则当前磁化方向为上一个处理周期确定的磁化方向。
S130、基于第一参考磁化方向对实测磁三分量数据中第三方向磁分量数据进行化极处理获得化极磁异常数据，第一参考磁化方向包括初始磁偏角和预设磁倾角。
具体的，基于第一参考磁化方向确定波数域转换因子，基于波数域转换因子，对第三方向磁分量数据进行化极处理获得化极磁异常数据ΔBz⊥，且此后循环过程中化极磁异常数据一直保持不变，其中，波数域转换因子可如下表示：
S140、基于当前磁化方向和实测磁三分量数据确定当前磁模量，并基于当前磁模量与化极磁异常数据确定下一磁化方向，以供进入下一次处理周期加载使用下一磁化方向。
其中，下一磁化方向包括初始磁偏角和下一磁倾角。
具体的，基于当前磁化方向(d0,i)对实测磁三分量数据进行垂直磁化处理得到垂直磁化磁三分量数据(和/>)，对垂直磁化磁三分量数据求平方和后开方获得当前磁模量/>具体可表示如下：
其中，和/>分别是基于当前磁化方向对ΔBx，ΔBy和ΔBz进行垂直磁化所得。
进一步的确定第一关联关系，基于化极磁异常数据、当前磁模量和第一关联关系确定当前相关系数，并将当前相关系数对应的磁化方向作为下一磁化方向，当前相关系数为第一相关系数中的最大值，第一关联关系为第一相关系数化极磁异常数据ΔBz⊥和当前磁模量/>之间的相关关系；其中，第一关联关系可如下表示：
S150、确定第一参数，判断所述第一参数是否满足预设条件，第一参数为所述当前磁倾角与下一磁倾角的差值的绝对值，预设条件为第一参数小于或等于预设值。
其中，预设值是根据实际情况进行设定的，且当前磁倾角与下一磁倾角的差值足够小才能说明现在的磁倾角精度已经足够，满足所需总磁化方向的确定。
S160、若第一参数不满足预设条件，则将下一磁化方向作为当前磁化方向，并基于当前磁化方向和实测磁三分量数据确定当前磁模量，以及基于化极磁异常数据与当前磁模量确定下一磁化方向，直至第一参数满足预设条件。
S170、若第一参数满足所述预设条件，则将下一磁化方向作为目标地质体的总磁化方向。
本发明实施例的技术方案，获取目标地质体的实测磁三分量数据，然后确定当前次处理周期使用的当前磁化方向，基于第一参考磁化方向对实测磁三分量数据中第三方向磁分量数据进行化极处理获得化极磁异常数据，基于当前磁化方向和实测磁三分量数据确定当前磁模量，并基于当前磁模量与化极磁异常数据确定下一磁化方向，进一步判断第一参数是否满足预设条件，第一参数为当前磁倾角与下一磁倾角的差值的绝对值，预设条件为第一参数小于或等于预设值，从而确定是否对磁化方向做进一步的更新以准确确定目标地质体的总磁化方向，解决了因分量转化而导致磁异常模量计算不准确的问题，也避免了引入由观测点磁场方向不准确带来的系统误差，提高了地质体的总磁化方向的准确性。
实施例二
本实施例的技术方案是对上述实施例的应用举例说明。图2是根据本发明实施例提供的一种地质体的总磁化方向的确定方法的流程图。本实施例为结合上述实施例对本发明进行说明，建立理论球体模型，球体半径5km，中心埋深6km，总磁化强度大小为10A/m，磁偏角为90°(以正北方向为0°起算)，磁倾角为70°，记为(90°,70°)。正常地磁场方向固定：偏角0°，倾角60°。当飞行高度为100m，测量点距1km×1km时，正演模拟的航空实测磁三分量ΔBx，ΔBy和ΔBz如图3(a)(b)(c)所示。针对这个模型实例，对本发明进行详细说明。
a、由实测三分量磁数据计算磁异常模量Ta：
计算的磁异常模量如图4所示，其中，ΔBx，ΔBy和ΔBz分别是第一方向x、第二方向y、第三方向z的实测磁三分量数据。
b、第二参考磁化方向(d,i)由磁偏角d和磁倾角i确定。
预设磁偏角在0°到360°之间，预设磁倾角在-90°到90°之间，按分辨率1°取值做磁化方向网格，网格横坐标为磁偏角d，网格纵坐标为磁倾角i。磁化方向网格如图5所示，且其为网格局部放大图。
c、按第二参考磁化方向的网格扫描所有假定的第二参考磁化方向(d,i)，在波数域实现ΔBz分量的化极计算，得到化极数据ΔB0 z⊥，初始波数域转换因子为：
其中，ωx和ωy分别为x方向(正北方向)和y方向(正东方向)的波数， 为总磁化方向的单位矢量，即/>其中d、i为假定的磁化方向的预设磁偏角、预设磁倾角，即第二参考磁化方向，预设磁偏角在0°到360°之间，预设磁倾角在-90°到90°之间。与常规的总磁场异常ΔT的转换因子不同，该转换因子不需要输入环境地磁场的方向，即不受总磁场方向的影响。
d、计算如图4所示的磁异常模量Ta和第二参考磁化方向(d,i)下的ΔB0 z⊥的第二相关系数
第二相关系数分布图如图6所示。
e、第二相关系数中最大的相关系数对应的磁化方向(d,i)＝(90，80)，即(90，80)为计算的地质体的总磁化方向的初始磁化方向(d0,i0)＝(90，80)。
f、对上一步计算初始磁倾角i进行修正。用垂直磁化的代替Ta计算其与ΔBz⊥的第一相关系数，采用迭代法实现代替的目的以保证算法的稳定性。具体方法是：利用上一步计算的初始磁化方向求垂直磁化的磁三分量/> 和/>波数域转换因子为：
g、通过垂直磁化的磁三分量和/>计算垂直磁化的/>
h、固定磁偏角为初始磁偏角d0，遍历所有假定的预设磁倾角i(-90，90°)，在波数域实现ΔBz分量的化极计算，得到化极数据ΔBz⊥；
i、计算经过垂直磁化后的磁模量和假定的第一参考磁化方向(d0,i)的ΔBz⊥的第一相关系数/>
第一相关系数分布图如图7所示。
j、第一相关系数中最大相关系数对应的磁化方向(d0,i1)即为第1次迭代所求的地质体的总磁化方向；
k、重复步骤f到j，实现迭代计算(d0,i)，直到前后两次计算的磁倾角i的差值足够小，停止迭代，输出(d0,i)，对于本实施例中(d0,i)＝(90°,71°)。计算的磁化方向与理论磁化方向(90°,70°)的误差很小。
本发明实施例的技术方案，获取目标地质体的实测磁三分量数据，然后确定当前次处理周期使用的当前磁化方向，基于第一参考磁化方向对实测磁三分量数据中第三方向磁分量数据进行化极处理获得化极磁异常数据，基于当前磁化方向和实测磁三分量数据确定当前磁模量，并基于当前磁模量与化极磁异常数据确定下一磁化方向，进一步判断第一参数是否满足预设条件，第一参数为当前磁倾角与下一磁倾角的差值的绝对值，预设条件为第一参数小于或等于预设值，从而确定是否对磁化方向做进一步的更新以准确确定目标地质体的总磁化方向，解决了因分量转化而导致磁异常模量计算不准确的问题，也避免了引入由观测点磁场方向不准确带来的系统误差，提高了地质体的总磁化方向的准确性。
实施例三
图8为本发明实施例提供的一种地质体的总磁化方向的确定装置的结构示意图，本实施例可适用于对地质体的总磁化方向进行确定的情况，该地质体的总磁化方向的确定装置可以采用硬件和/或软件的形式实现，该地质体的总磁化方向的确定装置可配置于任何具有网络通信功能的电子设备中。如图8所示，该装置包括：
第一数据确定模块210，用于获取目标地质体的实测磁三分量数据，所述实测磁三分量数据包括第一方向磁分量数据、第二方向磁分量数据和第三方向磁分量数据，所述第一方向为实测时正北方向，第二方向为实测时正东方向，第三方向为实测时垂直向下的方向；
磁化方向确定模块220，用于确定当前次处理周期使用的当前磁化方向，所述当前磁化方向包括初始磁偏角和当前磁倾角，所述初始磁偏角为初始磁化方向包含的一个方位角度，所述初始磁化方向为初次根据所述实测磁三分量数据所确定的磁化方向；
第二数据确定模块230，用于基于第一参考磁化方向对所述实测磁三分量数据中第三方向磁分量数据进行化极处理获得化极磁异常数据，所述第一参考磁化方向包括所述初始磁偏角和预设磁倾角；
信息确定模块240，用于基于所述当前磁化方向和所述实测磁三分量数据确定当前磁模量，并基于所述当前磁模量与所述化极磁异常数据确定下一磁化方向，以供进入下一次处理周期加载使用所述下一磁化方向，所述下一磁化方向包括初始磁偏角和下一磁倾角；
判断模块250，用于确定第一参数，判断第一参数是否满足预设条件，所述第一参数为所述当前磁倾角与下一磁倾角的差值的绝对值，所述预设条件为所述第一参数小于或等于预设值；
第一条件确定模块260，用于若所述第一参数不满足所述预设条件，则将所述下一磁化方向作为当前磁化方向，并基于所述当前磁化方向和所述实测磁三分量数据确定当前磁模量，以及基于所述化极磁异常数据与所述当前磁模量确定下一磁化方向，直至所述第一参数满足预设条件；
第二条件确定模块270，用于若所述第一参数满足所述预设条件，则将下一磁化方向作为所述目标地质体的总磁化方向。
可选的，磁化方向确定模块，用于：
若当前次处理周期为第一个处理周期，则当前磁化方向为初始磁化方向，所述初始磁化方向包括初始磁偏角和初始磁倾角；
若当前次处理周期不为第一个处理周期，则当前磁化方向为上一个处理周期确定的磁化方向。
可选的，磁化方向确定模块包括初始磁化方向确定单元，包括：
第一信息确定单元，用于基于所述实测磁三分量数据确定磁异常模量，并基于第二参考磁化方向对所述第三方向磁分量数据进行化极处理获得初始化极磁异常数据，所述第二参考磁化方向包括预设磁偏角和预设磁倾角；
第二信息确定单元，用于基于所述磁异常模量和所述初始化极磁异常数据确定初始磁化方向，所述初始磁化方向包括初始磁偏角和初始磁倾角。
可选的，第一信息确定单元包括初始化极磁异常数据确定单元，用于：
基于所述第二参考磁化方向确定初始波数域转换因子；
基于所述初始波数域转换因子，对所述第三方向磁分量数据进行化极处理获得初始化极磁异常数据。
可选的，第二数据确定模块，用于：
基于所述第一参考磁化方向确定波数域转换因子；
基于所述波数域转换因子，对所述第三方向磁分量数据进行化极处理获得所述化极磁异常数据。
可选的，信息确定模块包括当前磁模量确定单元，用于：
基于所述当前磁化方向对所述实测磁三分量数据进行垂直磁化处理获得所述当前磁模量。
可选的，信息确定模块包括下一磁化方向确定单元，用于：
确定第一关联关系，所述第一关联关系为第一相关系数、所述化极磁异常数据和所述当前磁模量之间的相关关系；
基于所述化极磁异常数据、所述当前磁模量和所述第一关联关系确定当前相关系数，并将所述当前相关系数对应的磁化方向作为所述下一磁化方向，所述当前相关系数为所述第一相关系数中的最大值。
本发明实施例所提供的地质体的总磁化方向的确定装置可执行本发明任意实施例所提供的地质体的总磁化方向的确定方法，具备执行方法相应的功能模块和有益效果。
实施例四
根据本公开的实施例，本公开还提供了一种电子设备、一种可读存储介质和一种计算机程序产品。
图9示出了可以用来实现本发明实施例的一种地质体的总磁化方向的确定方法的电子设备的结构示意图。电子设备旨在表示各种形式的数字计算机，诸如，膝上型计算机、台式计算机、工作台、个人数字助理、服务器、刀片式服务器、大型计算机、和其它适合的计算机。电子设备还可以表示各种形式的移动装置，诸如，个人数字处理、蜂窝电话、智能电话、可穿戴设备(如头盔、眼镜、手表等)和其它类似的计算装置。本文所示的部件、它们的连接和关系、以及它们的功能仅仅作为示例，并且不意在限制本文中描述的和/或者要求的本发明的实现。
如图9所示，电子设备10包括至少一个处理器11，以及与至少一个处理器11通信连接的存储器，如只读存储器(ROM)12、随机访问存储器(RAM)13等，其中，存储器存储有可被至少一个处理器执行的计算机程序，处理器11可以根据存储在只读存储器(ROM)12中的计算机程序或者从存储单元18加载到随机访问存储器(RAM)13中的计算机程序，来执行各种适当的动作和处理。在RAM 13中，还可存储电子设备10操作所需的各种程序和数据。处理器11、ROM 12以及RAM 13通过总线14彼此相连。输入/输出(I/O)接口15也连接至总线14。
电子设备10中的多个部件连接至I/O接口15，包括：输入单元16，例如键盘、鼠标等；输出单元17，例如各种类型的显示器、扬声器等；存储单元18，例如磁盘、光盘等；以及通信单元19，例如网卡、调制解调器、无线通信收发机等。通信单元19允许电子设备10通过诸如因特网的计算机网络和/或各种电信网络与其他设备交换信息/数据。
处理器11可以是各种具有处理和计算能力的通用和/或专用处理组件。处理器11的一些示例包括但不限于中央处理单元(CPU)、图形处理单元(GPU)、各种专用的人工智能(AI)计算芯片、各种运行机器学习模型算法的处理器、数字信号处理器(DSP)、以及任何适当的处理器、控制器、微控制器等。处理器11执行上文所描述的各个方法和处理，例如一种地质体的总磁化方向的确定方法。
在一些实施例中，一种地质体的总磁化方向的确定方法可被实现为计算机程序，其被有形地包含于计算机可读存储介质，例如存储单元18。在一些实施例中，计算机程序的部分或者全部可以经由ROM 12和/或通信单元19而被载入和/或安装到电子设备10上。当计算机程序加载到RAM 13并由处理器11执行时，可以执行上文描述的一种地质体的总磁化方向的确定方法的一个或多个步骤。备选地，在其他实施例中，处理器11可以通过其他任何适当的方式(例如，借助于固件)而被配置为执行一种地质体的总磁化方向的确定方法。
本文中以上描述的系统和技术的各种实施方式可以在数字电子电路系统、集成电路系统、现场可编程门阵列(FPGA)、专用集成电路(ASIC)、专用标准产品(ASSP)、芯片上系统的系统(SOC)、复杂可编程逻辑设备(CPLD)、计算机硬件、固件、软件、和/或它们的组合中实现。这些各种实施方式可以包括：实施在一个或者多个计算机程序中，该一个或者多个计算机程序可在包括至少一个可编程处理器的可编程系统上执行和/或解释，该可编程处理器可以是专用或者通用可编程处理器，可以从存储系统、至少一个输入装置、和至少一个输出装置接收数据和指令，并且将数据和指令传输至该存储系统、该至少一个输入装置、和该至少一个输出装置。
用于实施本发明的方法的计算机程序可以采用一个或多个编程语言的任何组合来编写。这些计算机程序可以提供给通用计算机、专用计算机或其他可编程数据处理装置的处理器，使得计算机程序当由处理器执行时使流程图和/或框图中所规定的功能/操作被实施。计算机程序可以完全在机器上执行、部分地在机器上执行，作为独立软件包部分地在机器上执行且部分地在远程机器上执行或完全在远程机器或服务器上执行。
在本发明的上下文中，计算机可读存储介质可以是有形的介质，其可以包含或存储以供指令执行系统、装置或设备使用或与指令执行系统、装置或设备结合地使用的计算机程序。计算机可读存储介质可以包括但不限于电子的、磁性的、光学的、电磁的、红外的、或半导体系统、装置或设备，或者上述内容的任何合适组合。备选地，计算机可读存储介质可以是机器可读信号介质。机器可读存储介质的更具体示例会包括基于一个或多个线的电气连接、便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦除可编程只读存储器(EPROM或快闪存储器)、光纤、便捷式紧凑盘只读存储器(CD-ROM)、光学储存设备、磁储存设备、或上述内容的任何合适组合。
为了提供与用户的交互，可以在电子设备上实施此处描述的系统和技术，该电子设备具有：用于向用户显示信息的显示装置(例如，CRT(阴极射线管)或者LCD(液晶显示器)监视器)；以及键盘和指向装置(例如，鼠标或者轨迹球)，用户可以通过该键盘和该指向装置来将输入提供给电子设备。其它种类的装置还可以用于提供与用户的交互；例如，提供给用户的反馈可以是任何形式的传感反馈(例如，视觉反馈、听觉反馈、或者触觉反馈)；并且可以用任何形式(包括声输入、语音输入或者、触觉输入)来接收来自用户的输入。
可以将此处描述的系统和技术实施在包括后台部件的计算系统(例如，作为数据服务器)、或者包括中间件部件的计算系统(例如，应用服务器)、或者包括前端部件的计算系统(例如，具有图形用户界面或者网络浏览器的用户计算机，用户可以通过该图形用户界面或者该网络浏览器来与此处描述的系统和技术的实施方式交互)、或者包括这种后台部件、中间件部件、或者前端部件的任何组合的计算系统中。可以通过任何形式或者介质的数字数据通信(例如，通信网络)来将系统的部件相互连接。通信网络的示例包括：局域网(LAN)、广域网(WAN)、区块链网络和互联网。
计算系统可以包括客户端和服务器。客户端和服务器一般远离彼此并且通常通过通信网络进行交互。通过在相应的计算机上运行并且彼此具有客户端-服务器关系的计算机程序来产生客户端和服务器的关系。服务器可以是云服务器，又称为云计算服务器或云主机，是云计算服务体系中的一项主机产品，以解决了传统物理主机与VPS服务中，存在的管理难度大，业务扩展性弱的缺陷。
应该理解，可以使用上面所示的各种形式的流程，重新排序、增加或删除步骤。例如，本发明中记载的各步骤可以并行地执行也可以顺序地执行也可以不同的次序执行，只要能够实现本发明的技术方案所期望的结果，本文在此不进行限制。
上述具体实施方式，并不构成对本发明保护范围的限制。本领域技术人员应该明白的是，根据设计要求和其他因素，可以进行各种修改、组合、子组合和替代。任何在本发明的精神和原则之内所作的修改、等同替换和改进等，均应包含在本发明保护范围之内。

Claims (10)

1.一种地质体的总磁化方向的确定方法，其特征在于，包括：
获取目标地质体的实测磁三分量数据，所述实测磁三分量数据包括第一方向磁分量数据、第二方向磁分量数据和第三方向磁分量数据，所述第一方向为实测时正北方向，所述第二方向为实测时正东方向，所述第三方向为实测时垂直向下的方向；
确定当前次处理周期使用的当前磁化方向，所述当前磁化方向包括初始磁偏角和当前磁倾角，所述初始磁偏角为初始磁化方向包含的一个方位角度，所述初始磁化方向为初次根据所述实测磁三分量数据所确定的磁化方向；
基于第一参考磁化方向对所述实测磁三分量数据中第三方向磁分量数据进行化极处理获得化极磁异常数据，所述第一参考磁化方向包括所述初始磁偏角和预设磁倾角；
基于所述当前磁化方向和所述实测磁三分量数据确定当前磁模量，并基于所述当前磁模量与所述化极磁异常数据确定下一磁化方向，以供进入下一次处理周期加载使用所述下一磁化方向，所述下一磁化方向包括初始磁偏角和下一磁倾角；
确定第一参数，判断所述第一参数是否满足预设条件，所述第一参数为所述当前磁倾角与下一磁倾角的差值的绝对值，所述预设条件为所述第一参数小于或等于预设值；
若所述第一参数不满足所述预设条件，则将所述下一磁化方向作为当前磁化方向，并基于所述当前磁化方向和所述实测磁三分量数据确定当前磁模量，以及基于所述化极磁异常数据与所述当前磁模量确定下一磁化方向，直至所述第一参数满足预设条件；
若所述第一参数满足所述预设条件，则将下一磁化方向作为所述目标地质体的总磁化方向。
2.根据权利要求1所述的方法，其特征在于，确定当前次处理周期使用的当前磁化方向，包括：
若当前次处理周期为第一个处理周期，则当前磁化方向为初始磁化方向，所述初始磁化方向包括初始磁偏角和初始磁倾角；
若当前次处理周期不为第一个处理周期，则当前磁化方向为上一个处理周期确定的磁化方向。
3.根据权利要求1或2所述的方法，其特征在于，所述初始磁化方向的确定过程，包括：
基于所述实测磁三分量数据确定磁异常模量，并基于第二参考磁化方向对所述第三方向磁分量数据进行化极处理获得初始化极磁异常数据，所述第二参考磁化方向包括预设磁偏角和预设磁倾角；
基于所述磁异常模量和所述初始化极磁异常数据确定初始磁化方向，所述初始磁化方向包括初始磁偏角和初始磁倾角。
4.根据权利要求2所述的方法，其特征在于，基于第二参考磁化方向对所述第三方向磁分量数据进行化极处理获得初始化极磁异常数据，包括：
基于所述第二参考磁化方向确定初始波数域转换因子；
基于所述初始波数域转换因子，对所述第三方向磁分量数据进行化极处理获得初始化极磁异常数据。
5.根据权利要求1所述的方法，其特征在于，基于第一参考磁化方向对所述实测磁三分量数据中第三方向磁分量数据进行化极处理获得化极磁异常数据，包括：
基于所述第一参考磁化方向确定波数域转换因子；
基于所述波数域转换因子，对所述第三方向磁分量数据进行化极处理获得所述化极磁异常数据。
6.根据权利要求1所述的方法，其特征在于，基于所述当前磁化方向和所述实测磁三分量数据确定当前磁模量，包括：
基于所述当前磁化方向对所述实测磁三分量数据进行垂直磁化处理获得所述当前磁模量。
7.根据权利要求1所述的方法，其特征在于，基于所述当前磁模量与所述化极磁异常数据确定下一磁化方向，包括：
确定第一关联关系，所述第一关联关系为第一相关系数、所述化极磁异常数据和所述当前磁模量之间的相关关系；
基于所述化极磁异常数据、所述当前磁模量和所述第一关联关系确定当前相关系数，并将所述当前相关系数对应的磁化方向作为所述下一磁化方向，所述当前相关系数为所述第一相关系数中的最大值。
8.一种地质体的总磁化方向的确定装置，其特征在于，包括：
第一数据确定模块，用于获取目标地质体的实测磁三分量数据，所述实测磁三分量数据包括第一方向磁分量数据、第二方向磁分量数据和第三方向磁分量数据，所述第一方向为实测时正北方向，所述第二方向为实测时正东方向，所述第三方向为实测时垂直向下的方向；
磁化方向确定模块，用于确定当前次处理周期使用的当前磁化方向，所述当前磁化方向包括初始磁偏角和当前磁倾角，所述初始磁偏角为初始磁化方向包含的一个方位角度，所述初始磁化方向为初次根据所述实测磁三分量数据所确定的磁化方向；
第二数据确定模块，用于基于第一参考磁化方向对所述实测磁三分量数据中第三方向磁分量数据进行化极处理获得化极磁异常数据，所述第一参考磁化方向包括所述初始磁偏角和预设磁倾角；
信息确定模块，用于基于所述当前磁化方向和所述实测磁三分量数据确定当前磁模量，并基于所述当前磁模量与所述化极磁异常数据确定下一磁化方向，以供进入下一次处理周期加载使用所述下一磁化方向，所述下一磁化方向包括初始磁偏角和下一磁倾角；
判断模块，用于确定第一参数，判断所述第一参数是否满足预设条件，所述第一参数为所述当前磁倾角与下一磁倾角的差值的绝对值，所述预设条件为所述第一参数小于或等于预设值；
第一条件确定模块，用于若所述第一参数不满足所述预设条件，则将所述下一磁化方向作为当前磁化方向，并基于所述当前磁化方向和所述实测磁三分量数据确定当前磁模量，以及基于所述化极磁异常数据与所述当前磁模量确定下一磁化方向，直至所述第一参数满足预设条件；
第二条件确定模块，用于若所述第一参数满足所述预设条件，则将下一磁化方向作为所述目标地质体的总磁化方向。
9.一种电子设备，其特征在于，所述电子设备包括：
至少一个处理器；以及
与所述至少一个处理器通信连接的存储器；其中，
所述存储器存储有可被所述至少一个处理器执行的计算机程序，所述计算机程序被所述至少一个处理器执行，以使所述至少一个处理器能够执行权利要求1-7中任一项所述的一种地质体的总磁化方向的确定方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于使处理器执行时实现权利要求1-7中任一项所述的一种地质体的总磁化方向的确定方法。
CN202311426553.9A 2023-10-31 2023-10-31 一种地质体的总磁化方向的确定方法、装置、设备和介质 Active CN117434614B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311426553.9A CN117434614B (zh) 2023-10-31 2023-10-31 一种地质体的总磁化方向的确定方法、装置、设备和介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311426553.9A CN117434614B (zh) 2023-10-31 2023-10-31 一种地质体的总磁化方向的确定方法、装置、设备和介质

Publications (2)

Publication Number Publication Date
CN117434614A true CN117434614A (zh) 2024-01-23
CN117434614B CN117434614B (zh) 2024-08-02

Family

ID=89554900

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311426553.9A Active CN117434614B (zh) 2023-10-31 2023-10-31 一种地质体的总磁化方向的确定方法、装置、设备和介质

Country Status (1)

Country Link
CN (1) CN117434614B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005001494A2 (fr) * 2003-06-24 2005-01-06 Magnetic@Soft Procede de calcul du champ magnetique genere par un systeme materiel et dispositif de mise en oeuvre d’un tel procede de calcul
WO2013134410A1 (en) * 2012-03-07 2013-09-12 Crocus Technology Inc. Magnetic logic units configured to measure magnetic field direction
CN106405664A (zh) * 2016-08-25 2017-02-15 中国科学院地质与地球物理研究所 一种磁异常化极方法
CN114021408A (zh) * 2021-11-05 2022-02-08 中南大学 二维强磁场数值模拟方法、装置、设备及介质
CN114296143A (zh) * 2021-12-30 2022-04-08 中国地质大学（武汉） 基于广义多重相关的地下磁性体总磁化方向估计方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005001494A2 (fr) * 2003-06-24 2005-01-06 Magnetic@Soft Procede de calcul du champ magnetique genere par un systeme materiel et dispositif de mise en oeuvre d’un tel procede de calcul
WO2013134410A1 (en) * 2012-03-07 2013-09-12 Crocus Technology Inc. Magnetic logic units configured to measure magnetic field direction
CN106405664A (zh) * 2016-08-25 2017-02-15 中国科学院地质与地球物理研究所 一种磁异常化极方法
CN114021408A (zh) * 2021-11-05 2022-02-08 中南大学 二维强磁场数值模拟方法、装置、设备及介质
CN114296143A (zh) * 2021-12-30 2022-04-08 中国地质大学（武汉） 基于广义多重相关的地下磁性体总磁化方向估计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李金朋;张英堂;范红波;李志宁;尹刚;武禹陶;: "剩磁条件下磁性体磁化方向相关估计方法", 大地测量与地球动力学, no. 07, 15 July 2017 (2017-07-15) *

Also Published As

Publication number Publication date
CN117434614B (zh) 2024-08-02

Similar Documents

Publication Publication Date Title
US20190242952A1 (en) Compensating for Distortion in an Electromagnetic Tracking System
CN110146839A (zh) 一种移动平台磁梯度张量系统校正方法
CN109581383A (zh) 各向异性海面雷达后向散射模拟方法及系统
Zhao et al. Robust adaptive heading tracking fusion for polarization compass with uncertain dynamics and external disturbances
CN117434614B (zh) 一种地质体的总磁化方向的确定方法、装置、设备和介质
CN116368403A (zh) 使用交叉偶极声波数据表征方位各向异性的方法
Pace The solution of the Einstein’s equations in the vacuum region surrounding a spherically symmetric mass distribution
CN104316960B (zh) 一种基于vsp的储层裂缝预测方法及系统
CN115805593B (zh) 力传感器安装信息确定方法、装置、设备和介质
CN115984501A (zh) 一种地质三维模型建立方法、装置、设备及存储介质
CN114972511A (zh) 标定物位姿的确定方法、装置、电子设备及存储介质
CN115032704A (zh) 一种磁异常探测方法、系统、介质、设备及终端
CN111596346B (zh) 弹性波速度反演方法和装置
CN108983302A (zh) 一种极区卫星重力梯度张量无θ奇异性详细计算模型
CN116380056B (zh) 惯性定位方法、装置、电子设备及存储介质
CN116753899A (zh) 一种平面结构的空间姿态检测方法、装置、设备及介质
CN116358606B (zh) 惯导系统的初始粗对准方法、装置、设备及介质
CN117592381B (zh) 一种大气波导参数反演模型训练方法、装置、设备及介质
CN110207699A (zh) 一种定位方法和装置
CN116736374B (zh) 多方位地震数据的融合方法、装置、设备和介质
CN115861435A (zh) 雷达外参评估方法、装置、电子设备及存储介质
CN117908109A (zh) 一种盐岩地层速度的确定方法、装置、设备及介质
Jianjun et al. Recent Advances in the Geodesy Data Processing
CN118803572A (zh) Mr指纹定位方法及装置、电子设备和存储介质
CN118501940A (zh) 一种致密气藏储层含气饱和度预测方法、装置、电子设备及存储介质

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