一种快速反演地下地质体空间位置的方法与流程

文档序号:12359969阅读:来源:国知局

技术特征:

1.一种快速反演地下地质体空间位置的方法,依次包括以下步骤:

步骤一、在勘查区开展高精度磁法测量,获得磁测数据;

步骤二、对磁测数据进行日变校正、正常场梯度改正、高度改正及基点改正,得到测点异常△T值;

步骤三、将步骤二所得到的磁异常△T值进行数据平滑处理,消除高频干扰,得到光滑的磁异常曲线并对所有光滑处理后的磁异常△T值进行网格化处理;

步骤四、对步骤三中的处理后的磁异常△T值进行化极处理得到化极磁异常△T值网格,对化极磁异常△T值网格进行X、Y、Z三个方向求取导数并得到X、Y、Z三个方向导数网格;

步骤五、步骤四中化极磁异常△T值网格和步骤四中X、Y、Z三方向导数网格进行分割,把以上图件剖分成小网格;

步骤六、对每个小网格构造指数进行赋值,并且利用欧拉三维反褶积反演计算,得出所有小网格内的欧拉解集,设置解集误差容许范围指数为15%并筛选出符合条件的小网格内的欧拉解;

步骤七、提取出小网格欧拉解集,把所有欧拉解的坐标、深度数据编号并放在同一数据库;

步骤八、对所有解集数据主要是解的坐标和深度进行成图处理,对步骤六所述小网格求出的欧拉解的小地质单元汇聚起来并且在三维空间展示则可形成模拟的地质体,最终形成欧拉解集三维立体图,返回步骤六对小网格单元构造指数进行重新赋值,重复步骤七、步骤八,得出最后解译成果图件。

2.根据权利要求1所述的一种快速反演地下地质体空间位置的方法,其特征在于:所述的步骤一中磁测总精度应高于2nT。

3.根据权利要求2所述的一种快速反演地下地质体空间位置的方法,其特 征在于:所述的步骤五中,所述的小网格个数为n个,大小为10m×10m。

4.根据权利要求3所述的一种快速反演地下地质体空间位置的方法,其特征在于:所述的步骤六中,对每个小网格构造指数进行赋值,对为岩体接触带的网格赋值为0,对岩墙或岩席的网格赋值为1,对水平管状体的网格赋值2,对近似球体形状的岩体的网格赋值为3。

5.根据权利要求4所述的一种快速反演地下地质体空间位置的方法,其特征在于:欧拉搜索窗口大小为5m×5m。

6.根据权利要求5所述的一种快速反演地下地质体空间位置的方法,其特征在于:所述的欧拉三维反褶积运算采用下式进行计算:

<mrow> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <msub> <mi>x</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mfrac> <mrow> <mo>&PartialD;</mo> <mi>f</mi> </mrow> <mrow> <mo>&PartialD;</mo> <mi>x</mi> </mrow> </mfrac> <mo>+</mo> <mrow> <mo>(</mo> <mi>y</mi> <mo>-</mo> <msub> <mi>y</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mfrac> <mrow> <mo>&PartialD;</mo> <mi>f</mi> </mrow> <mrow> <mo>&PartialD;</mo> <mi>f</mi> </mrow> </mfrac> <mo>+</mo> <mrow> <mo>(</mo> <mi>z</mi> <mo>-</mo> <msub> <mi>z</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mfrac> <mrow> <mo>&PartialD;</mo> <mi>f</mi> </mrow> <mrow> <mo>&PartialD;</mo> <mi>z</mi> </mrow> </mfrac> <mo>=</mo> <mi>N</mi> <mrow> <mo>(</mo> <mi>B</mi> <mo>-</mo> <mi>f</mi> <mo>)</mo> </mrow> </mrow>

式中N为构造指数值,B为背景磁场值,f为地下地质体场源(x0,y0,z0)在观测点(x,y,z)处的磁场函数表达式。

当前第2页1 2 3 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1