[发明专利]钻孔瞬变电磁一维反演电性信息立体成像方法

权利要求书

1.一种钻孔瞬变电磁反演电性信息立体成像方法，其特征在于，包括以下步骤：

步骤一，对观测得到的钻孔瞬变电磁垂直分量进行一维反演得到径向探测距离r和电阻率，绘制反演电阻率断面图，从中求出各电阻率异常区域所占测点数平均值；

步骤二，以钻孔瞬变电磁水平分量中每个测点为插值窗口中心，步骤一中的测点数平均值作为水平方向插值窗口大小，提取异常场得到水平分量异常曲线，对异常场的异常数据进行正规化后提取特征点及特征点对应的正规化点号；

步骤三，根据正规化点号的整体异常分布确定每个异常场对应的水平分量异常曲线：正弦形态或反正弦形态；

步骤四，根据分类结果求取工具面角度：

步骤4.1，对同一测点同一测道的水平分量异常曲线类型进行组合，按照曲线分类图所示依据自动划分对应测道电性信息所处的象限；所述曲线分类图中，第一象限为X反正弦形态，Y反正弦形态，第二象限为X正弦形态，Y反正弦形态，第三象限为X正弦形态，Y正弦形态，第四象限为X反正弦形态，Y正弦形态；

步骤4.2，由如下公式计算所有测点所有测道的电性信息在对应象限中的偏转角度：

式中为第n个测道X分量的异常场值，为第n个测道Y分量的异常场值；

步骤4.3，根据对应测点测道电性信息对应的象限，求出所有电性信息对应的工具面角α：a)异常体在第一象限：α＝6；b)异常体在第二象限：α＝π-θ；c)异常体在第三象限：α＝π+θ；d)异常体在第四象限：α＝2π-θ；

步骤五，步骤四得到的工具面角与测道一一对应，通过建立径向探测距离-电阻率与测道的映射关系，进而得到径向探测距离-电阻率与工具面角的对应关系，实现将径向探测距离-电阻率分配给对应工具面角：

步骤5.1，引入径向探测距离计算公式

式中，r为径向探测距离，ρ为电阻率，t为测道，C为深度系数；

步骤5.2，计算大量一维正演模型，对比不同电阻率参数情况下的层状模型分界面的计算深度与实际深度之间的差别，绘制出深度系数C与电阻率ρ的关系图；

步骤5.3，在对数坐标系下，深度系数与电阻率呈现线性关系，通过拟合可以得到直线方程为

log10(C)＝-0.50009365×log10(ρ)-0.00137449

步骤5.4，根据映射关系公式逐层推导出层状地层情况下径向探测距离与测道之间的计算公式，建立二者之间的映射关系，进而建立反演径向探测距离与每一测道对应的工具面角的一一对应关系；

具体计算方式如下：

a.在地下第一层有如下关系成立，C₁为第一层深度系数

b.根据如下公式计算每一层对应的标准系数C_i

log10(C_i)＝-0.50009365×log10(ρ_i)-0.00137449

c.根据如下公式计算每一层对应的相对系数Cc_i

A_i＝log10(C_i)/log10(C₁)

d.除第一层外的其他地层真深度与测道可建立如下关系

即在第一层外的其他地层的测道可以按如下公式计算

式中，i为层状地层中的第i层；由径向探测距离和电阻率可反推出对应的测道，再通过插值可以获得实际采样测道对应的径向探测距离和电阻率；

步骤六，由工具面角与对应测点测道的径向探测距离计算反演电性信息对应的平面坐标

式中为第m个测点第n个测道视电阻率对应的X坐标，为第m个测点第n个测道视电阻率对应的Y坐标，为第m个测点第n个测道视电阻率对应的钻孔径向探测距离；

步骤七，以平面坐标分别为X、Y，以钻孔深度为Z，绘制钻孔周围一维反演电性信息立体成像结果图，以将低值区域划分为可能的隐蔽致灾水体。