一种电磁数据的保幅静校正方法
一种电磁数据的保幅静校正方法
【技术领域】
[0001] 本发明涉及地球物理勘探方法，是一种电磁勘探数据，包括大地电磁（MT)、音频大 地电磁（AMT)和连续剖面法大地电磁（CEMP)的数据处理技术。
【背景技术】
[0002] 电磁法勘探中的静态效应是指由于近表存在局部电性不均匀体时，电阻率分界面 上极化电荷的堆集引起电场的畸变，由此产生一个与外电场成正比的附加电场，且与频率 无关。对于大地电磁数据（MT)，表现在单点曲线上，就是对数域电阻率曲线沿纵轴产生平 移，相位曲线基本上不受影响。表现在视电阻率断面图上，就是电阻率值出现直立的陡变 带，俗称"挂面条"现象，而相应的相位拟断面图上就没有这种现象。表现在单频点的视电 阻率平面图上，就是存在很多"畸变点"现象。静态效应的强度可达两三个数量级，在推断 深度时可引起很大误差，并使构造解释复杂化，因此在反演解释之前需进行静校正。
[0003] 目前，静校正的方法很多，为了压制静态效应，不断提高MT的应用效果，国内外很 多专家发表了相关论文，提出了滤波法、理论计算法、瞬变电磁校正法、地表电阻率校正法、 高频电磁函数向上延拓法等。这些方法大部分都是基于二维MT静校正，各种校正方法的原 理和对资料的要求不同，在实际应用中都受到一定程度的限制。由于滤波法易于实现，处理 速度快，在MT/CEMP资料的静校正处理中应用最为广泛。但是经过滤波法处理后，很多测点 的视电阻率曲线形态会被改变，尤其是在山前带、地表岩性突变带等地方表现明显，实践中 静校正效果不好。
【发明内容】
[0004] 本发明针对采用传统的滤波法进行电磁数据静校正后效果不好的缺陷，提供一种 可以克服视电阻率曲线形态会被改变的缺陷、达到保幅的目的电磁数据的保幅静校正方 法。
[0005] 本发明采用以下步骤实现：
[0006] 1)根据工区采集的大地电磁数据，确定一个深度和电阻率稳定的电性层，得到电 磁数据曲线上对应的频段，对频段上电磁数据测线上所有测点的两支视电阻率曲线做平 均，计算平均视电阻率数据<;
[0007] 所述的电磁数据包括大地电磁数据（MT)、音频大地电磁数据（AMT)、连续剖面法 大地电磁数据（CEMP)。
[0008] 所述的平均包括算术平均和几何平均。
[0009] 所述的计算平均视电阻率数据￥是：
[0010]
[0011] 其中，<(/i)、分别为第i个测点第j个频率A的XY和YX模式实测视 电阻率值，n和m为选定的频点号。
[0012] 2)利用以下公式求取各个测点两支曲线的初步校正因子；
[0014] 其中，f和分别为第i个测点XY和YX模式数据的初步校正因子，Au为由（4) 式获得的计算平均值，其它参数的含义与（4)式相同。
[0015] 3)按照以下公式进行初步校正，得到校正后的两支视电阻率曲线数据/^f和
[0016] 再对Af/f分别进行滤波，获得滤波后的视电阻率曲线数据pf和/其它参 数的含义与（4)、（5)式相同。
[0017] 所述的滤波包括二维中值滤波、二维空间滤波、三维中值滤波和三维空间滤波。
[0018] 4)重新选择一个频率段，按照步骤1)至3)分别计算出滤波前和滤波后各个测点 两支曲线在该频段范围内视电阻率的算术或者几何平均值PfT，和
[0021] 其中，nfl和nf2分别为选定的频点号，Af (乃)分别为第si个测点第j 个频率&的XY和YX模式实测视电阻率值，PfX分别为计算的滤波前XY和YX模式 的视电阻率算术或者几何平均值，和P；f分别为计算的滤波后XY和YX模式的视电阻 率算术或者几何平均值。
[0022] 5)按照以下公式求取各个测点两支曲线的保幅校正因子：
[0024] 其中，#和〇别为第i个测点的XY和YX曲线的保幅校正因子。
[0025] 6)按照以下公式计算获得最终的保幅校正后两支曲线视电阻率数据和 完成静校正。
[0027] 其中，和就是最终的XY和YX曲线的保幅校正后视电阻率数据。
[0028] 与【背景技术】相比，本发明可以克服应用传统的滤波方法后视电阻率曲线形态会被 改变的缺陷，尤其是在山前带、地表岩性突变带等地方，效果改善明显，达到保幅的目的。
【附图说明】
[0029] 图1为电磁数据保幅静校正方法基本流程图。
【具体实施方式】
[0030] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0031] 以下结合附图和实验实例详细说明本发明。
[0032] 图1是保幅静校正方法实施的流程图。具体方法是：
[0033] 1)根据工区采集的大地电磁数据，确定一个深度和电阻率稳定的电性层，得到电 磁数据曲线上对应的频段，如选取第5号到第9号频率段，则公式（4)中n= 5，m= 9。对 频段上电磁数据测线上所有测点的两支视电阻率曲线做平均，计算平均视电阻率数据;
[0034] 所述的电磁数据包括大地电磁数据（MT)、音频大地电磁数据（AMT)、连续剖面法 大地电磁数据（CEMP)。
[0035] 所述的平均包括算术平均和几何平均。
[0036] 所述的计算平均视电阻率数据Pf是：
[0038] 其中，W(/,)、/f(/,)分别为第i个测点第j个频率fj的XY和YX模式实测视 电阻率值，n和m为选定的频点号，在本例子中分别为5和9。
[0039] 2)利用以下公式求取各个测点两支曲线的初步校正因子；
[0041] 其中，f和f分别为第i个测点XY和YX模式数据的初步校正因子，pf为由（4) 式获得的计算平均值，其它参数的含义与（4)式相同。
[0042] 3)按照以下公式进行初步校正，得到校正后的两支视电阻率曲线数据pf和Pf:
[0044] 再对Pf分别进行滤波，如在本例子中，采用二维中值滤波，中值滤波窗口选 择为（5, 3)，获得滤波后的视电阻率曲线数据pjf和其它参数的含义与（4)、（5)式相 同。
[0045] 所述的滤波包括二维中值滤波、二维空间滤波、三维中值滤波和三维空间滤波。
[0046] 4)重新选择一个频率段，如第7号频率到第13号频率段，即公式（7)和公式⑶ 中的nfl= 7，nf2 = 13,按照步骤1)至3)分别计算出滤波前和滤波后各个测点两支曲线 在该频段范围内视电阻率的算术或者几何平均值:
[0049] 其中，nfl和nf2分别为选定的频点号，在本例子中分别为7和13， Af(/,_)、Af(/;)分别为第si个测点第j个频率&的XY和YX模式实测视电阻率值， 、y<}T分别为计算的滤波前XY和YX模式的视电阻率算术或者几何平均值，和 分别为计算的滤波后XY和YX模式的视电阻率算术或者几何平均值。
[0050] 5)按照以公式（9)求取各个测点两支曲线的保幅校正因子：
[0052] 其中，句厂和片分别为第i个测点的XY和YX曲线的保幅校正因子。
[0053] 6)按照公式（10)计算获得最终的保幅校正后两支曲线视电阻率数据和， 完成保幅静校正。
[0055] 其中，和就是最终的XY和YX曲线的保幅校正后视电阻率数据。
[0056] 以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详 细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保 护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本 发明的保护范围之内。
【主权项】
1. 一种电磁数据的保幅静校正方法，其特征在于采用以下步骤实现： 1) 根据工区采集的大地电磁数据，确定一个深度和电阻率稳定的电性层，得到电磁数 据曲线上对应的频段，对频段上电磁数据测线上所有测点的两支视电阻率曲线做平均，计 算平均视电阻率数据<; 所述的计算平均视电阻率数据P是：其中，αλτ(Λ)、AnC/;)分别为第i个测点第j个频率A的XY和YX模式实测视电阻 率值，η和m为选定的频点号； 2) 利用以下公式求取各个测点两支曲线的初步校正因子；其中，和Af分别为第i个测点XY和YX模式数据的初步校正因子，Pf为由（4)式 获得的计算平均值，其它参数的含义与（4)式相同。 3) 按照以下公式进行初步校正，得到校正后的两支视电阻率曲线数据/jf和Pf :(6) 再对Af 分别进行滤波，获得滤波后的视电阻率曲线数据/7=和其它参数的 含义与（4)、（5)式相同。 4) 重新选择一个频率段，按照步骤1)至3)分别计算出滤波前和滤波后各个测点两支 曲线在该频段范围内视电阻率的算术或者几何平均值:(7) (8) 其中，nfl和nf2分别为选定的频点号，Af Pf (/v)分别为第si个测点第j个 频率&的XY和YX模式实测视电阻率值，Pfy、分别为计算的滤波前XY和YX模式的 视电阻率算术或者几何平均值，和/分别为计算的滤波后XY和YX模式的视电阻率 算术或者几何平均值。 5) 按照以下公式求取各个测点两支曲线的保幅校正因子：其中，和#分别为第i个测点的XY和YX曲线的保幅校正因子。 YY 6) 按照以下公式计算获得最终的保幅校正后两支曲线视电阻率数据/和户#，完成 静校正；(10) 其中，/^和就是最终的XY和YX曲线的保幅校正后视电阻率数据。2. 根据权利要求1所述的一种电磁数据的保幅静校正方法，其特征在于步骤1)所述 的电磁数据包括大地电磁数据（MT)、音频大地电磁数据（AMT)、连续剖面法大地电磁数据 (CEMP)〇3. 根据权利要求1所述的一种电磁数据的保幅静校正方法，其特征在于步骤1)所述的 平均包括算术平均和几何平均。4. 根据权利要求1所述的一种电磁数据的保幅静校正方法，其特征在于步骤3)所述的 滤波包括二维中值滤波、二维空间滤波、三维中值滤波和三维空间滤波，其中：μ是土体介质的泊松比。
【专利摘要】一种电磁数据的保幅静校正方法，采用以下步骤实现：1)根据工区采集的大地电磁数据，确定一个深度和电阻率稳定的电性层，得到电磁数据曲线上对应的频段，对频段上电磁数据测线上所有测点的两支视电阻率曲线做平均，计算平均视电阻率数据；2)求取各个测点两支曲线的初步校正因子；3)进行初步校正，得到校正后的两支视电阻率曲线数据和；4)重新选择一个频率段，按照步骤1)至3)分别计算出滤波前和滤波后各个测点两支曲线在该频段范围内视电阻率的算术或者几何平均值和；5)求取各个测点两支曲线的保幅校正因子；6)计算获得最终的保幅校正后两支曲线视电阻率数据和，完成静校正。
【IPC分类】G01V3/38
【公开号】CN104903750
【申请号】CN201480002779
【发明人】王雅苹
【申请人】王雅苹
【公开日】2015年9月9日
【申请日】2014年7月23日
【技术领域】
[0001] 本发明涉及地球物理勘探方法，是一种电磁勘探数据，包括大地电磁（MT)、音频大 地电磁（AMT)和连续剖面法大地电磁（CEMP)的数据处理技术。
【背景技术】
[0002] 电磁法勘探中的静态效应是指由于近表存在局部电性不均匀体时，电阻率分界面 上极化电荷的堆集引起电场的畸变，由此产生一个与外电场成正比的附加电场，且与频率 无关。对于大地电磁数据（MT)，表现在单点曲线上，就是对数域电阻率曲线沿纵轴产生平 移，相位曲线基本上不受影响。表现在视电阻率断面图上，就是电阻率值出现直立的陡变 带，俗称"挂面条"现象，而相应的相位拟断面图上就没有这种现象。表现在单频点的视电 阻率平面图上，就是存在很多"畸变点"现象。静态效应的强度可达两三个数量级，在推断 深度时可引起很大误差，并使构造解释复杂化，因此在反演解释之前需进行静校正。
[0003] 目前，静校正的方法很多，为了压制静态效应，不断提高MT的应用效果，国内外很 多专家发表了相关论文，提出了滤波法、理论计算法、瞬变电磁校正法、地表电阻率校正法、 高频电磁函数向上延拓法等。这些方法大部分都是基于二维MT静校正，各种校正方法的原 理和对资料的要求不同，在实际应用中都受到一定程度的限制。由于滤波法易于实现，处理 速度快，在MT/CEMP资料的静校正处理中应用最为广泛。但是经过滤波法处理后，很多测点 的视电阻率曲线形态会被改变，尤其是在山前带、地表岩性突变带等地方表现明显，实践中 静校正效果不好。
【发明内容】
[0004] 本发明针对采用传统的滤波法进行电磁数据静校正后效果不好的缺陷，提供一种 可以克服视电阻率曲线形态会被改变的缺陷、达到保幅的目的电磁数据的保幅静校正方 法。
[0005] 本发明采用以下步骤实现：
[0006] 1)根据工区采集的大地电磁数据，确定一个深度和电阻率稳定的电性层，得到电 磁数据曲线上对应的频段，对频段上电磁数据测线上所有测点的两支视电阻率曲线做平 均，计算平均视电阻率数据<;
[0007] 所述的电磁数据包括大地电磁数据（MT)、音频大地电磁数据（AMT)、连续剖面法 大地电磁数据（CEMP)。
[0008] 所述的平均包括算术平均和几何平均。
[0009] 所述的计算平均视电阻率数据￥是：
[0010]
[0011] 其中，<(/i)、分别为第i个测点第j个频率A的XY和YX模式实测视 电阻率值，n和m为选定的频点号。
[0012] 2)利用以下公式求取各个测点两支曲线的初步校正因子；
[0014] 其中，f和分别为第i个测点XY和YX模式数据的初步校正因子，Au为由（4) 式获得的计算平均值，其它参数的含义与（4)式相同。
[0015] 3)按照以下公式进行初步校正，得到校正后的两支视电阻率曲线数据/^f和
[0016] 再对Af/f分别进行滤波，获得滤波后的视电阻率曲线数据pf和/其它参 数的含义与（4)、（5)式相同。
[0017] 所述的滤波包括二维中值滤波、二维空间滤波、三维中值滤波和三维空间滤波。
[0018] 4)重新选择一个频率段，按照步骤1)至3)分别计算出滤波前和滤波后各个测点 两支曲线在该频段范围内视电阻率的算术或者几何平均值PfT，和
[0021] 其中，nfl和nf2分别为选定的频点号，Af (乃)分别为第si个测点第j 个频率&的XY和YX模式实测视电阻率值，PfX分别为计算的滤波前XY和YX模式 的视电阻率算术或者几何平均值，和P；f分别为计算的滤波后XY和YX模式的视电阻 率算术或者几何平均值。
[0022] 5)按照以下公式求取各个测点两支曲线的保幅校正因子：
[0024] 其中，#和〇别为第i个测点的XY和YX曲线的保幅校正因子。
[0025] 6)按照以下公式计算获得最终的保幅校正后两支曲线视电阻率数据和 完成静校正。
[0027] 其中，和就是最终的XY和YX曲线的保幅校正后视电阻率数据。
[0028] 与【背景技术】相比，本发明可以克服应用传统的滤波方法后视电阻率曲线形态会被 改变的缺陷，尤其是在山前带、地表岩性突变带等地方，效果改善明显，达到保幅的目的。
【附图说明】
[0029] 图1为电磁数据保幅静校正方法基本流程图。
【具体实施方式】
[0030] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0031] 以下结合附图和实验实例详细说明本发明。
[0032] 图1是保幅静校正方法实施的流程图。具体方法是：
[0033] 1)根据工区采集的大地电磁数据，确定一个深度和电阻率稳定的电性层，得到电 磁数据曲线上对应的频段，如选取第5号到第9号频率段，则公式（4)中n= 5，m= 9。对 频段上电磁数据测线上所有测点的两支视电阻率曲线做平均，计算平均视电阻率数据;
[0034] 所述的电磁数据包括大地电磁数据（MT)、音频大地电磁数据（AMT)、连续剖面法 大地电磁数据（CEMP)。
[0035] 所述的平均包括算术平均和几何平均。
[0036] 所述的计算平均视电阻率数据Pf是：
[0038] 其中，W(/,)、/f(/,)分别为第i个测点第j个频率fj的XY和YX模式实测视 电阻率值，n和m为选定的频点号，在本例子中分别为5和9。
[0039] 2)利用以下公式求取各个测点两支曲线的初步校正因子；
[0041] 其中，f和f分别为第i个测点XY和YX模式数据的初步校正因子，pf为由（4) 式获得的计算平均值，其它参数的含义与（4)式相同。
[0042] 3)按照以下公式进行初步校正，得到校正后的两支视电阻率曲线数据pf和Pf:
[0044] 再对Pf分别进行滤波，如在本例子中，采用二维中值滤波，中值滤波窗口选 择为（5, 3)，获得滤波后的视电阻率曲线数据pjf和其它参数的含义与（4)、（5)式相 同。
[0045] 所述的滤波包括二维中值滤波、二维空间滤波、三维中值滤波和三维空间滤波。
[0046] 4)重新选择一个频率段，如第7号频率到第13号频率段，即公式（7)和公式⑶ 中的nfl= 7，nf2 = 13,按照步骤1)至3)分别计算出滤波前和滤波后各个测点两支曲线 在该频段范围内视电阻率的算术或者几何平均值:
[0049] 其中，nfl和nf2分别为选定的频点号，在本例子中分别为7和13， Af(/,_)、Af(/;)分别为第si个测点第j个频率&的XY和YX模式实测视电阻率值， 、y<}T分别为计算的滤波前XY和YX模式的视电阻率算术或者几何平均值，和 分别为计算的滤波后XY和YX模式的视电阻率算术或者几何平均值。
[0050] 5)按照以公式（9)求取各个测点两支曲线的保幅校正因子：
[0052] 其中，句厂和片分别为第i个测点的XY和YX曲线的保幅校正因子。
[0053] 6)按照公式（10)计算获得最终的保幅校正后两支曲线视电阻率数据和， 完成保幅静校正。
[0055] 其中，和就是最终的XY和YX曲线的保幅校正后视电阻率数据。
[0056] 以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详 细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保 护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本 发明的保护范围之内。
【主权项】
1. 一种电磁数据的保幅静校正方法，其特征在于采用以下步骤实现： 1) 根据工区采集的大地电磁数据，确定一个深度和电阻率稳定的电性层，得到电磁数 据曲线上对应的频段，对频段上电磁数据测线上所有测点的两支视电阻率曲线做平均，计 算平均视电阻率数据<; 所述的计算平均视电阻率数据P是：其中，αλτ(Λ)、AnC/;)分别为第i个测点第j个频率A的XY和YX模式实测视电阻 率值，η和m为选定的频点号； 2) 利用以下公式求取各个测点两支曲线的初步校正因子；其中，和Af分别为第i个测点XY和YX模式数据的初步校正因子，Pf为由（4)式 获得的计算平均值，其它参数的含义与（4)式相同。 3) 按照以下公式进行初步校正，得到校正后的两支视电阻率曲线数据/jf和Pf :(6) 再对Af 分别进行滤波，获得滤波后的视电阻率曲线数据/7=和其它参数的 含义与（4)、（5)式相同。 4) 重新选择一个频率段，按照步骤1)至3)分别计算出滤波前和滤波后各个测点两支 曲线在该频段范围内视电阻率的算术或者几何平均值:(7) (8) 其中，nfl和nf2分别为选定的频点号，Af Pf (/v)分别为第si个测点第j个 频率&的XY和YX模式实测视电阻率值，Pfy、分别为计算的滤波前XY和YX模式的 视电阻率算术或者几何平均值，和/分别为计算的滤波后XY和YX模式的视电阻率 算术或者几何平均值。 5) 按照以下公式求取各个测点两支曲线的保幅校正因子：其中，和#分别为第i个测点的XY和YX曲线的保幅校正因子。 YY 6) 按照以下公式计算获得最终的保幅校正后两支曲线视电阻率数据/和户#，完成 静校正；(10) 其中，/^和就是最终的XY和YX曲线的保幅校正后视电阻率数据。2. 根据权利要求1所述的一种电磁数据的保幅静校正方法，其特征在于步骤1)所述 的电磁数据包括大地电磁数据（MT)、音频大地电磁数据（AMT)、连续剖面法大地电磁数据 (CEMP)〇3. 根据权利要求1所述的一种电磁数据的保幅静校正方法，其特征在于步骤1)所述的 平均包括算术平均和几何平均。4. 根据权利要求1所述的一种电磁数据的保幅静校正方法，其特征在于步骤3)所述的 滤波包括二维中值滤波、二维空间滤波、三维中值滤波和三维空间滤波，其中：μ是土体介质的泊松比。
【专利摘要】一种电磁数据的保幅静校正方法，采用以下步骤实现：1)根据工区采集的大地电磁数据，确定一个深度和电阻率稳定的电性层，得到电磁数据曲线上对应的频段，对频段上电磁数据测线上所有测点的两支视电阻率曲线做平均，计算平均视电阻率数据；2)求取各个测点两支曲线的初步校正因子；3)进行初步校正，得到校正后的两支视电阻率曲线数据和；4)重新选择一个频率段，按照步骤1)至3)分别计算出滤波前和滤波后各个测点两支曲线在该频段范围内视电阻率的算术或者几何平均值和；5)求取各个测点两支曲线的保幅校正因子；6)计算获得最终的保幅校正后两支曲线视电阻率数据和，完成静校正。
【IPC分类】G01V3/38
【公开号】CN104903750
【申请号】CN201480002779
【发明人】王雅苹
【申请人】王雅苹
【公开日】2015年9月9日
【申请日】2014年7月23日
最新回复(0)