山野莽夫

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

三度体球体磁场正演

2016年1月28日 3006点热度 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日

小菜菜

菜鸟

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

文章评论

取消回复

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

标签聚合
wordpress 地震学程序 虚拟机 宝塔面板 c语言 模板 onedrive ppt
最新 热点 随机
最新 热点 随机
Azure Student 微软云 学生订阅 免费12个月用量避坑注意点集合 MP3音频文件格式详细解析 python按固定采样点个数分割wav格式音频 愉快使用谷歌免费人工智能平台colab,训练你的神经网络模型,为你的学术生活添砖加瓦 华为云版轻量应用服务器-云耀云服务器简单体验评测 Cloudflare 免费CDN自定义节点ip之自选cloudflare 高速节点ip工具分享
c语言中函数二维数组作为参数传递方法 层速度,平均速度,叠加速度,均方根速度的关系 多线程bt、http、ftp下载工具Aria2开启https模式以及面板的配置 谷歌云服务Google Cloud(GCP)不换号重新开启300美元试用并且保留老号虚拟机实例的方法 nginx子目录绑定二级域名的方法 吐槽电信校园宽带的坑

COPYRIGHT © 2021 shanyemangfu.com. ALL RIGHTS RESERVED.

Theme Kratos Made By Seaton Jiang

蜀ICP备15031791号-2