山野莽夫

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

三度体球体磁场正演

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

三度体球体磁场正演

一 、 试验要求:

根据正演公式,计算给定参数球体磁场水平X分量、水平Y分量,垂直z分量,总磁异常,并总结异常特征。

二 、实验课时:2课时

三 、实验成果提交:

各磁场分量及总场平面等值线图、主剖面图。

四 、球体正演公式

五 、球体正演参数

埋深(H)

半径(R)

计算范围(x,y)

磁化强度(M)

磁倾角(I)

磁偏角(A)

200

100

(-500,500)

1

45

45

六 、球体正演图(平面等值线图和主剖面图)

  1. 水平X分量图
    1. 平面等值线图

  1. 主剖面

  1. 水平Y分量
    1. 平面等值线图

    1. 主剖面

  1. 垂直Z分量图
    1. 平面等值线图
    2. 主剖面

  1. 总场异常图
    1. 平面等值线图

  1. 主剖面图

异常特征分析:Hax,Hay,Za,主剖面过正负异常中心。 Za等值线呈现等轴状,负异常包着正异常:极大值和极小值的连线对应磁化强度M在平面上的投影方向:极小值位于正异常的北侧,极大值位于坐标原点的南侧,用Za平面等值线图中的极大值和极小值的连线确定主剖面。由图可知,受磁化倾角的影响比Za大,在相同的磁化倾角其负值较大。

七、源代码

#include<stdio.h>

#include<math.h>

#include<stdlib.h>

#define dz1 "file1.dat"

main()

{

double m,M=1;int x,y,z=0,R=200,R0=100;float I=0.785,A=0.785,Ha,Hb,Za,T;int i,j;double s1,s2,s3;

FILE *fp;

fp=fopen(dz1,"w+");

fclose(fp);

fp=fopen(dz1,"a+");//文件操作

m=3.14*R0*R0*R0*4/3;

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

{

for(y=-500;y<=500;y=y+10)//计算开始

{

Ha=100*m/(pow((x*x+y*y+R*R),(2.5)))*((2*x*x-y*y-R*R)*cos(I)*cos(A)-3*R*x*sin(I)+3*x*y*cos(I)*sin(A));

Hb=100*m/(pow((x*x+y*y+R*R),(2.5)))*((2*y*y-x*x-R*R)*cos(I)*sin(A)-3*R*y*sin(I)+3*x*y*cos(I)*cos(A));

Za=100*m/(pow((x*x+y*y+R*R),(2.5)))*((2*R*R-y*y-x*x)*sin(I)-3*R*x*cos(I)*cos(A)-3*R*y*cos(I)*sin(A));

T =100*m/(pow((x*x+y*y+R*R),(2.5)))*((2*R*R-y*y-x*x)*sin(I)*sin(I)+(2*x*x-y*y-R*R)*cos(I)*cos(I)*cos(A)*

cos(A)+(2*y*y-x*x-R*R)*cos(I)*cos(I)*sin(A)*sin(A)-3*R*x*sin(2*I)*cos(A)+3*x*y*cos(I)*cos(I)*cos(I)*sin(2*A)

-3*R*y*sin(2*I)*sin(A) );//计算

printf("%f\t%f\t%f\t%f\n",Ha,Hb,Za,T);//输出到屏幕

fprintf(fp,"%d\t%d\t%f\t%f\t%f\t%f\n",x,y,Ha,Hb,Za,T);//输出到文件

}

}

fclose(fp);

}

 

 

标签: 三度体 地磁学 地磁学实验报告 地磁学程序 正演 磁异常
最后更新: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工具分享
为自己的文件多准备几个坑-借助使用免费Mover.io文件备份同步迁移工具安心的使用office 365 5T网盘(支持多种存储工具如onedrive和谷歌网盘) 给wordpress博客添加博客统计小工具 太子妃升职记360云盘资源 excel工作表(表格xlsx)去除工作表视图保护和工作簿保护 请注意lycamobile已经凉了 斯奈尔定律和Zoeppritz方程
标签聚合
模板 ppt c语言 地震学程序 wordpress 宝塔面板 虚拟机 onedrive
最近评论
小菜菜 发布于 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