【技术实现步骤摘要】
一种基于网格点格架的磁张量异常空间域快速正演算法
本专利技术涉及重磁勘探的正演计算，特别是一种基于网格点格架的磁张量异常空间域快速正演算法。
技术介绍
重磁勘探的正演计算中，常用等距网格线将模型剖分成若干个长方体单元，然后计算每一个长方体单元对观测点的异常，再将所有长方体单元对观测点的异常求和，即为整个模型体对观测点的异常。这其中存在大量重复的计算，导致数值模拟及反演效率大大降低。姚长利，郝天珧，管志宁，张聿文在地球物理学报，2003，46(2)：252-258.发表的重磁遗传算法三维反演中高速计算及有效存储方法技术[J].中利用三维密度模型的正演公式说明剖分模型的正演情况，提出基于几何格架等效压缩存储技术，利用平移等效性简化磁异常类型中的倾斜磁化问题。该方法先计算某个长方体单元对某一层所有网格点的异常，然后存储再调用，提高了计算效率，但由于几何格架存在重复的计算，因此其计算效率可以进一步提高。
技术实现思路
本专利技术的专利技术目的在于：针对现有技术存在的几何格架存在重复计算的问题，提供一种基于...
【技术保护点】
1.一种基于网格点格架的磁张量异常空间域快速正演算法，其特征在于，包括如下步骤：/nS100根据实际里程建立模型体，以设置的网格间距将模型体剖分成网格，即模型体由多个长方体单元组合而成，在网格中异常体的位置相应设置异常体的剩余磁化率；/nS200分别计算模型空间内所有网格点对观测点P1、P2、P3、P4的磁张量异常的格架函数，然后存储为格架函数数据库；其中，P1、P2、P3、P4是模型体上表面的四个顶点；/nS300将模型体上表面的网格点作为观测点，计算整个模型体对一个观测点的磁张量异常；/nS310任一长方体单元有8个顶点，每个顶点与所述观测点有相对位置关系，根据平移等效...
【技术特征摘要】
1.一种基于网格点格架的磁张量异常空间域快速正演算法，其特征在于，包括如下步骤：
S100根据实际里程建立模型体，以设置的网格间距将模型体剖分成网格，即模型体由多个长方体单元组合而成，在网格中异常体的位置相应设置异常体的剩余磁化率；
S200分别计算模型空间内所有网格点对观测点P1、P2、P3、P4的磁张量异常的格架函数，然后存储为格架函数数据库；其中，P1、P2、P3、P4是模型体上表面的四个顶点；
S300将模型体上表面的网格点作为观测点，计算整个模型体对一个观测点的磁张量异常；
S310任一长方体单元有8个顶点，每个顶点与所述观测点有相对位置关系，根据平移等效性，分别调用格架函数数据库中所述相对位置关系对应的格架函数，将8个格架函数求和后乘以物性参数得到该长方体单元对所述观测点的磁张量异常；
S320重复步骤S310，直到得到所有长方体单元对该观测点的磁张量异常；
S330将所有长方体单元对该观测点的磁张量异常求和，即为整个模型体对该观测点的磁张量异常；
S400重复步骤S300，直到得出整个模型体对工作区平面各观测点的磁张量异常，即全张量磁梯度数据Bxx、Bxy、Bxz、Byy、Byz、Bzz。
2.根据权利要求1所述的算法，其特征在于，
所述步骤S100中，异常体的剩余磁化率为：去背景场的观测数据异常，即异常体的物性差。
3.根据权利要求2所述的算法，其特征在于，
所述步骤S200中，所述格架函数数据库包括所有网格点对观测点P1、P2、P3、P4的磁张量异常的格架函数集合MagneticP1、MagneticP2、MagneticP3、MagneticP4；观测点P1坐标为(1,1)，观测点P2坐标为(NX+1,1)，观测点P3坐标为(NX+1,NY+1)，观测点P4坐标为(1,NY+1),NX、NY分别表示模型体X向、Y向长方体单元个数。
4.根据权利要求3所述的算法，其特征在于，
所述步骤S310中，调用格架函数数据库中所述相对位置关系对应的格架函数，具体包括：
S311判断观测点和长方体单元的相对位置关系类型；
S...
【专利技术属性】
技术研发人员：王科，张志厚，寇川，范祥泰，赵思为，路润琪，尹小康，姚禹，
申请(专利权)人：中铁二院工程集团有限责任公司，西南交通大学，
类型：发明
国别省市：四川;51
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。