山野莽夫

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

三度体球体磁场正演

2016年1月28日 3535点热度 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来减少垃圾评论。了解我们如何处理您的评论数据。

标签聚合
c语言 模板 ppt 地震学程序 虚拟机 wordpress onedrive 宝塔面板
最新 热点 随机
最新 热点 随机
Azure Student 微软云 学生订阅 免费12个月用量避坑注意点集合 MP3音频文件格式详细解析 python按固定采样点个数分割wav格式音频 愉快使用谷歌免费人工智能平台colab,训练你的神经网络模型,为你的学术生活添砖加瓦 华为云版轻量应用服务器-云耀云服务器简单体验评测 Cloudflare 免费CDN自定义节点ip之自选cloudflare 高速节点ip工具分享
Vmware 虚拟机安装windows 10 没有网络的解决办法 指环王1、2、3下载链接 wordpress给评论者链接添加nofollow、新窗口打开以及go跳转 地震学和地球结构(震源追踪) 关于 Winrar 压缩软件 去弹窗广告注册纯净版含注册key文件

COPYRIGHT © 2021 shanyemangfu.com. ALL RIGHTS RESERVED.

Theme Kratos Made By Seaton Jiang

蜀ICP备15031791号-2