山野莽夫

  • 归档
    • 随笔
    • 建站资源
    • 分享
    • 代码
  • 地球物理学
    • 专业课
    • 概念解释
  • 计算机
  • 互联网
  • 教程
  • 规划
  • 实验室
    • 珍藏的软件
    • 贴吧云签到
    • A1账号自助申请
山野莽夫
小学生的挣扎的点点滴滴
  1. 首页
  2. 地球物理学
  3. 专业课
  4. 地磁学
  5. 正文

试验项目二 二度体磁异常正演

2016年1月28日 3640点热度 0人点赞 0条评论

一、试验要求
根据水平圆柱体和有限延伸厚板正演公式,计算给定参数二度体磁场X分量,Z分量,总场异常,并总结磁异常特征
二、实验课时:4课时
三、实验成果提交 :各磁场分量及总场主剖面图
四、水平圆柱体及斜交磁化厚板正演公式

 



 
五、正演参数

水平圆

柱体

埋深

计算范围(X,Y)

有效磁矩ms

磁倾角(I)

磁偏角(A)

300

(-500,500)

314

45

45

无限延伸厚板

埋深

计算范围(X,Y)

板宽

(2b)

有效磁化强度

板体倾角

磁倾角(I)

磁偏角

300

(-500,

500)

20

1

45

45

45

 

  1. 水平圆柱体磁异常图

(1)水平圆柱体水平X分量图(倾角0-90变化规律)


(2)水平圆柱体垂直Z分量图(倾角0-90变化规律)


 

 

(3)总磁异常图(倾角0-90变化规律)


磁异常特征:当垂直磁化时Za为两边有负值的对称曲线,而Ha为原点对称曲线;若为水平磁化,可以知道受磁化倾角的影响比Za大,若水平圆柱体斜磁化,则Ha,Za,均为非对称曲线一般正值小于Za,而负值大于Za。水平磁化的Z分量曲线和垂直磁化时的X分量曲线重合。

  1. 厚板磁异常图

(1)、厚板水平X分量(角0-90变化规律)


(2)、厚板垂直Z分量(角0-90变化规律)


(3)厚板总磁异常图(角0-90变化规律)


磁异常特征分析:顺层磁化时,Za剖曲线对称于Z轴,且没有负值,Hax曲线对称于原点。相当于Za和Ha分别和不同角函数的合成。当角从零度变到九十度时,Za曲线有轴对称变成原点对称,Hax则由原点对称变为反轴向对称。Za为非对称曲线,其不对称性随角增大而增大。

七、源代码

1、水平圆柱体

#define
_CRT_SECURE_NO_WARNINGS

#include
"stdafx.h"

#include<stdio.h>

#include<math.h>

#include<stdlib.h>

#define
dcx1
"file1.dat"

#define
dcx2
"file2.dat"

#define
dcx3
"file3.dat"

int main()

{

    

    double Za, Ha, T, M = 1, is = 0.785,I=0.785,r=10, R = 300,A=0.785,is0,x; int i;

    FILE *fp1,*fp2,*fp3;

    

    is0 = atan(1.0 / cos(A) / tan(I));

    fopen_s(&fp1,dcx1, "a+");

    fopen_s(&fp2, dcx2, "a+");

    fopen_s(&fp3, dcx3, "a+");

    

    for (x = -500; x <= 500; x = x + 10)

    {

        fprintf(fp1, "%f\t", x);

        fprintf(fp2, "%f\t", x);

        fprintf(fp3, "%f\t", x);

 

        for (is = 0; is <= 90; is = is + 30)

        {

        

            Za = 1256 * r * r / 2.0 / pow((x*x + R*R), 2.0)*((R*R - x*x)*sin(is / 180.0 * 3.14) - 2 * R*x*cos(is / 180.0 * 3.14));

 

            Ha = -1256 * r * r / 2.0 / pow((x*x + R*R), 2.0)*((R*R - x*x)*cos(is / 180.0 * 3.14) + 2 * R*x*sin(is / 180.0 * 3.14));

            

            T = -1256 * r * r / 2.0 / pow((x*x + R*R), 2.0)*sin(I) / sin(is0)*((R*R - x*x)*cos(is / 180.0 * 3.14 + is0) + 2 * R*x*

                

                sin(is / 180.0 * 3.14));

            fprintf(fp1, "%f\t", Za);

            fprintf(fp2, "%f\t", Ha);

            fprintf(fp3, "%f\t", T);

            //printf("%f\t%f\t%f\t%f\n", x, Za, Ha, T);

 

            //fprintf(fp1, "%f\t%f\t%f\t%f\n", x, Za, Ha, T);

            

        }

        fprintf(fp1, "\n");

        fprintf(fp2, "\n");

        fprintf(fp3, "\n");

    }

    fclose(fp1);

    fclose(fp2);

    fclose(fp3);

    return 0;

}

 

2、厚板

#define
_CRT_SECURE_NO_WARNINGS

#include
"stdafx.h"

#include<stdio.h>

#include<stdlib.h>

#include<math.h>

#define
dcx1
"file1.dat"

#define
dcx2
"file2.dat"

#define
dcx3
"file3.dat"

int main( )

{

    float u0 = 1256, Ms = 1,a=45,x,r,b=10,h=300,Z,H,T,I=0.785,is;

    FILE *fp1, *fp2, *fp3;

 

    fopen_s(&fp1, dcx1, "a+");

    fopen_s(&fp2, dcx2, "a+");

    fopen_s(&fp3, dcx3, "a+");//非标准C库函数

    
 

        

        for (x = -500; x <= 500; x = x + 10)

        {

            fprintf(fp1, "%f\t", x);

            fprintf(fp2, "%f\t", x);

            fprintf(fp3, "%f\t", x);

            

            for (r = 0; r <= 90; r = r + 30)

            {

 

 

                is = a - r;

            Z = u0*Ms*sin(a / 180.0*3.14) / 2.0 / 3.14*(0.5*sin(r / 180.0*3.14)*log(((x - b)*(x - b) + h*h) / ((x + b)*(x + b) + h*h)) + cos(r / 180.0*3.14)*(atan((x + b) / h) - atan((x - b) / h)));

 

            H = u0*Ms*sin(a / 180.0*3.14) / 2.0 / 3.14*(0.5*cos(r / 180.0*3.14)*log(((x - b)*(x - b) + h*h) / ((x + b)*(x + b) + h*h)) - sin(r / 180.0*3.14)*(atan((x + b) / h) - atan((x - b) / h)));

 

            T = u0*Ms*sin(a / 180.0*3.14) / 2.0 / 3.14*sin(I) / sin(is / 180.0*3.14)*(0.5*cos(a / 180.0*3.14 - 2 * is / 180.0*3.14)*log(((x - b)*(x - b) + h*h) / ((x + b)*(x + b) + h*h)) - sin(1 - 2 * is / 180.0*3.14)*

                

                (atan((x + b) / h) - atan((x - b) / h)));//计算

            

            fprintf(fp1, "%f\t", Z);

            fprintf(fp2, "%f\t", H);

            fprintf(fp3, "%f\t", T);

             //printf("%f\t%f\t%f\t%f\n", x, Z, H, T);

 

            //fprintf(fp, "%f\t%f\t%f\t%f\n", x, Z, H, T);//输出到文件

        }

 

        fprintf(fp1, "\n");

        fprintf(fp2, "\n");

        fprintf(fp3, "\n");

    }

    fclose(fp1);

    fclose(fp2);

    fclose(fp3);

    return 0;

}

 

 

标签: 二度体 地磁学 地磁学实验报告 地磁学程序 地磁学试验 水平圆柱体 水平无限厚板 磁异常
最后更新:2016年1月28日

小菜菜

菜鸟

打赏 点赞
< 上一篇
下一篇 >

文章评论

razz evil exclaim smile redface biggrin eek confused idea lol mad twisted rolleyes wink cool arrow neutral cry mrgreen drooling persevering
取消回复

此站点使用Akismet来减少垃圾评论。了解我们如何处理您的评论数据。

联系方式

QQ群 | TG群 | 邮箱

最新 热点 随机
最新 热点 随机
Azure Student 微软云 学生订阅 免费12个月用量避坑注意点集合 MP3音频文件格式详细解析 python按固定采样点个数分割wav格式音频 愉快使用谷歌免费人工智能平台colab,训练你的神经网络模型,为你的学术生活添砖加瓦 华为云版轻量应用服务器-云耀云服务器简单体验评测 Cloudflare 免费CDN自定义节点ip之自选cloudflare 高速节点ip工具分享
本站提供百度贴吧自动签到工具 腾讯云服务器安装lnmp一键包 LOLHunter - 亨特儿 最好用的英雄联盟免费刷金币脚本 二层电测深电阻率正演模拟 谷歌云平台gcp300美元试用失败经历。 素数算法-试除法和筛选法的不同境界与C语言实现
标签聚合
wordpress 地震学程序 虚拟机 onedrive ppt 宝塔面板 模板 c语言
最近评论
小菜菜 发布于 6 个月前(11月24日) 这玩意已经废了,成收割工具了,不能再用了。
eamon 发布于 7 个月前(11月07日) 我一年不用了才发现这个休眠管理费每月15,一共扣了我135元,然后我消费还消费不了,我宁愿消费掉也不...
magic 发布于 11 个月前(07月03日) 请问账号不注销会有什么影响吗?
magic 发布于 11 个月前(07月01日) 我想问一下 如果不注销账号就留着会怎么样
qwp6601 发布于 12 个月前(06月04日) 有没有方法改为bing

COPYRIGHT © 2021 shanyemangfu.com. ALL RIGHTS RESERVED.

Theme Kratos Made By Seaton Jiang

蜀ICP备15031791号-2