一种电磁数据的保幅静校正方法

摘要

本发明是电磁勘探数据处理技术，采集电磁数据，确定一个深度和电阻率稳定的电性层，对对应频段上电磁数据测线上所有测点的两支视电阻率曲线做平均，计算平均视电阻率数据；求取各个测点两支曲线的初步校正因子校正得到两支视电阻率曲线数据，重新选择频率段分别计算出滤波前和滤波后各个测点两支曲线在该频段范围内视电阻率的平均值，求取各测点两支曲线的保幅校正因子，计算获得最终的保幅校正后两支曲线视电阻率数据。本发明可以克服应用传统的滤波方法后视电阻率曲线形态会被改变的缺陷，尤其是在山前带、地表岩性突变带等地方，效果改善明显，达到保幅的目的。

主权项

一种电磁数据的保幅静校正方法，特点是采用以下步骤实现：1）根据工区采集的大地电磁数据，确定一个深度和电阻率稳定的电性层，得到电磁数据曲线上对应的频段，对频段上电磁数据测线上所有测点的两支视电阻率曲线做平均，计算平均视电阻率数据2）利用以下公式求取各个测点两支曲线的初步校正因子；$k_i^{\mathrm{XY}}=\log {\rho }_i^a-\log \frac{\underset{j=n}{\overset{m}{\Sigma }}{\rho }_i^{\mathrm{XY}}\left(f_j\right)}{m-n+1},$$k_i^{\mathrm{YX}}=\log {\rho }_i^a-\log \frac{\underset{j=n}{\overset{m}{\Sigma }}{\rho }_i^{\mathrm{YX}}\left(f_j\right)}{m-n+1}---\left(5\right)$其中，和分别为第i个测点XY和YX模式数据的初步校正因子，为由(4)式获得的计算平均值，其它参数的含义与(4)式相同；3）按照以下公式进行初步校正，得到校正后的两支视电阻率曲线数据和${\rho }_{\mathrm{si}}^{\mathrm{XY}}=10**(\log {\rho }_i^{\mathrm{XY}}+k_i^{\mathrm{XY}}),$${\rho }_{\mathrm{si}}^{\mathrm{YX}}=10**(\log {\rho }_i^{\mathrm{YX}}+k_i^{\mathrm{YX}})---\left(6\right)$再对分别进行滤波，获得滤波后的视电阻率曲线数据和其它参数的含义与(4)、(5)式相同；4）重新选择一个频率段，按照步骤1）至3）分别计算出滤波前和滤波后各个测点两支曲线在该频段范围内视电阻率的算术或者几何平均值和${\rho }_{\mathrm{si}}^{\mathrm{aXY}}=\frac{\underset{j=\mathrm{nf}1}{\overset{\mathrm{nf}2}{\Sigma }}{\rho }_{\mathrm{si}}^{\mathrm{XY}}\left(f_j\right)}{\mathrm{nf}2-\mathrm{nf}1+1},$${\rho }_{\mathrm{si}}^{\mathrm{aYX}}=\frac{\underset{j=\mathrm{nf}1}{\overset{\mathrm{nf}2}{\Sigma }}{\rho }_{\mathrm{si}}^{\mathrm{YX}}\left(f_j\right)}{\mathrm{nf}2-\mathrm{nf}1+1}---\left(7\right)$${\rho }_{\mathrm{fsi}}^{\mathrm{aXY}}=\frac{\underset{j=\mathrm{nf}1}{\overset{\mathrm{nf}2}{\Sigma }}{\rho }_{\mathrm{fsi}}^{\mathrm{XY}}\left(f_j\right)}{\mathrm{nf}2-\mathrm{nf}1+1},$${\rho }_{\mathrm{fsi}}^{\mathrm{aYX}}=\frac{\underset{j=\mathrm{nf}1}{\overset{\mathrm{nf}2}{\Sigma }}{\rho }_{\mathrm{fsi}}^{\mathrm{YX}}\left(f_j\right)}{\mathrm{nf}2-\mathrm{nf}1+1}---\left(8\right)$其中，nf1和nf2分别为选定的频点号，分别为第si个测点第j个频率f_j的XY和YX模式实测视电阻率值，分别为计算的滤波前XY和YX模式的视电阻率算术或者几何平均值，和分别为计算的滤波后XY和YX模式的视电阻率算术或者几何平均值；5）按照以下公式求取各个测点两支曲线的保幅校正因子：$k_{\mathrm{fi}}^{\mathrm{XY}}=\log {\rho }_{\mathrm{si}}^{\mathrm{aXY}}-\log {\rho }_{\mathrm{fsi}}^{\mathrm{aXY}},$$k_{\mathrm{fi}}^{\mathrm{YX}}=\log {\rho }_i^{\mathrm{aYX}}-\log {\rho }_{\mathrm{fsi}}^{\mathrm{aYX}}---\left(9\right)$其中，和分别为第i个测点的XY和YX曲线的保幅校正因子；6）按照以下公式计算获得最终的保幅校正后两支曲线视电阻率数据和完成静校正；${\rho }_{\mathrm{psi}}^{\mathrm{XY}}=10**(\log {\rho }_{\mathrm{si}}^{\mathrm{XY}}+k_{\mathrm{fi}}^{\mathrm{XY}}),$${\rho }_{\mathrm{psi}}^{\mathrm{YX}}=10**(\log {\rho }_{\mathrm{si}}^{\mathrm{YX}}+k_{\mathrm{fi}}^{\mathrm{YX}})---\left(10\right)$其中，和就是最终的XY和YX曲线的保幅校正后视电阻率数据。