[发明专利]一种井中磁力异常估计有效磁化倾角的方法

权利要求书

1.一种井中磁力异常估计有效磁化倾角的方法，其特征在于，包括以下步骤：

S1、基于井中磁异常三分量包括垂直分量Z_a、水平分量H_ax和水平分量H_ay数据计算得到磁异常模量T_a；

S2、对Z_a、H_ax和H_ay三个分量进行傅里叶正变换分别得到相应的频谱；

S3、设置初始有效磁化倾角度数；

S4、基于S3设置的有效磁化倾角度数以及井与异常体的相对位置关系计算出垂直磁化方向转换因子以及水平磁化方向转换因子；

S5、Z_a频谱乘以垂直磁化方向转换因子，H_ax和H_ay的频谱分别乘以水平磁化方向转换因子；

S6、将乘以转换因子的频谱分别进行傅里叶逆变换得到磁化方向转到垂直方向的Z_a数据以及磁化方向转到水平方向的H_ax和H_ay数据；

S7、将S6中转换后的Z_a、H_ax和H_ay分别与S1中磁异常模量T_a进行相关性计算，记录得到的相关性数值；

S8、将有效磁化倾角增加n°，0°≤n°≤360°，循环S4至S8；

S9、当有效磁化倾角增加到指定的度数，停止循环，基于记录的相关性数值绘制出Z_a、H_ax和H_ay的三条相关性曲线，寻找Z_a相关性曲线的最小值点以及H_ax和H_ay相关性曲线的最大值点，得到的最值点对应的度数即为估计的有效磁化倾角；

S4中频率域转换因子与井和异常体的相对位置有关，具体关系为：在建立的坐标系下，所研究剖面中，当井位于异常体正方向即在对应横坐标轴正方向时，磁化方向转换因子为：

α₂＝cosI_s2γ₂＝sinI_s2α₁＝cosI_s1γ₁＝sinI_s1

其中，w表示磁化方向转化因子，I_s2为转换后磁化倾角，I_s1为原磁化倾角；

当井位于异常体负方向即横坐标轴负方向时，转换因子为：

α₂＝cosI_s2γ₂＝sinI_s2α₁＝cosI_s1γ₁＝sinI_s1

其中，I_s2为转换后磁化倾角，I_s1为原磁化倾角，α₁，γ₁为原磁化方向余弦，α₂，γ₂为新磁化方向余弦；

当I_s2为90°时，转换因子为垂直磁化方向转换因子；当I_s2为0°时，转换因子为水平磁化方向转换因子。

2.根据权利要求1所述的一种井中磁力异常估计有效磁化倾角的方法，其特征在于，S1中磁异常模量T_a的计算公式如下：

其中，Ta表示磁异常模量，Z_a表示垂直分量，H_ax表示水平分量，H_ay表示水平分量。

3.根据权利要求1所述的一种井中磁力异常估计有效磁化倾角的方法，其特征在于，S3中初始有效磁化倾角度数设置为x°，其中0°≤x°≤360°。